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Preface 

The present notes cover a series of three lectures on the quantum Hall effect 
given at the Singapore session "Ultracold Gases and Quantum Information" at 
Les Houches Summer School 2009. Almost 30 years after the discovery of the 
quantum Hall effect, the research subject of quantum Hall physics has definitely 
acquired a high degree of maturity that is reflected by a certain number of 
excellent reviews and books, of which we can cite only a few [HOG] for possible 
further or complementary reading. Also the different sessions of Les Houches 
Summer School have covered in several aspects quantum Hall physics, and S. 
M. Girvin's series of lectures in 1998 [3] have certainly become a reference in 
the fieldH Girvin's lecture notes were indeed extremely useful for myself when 
I started to study the quantum Hall effect at the beginning of my Master and 
PhD studies. 

The present lecture notes are complementary to the existing literature in 
several aspects. One should first mention its introductory character to the field, 
which is in no way exhaustive. As a consequence, the presentation of one-particle 
physics and a detailed discussion of the integer quantum Hall effect occupy the 
major part of these lecture notes, whereas the - certainly more interesting 
fractional quantum Hall effect, with its relation to strongly-correlated electrons, 
its fractionally charged quasi-particles and fractional statistics, is only briefly 
introduced. 

Furthermore, we have tried to avoid as much as possible the formal aspects 
of the fractional quantum Hall effect, which is discussed only in the framework of 
trial wave functions a la Laughlin. We have thus omitted, e.g., a presentation of 
Chern-Simons theories and related quantum-field theoretical approaches, such 
as the Hamiltonian theory of the fractional quantum Hall effect [5], as much 
as the relation between the quantum Hall effect and conformal field theories. 
Although these theories are extremely fruitful and still promising for a deeper 
understanding of quantum Hall physics, a detailed discussion of them would 
require more space than these lecture notes with their introductory character 
can provide. 

Another complementary aspect of the present lecture notes as compared 
to existing textbooks consists of an introduction to Landau-level quantisation 
that treats in a parallel manner the usual non-relativistic electrons in semicon- 
ductor heterostructures and relativistic electrons in graphene (two-dimensional 
graphite). Indeed, the 2005 discovery of a quantum Hall effect in this amazing 
material [5J [7] has given a novel and unexpected boost to research in quantum 
Hall physics. 

As compared to the (oral) lectures, the present notes contain slightly more 
information. An example is Laughlin's plasma analogy, which is described in 
Sec. 14.2.51 although it was not discussed in the oral lectures. Furthermore, 
I have decided to add a chapter on multi-component quantum Hall systems, 
which, for completeness, needed to be at least briefly discussed. 

1 These lectures are also available on the preprint server, 

http: / / arxiv.org /abs /cond-mat 79907002] 
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Before the Singapore session of Les Houches Summer School, this series of 
lectures had been presented in a similar format at the (French) Summer School of 
the Research Grouping "Physique Mesoscopique" at the Institute of Scientific 
Research, Cargese, Corsica, in 2008. Furthermore, a longer series of lectures 
on the quantum Hall effect was prepared in collaboration with my colleague 
and former PhD advisor Pascal Lederer (Orsay, 2006). Its aim was somewhat 
different, with an introduction to the Hamiltonian theories of the fractional 
quantum Hall effect and correlation effects in multi-component systems. As 
already mentioned above, the latter aspect is only briefly introduced within the 
present lecture notes and a discussion of Hamiltonian theories is completely 
absent. The Orsay series of lectures was repeated by Pascal Lederer at the 
Ecole Polytechnique Federate in Lausanne Switzerland, in 2006, and at the 
University of Recife, Brazil, in 2007. The finalisation of these longer and more 
detailed lecture notes (in French) is currently in progress. The graphene- related 
aspects of the quantum Hall effect have furthermore been presented in a series 
of lectures on graphene (Orsay, 2008) prepared in collaboration with Jean-Noel 
Fuchs, whom I would like to thank for a careful reading of the present notes. 
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Chapter 1 

Introduction 



Quantum Hall physics - the study of two-dimensional (2D) electrons in a strong 
perpendicular magnetic field [see Fig. Il.lf a)] - has become an extremely im- 
portant research subject during the last two and a half decades. The interest 
for quantum Hall physics stems from its position at the borderline between low- 
dimensional quantum systems and systems with strong electronic correlations, 
probably the major issues of modern condensed-matter physics. From a theoret- 
ical point of view, the study of quantum Hall systems required the elaboration 
of novel concepts some of which were better known in quantum-field theories 
used in high-energy rather than in condensed-matter physics, such e.g. charge 
fractionalisation, non-commutative geometries and topological field theories. 

The motivation of the present lecture notes is to provide in an accessible man- 
ner the basic knowledge of quantum Hall physics and to enable thus interested 
graduate students to pursue on her or his own further studies in this subject. We 
have therefore tried, whereever we feel that a more detailed discussion of some 
aspects in this large field of physics would go beyond the introductory character 
of these notes, to provide references to detailed and pedagogical references or 
complementary textbooks. 

1.1 History of the (Quantum) Hall Effect 
1.1.1 The physical system 

Our main knowledge of quantum Hall systems, i.e. a system of 2D electrons in 
a perpendicular magnetic field, stems from electronic transport measurements, 
where one drives a current / through the sample and where one measures both 
the longitudinal and the transverse resistance (also called Hall resistance). The 
difference between these two resistances is essential and may be defined topo- 
logically: consider a current that is driven through the sample via two arbitrary 
contacts [CI and C4 in Fig. 11.1( a)] and draw (in your mind) a line between 
these two contacts. A longitudinal resistance is a resistance measured between 
two (other) contacts that may be connected by a line that does not cross the 
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resistance resistance magnetic field B 



Figure 1.1: (a) 2D electrons in a perpendicular magnetic field (quantum Hall 
system). In a typical transport measurement, a current / is driven through 
the system via the contacts CI and C4. The longitudinal resistance may be 
measured between the contacts C5 and C6 (or alternatively between C2 and 
C3). The transverse (or Hall) resistance is measured, e.g., between the contacts 
C3 and C5. (b) Classical Hall resistance as a function of the magnetic field. 



line connecting CI and C4. In Fig. Il.lf a). we have chosen the contacts C5 
and C6 for a possible longitudinal resistance measurement. The transverse re- 
sistance is measured between two contacts that are connected by an imaginary 
line that necessarily crosses the line connecting CI and C4 [e.g. C3 and C5 in 
Fig. EKb)]. 



1.1.2 Classical Hall effect 

Evidently, if there is a quantum Hall effect, it is most natural to expect that 
there exists also a classical Hall effect. This is indeed the case, and its history 
goes back to 1879 when Hall showed that the transverse resistance Rh of a thin 
metallic plate varies linearly with the strength B of the perpendicular magnetic 
field [Fig. ITTTbl]. 

B , , 

Rh = , 1.1 

qn e i 

where q is the carrier charge (q = —e for electrons in terms of the elementary 
charge e that we define positive in the remainder of these lectures) and n e i is 
the 2D carrier density. Intuitively, one may understand the effect as due to 
the Lorentz force, which bends the trajectory of a charged particle such that 
a density gradient is built up between the two opposite sample sides that are 
separated by the contacts CI and C4. Notice that the classical Hall resistance is 
still used today to determine, in material science, the carrier charge and density 
of a conducting material. 

More quantitatively, the classical Hall effect may be understood within the 
Drude model for diffusive transport in a metal. Within this model, one considers 
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independent charge carriers of momentum p described by the equation of motion 



dt \ TUb J T 

where E and B are the electric and magnetic fields, respectively. Here, we con- 
sider transport of negatively charged particles (i.e. electrons with q = — e) with 
band mass irib. The last term takes into account relaxation processes due to the 
diffusion of electrons by generic impurities, with a characteristic relaxation time 
t. The macroscopic transport characteristics, i.e. the resistivity or conductivity 
of the system, are obtained from the static solution of the equation of motion, 
dp/dt = 0, and one finds for 2D electrons with p = (p x ,p y ) 

p, cB p x 
eE x = p y , 

Wib T 

tp eB Py 
ehy = — p x , 

m b t 

where we have chosen the magnetic field in the z-dircction. In the above ex- 
pressions, one notices the appcarcncc of a characteristic frequency, 

= — , (1-2) 

m b 

which is called cyclotron frequency because it characterises the cyclotron mo- 
tion of a charged particle in a magnetic field. With the help of the Drudc 
conductivity, 

n eI e 2 T 

oo = , (1.3) 

m b 

one may rewrite the above equations as 

<j E x = -en e i— - en e i^-(uj C T), 
m b m b 

<JoEy = eriei — ivcr) — en e i — , 
m b m b 



or, in terms of the current density 



j = -en ei — , (1.4) 
m b 



in matrix form as E = pj, with the resistivity tensor 

CTO V _W CT 1 J cr V ~L lB 1 

where we have introduced, in the last step, the mobility 

CT 



(1.5) 



M =— . (1.6) 

m b 
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From the above expression, one may immediately read off the Hall resistivity 
(the off-diagonal terms of the resistivity tensor p) 

uj c t eB m b B 
pH = = — T x — = . 1.7 

cr m b n e ie z T en el 



Furthermore, the conductivity tensor is obtained from the resistivity (|1.5j) . by 
matrix inversion, 

' = = <L8) 

with o l = co/(l + loqT 2 ) and oh = <^o u) c T /{^ + uj c t2 )- ^ ^ s instructive to 
discuss, based on these expressions, the theoretical limit of vanishing impurities, 
i.e. the limit wpT — > oo of very long scattering times. In this case the resistivity 
and conductivity tensors read 

— \ / -Z** \ 

en Q< ) and a = ( ^ * ) , (1.9) 



en c i 



respectively. Notice that if we had put under the carpet the matrix character of 
the conductivity and resistivity and if we had only considered the longitudinal 
components, we would have come to the counter-intuitive conclusion that the 
(longitudinal) resistivity would vanish at the same time as the (longitudinal) 
conductivity. The transport properties in the clean limit loqt — * oo are therefore 
entirely governed, in the presence of a magnetic field, by the off-diagonal, i.e. 
transverse, components of the conductivity/resistivity. We will come back to 
this particular feature of quantum Hall systems when discussing the integer 
quantum Hall effect below. 



Resistivity and resistance 

The above treatment of electronic transport in the framework of the Drudc 
model allowed us to calculate the conductivity or resistivity of classical diffusive 
2D electrons in a magnetic field. However, an experimentalist does not measure 
a conductivity or resistivity, i.e. quantities that are easier to calculate for a 
theoretician, but a conductance or a resistance. Usually, these quantities are 
related to one another but depend on the geometry of the conductor - the 
resistance R is thus related to the resistivity p by R = (L/A)p, where L is the 
length of the conductor and A its cross section. From the scaling point of view 
of a G?-dimensional conductor, the cross section scales as L d_1 , such that the 
scaling relation between the resistance and the resistivity is 

R~pL 2 - d , (1.10) 

and one immediately notices that a 2D conductor is a special case. From the 
dimensional point of view, resistance and resistivity are the same in 2D, and the 
resistance is scale-invariant. Naturally, this scaling argument neglects the fact 
that the length L and the width W (the 2D cross section) do not necessarily 
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Figure 1.2: (a) Sketch of the Shubnikov-de Haas effect. Above a critical field 
B c , the longitudinal resistance (grey) starts to oscillate as a function of the 
magnetic field. The Hall resistance remains linear in B. (b) Density of states 
(DOS). In a clean system, the DOS consists of equidistant delta peaks (grey) at 
the energies e n = hujc(n + 1/2), whereas in a sample with a stronger impurity 
concentration, the peaks are broadened (dashed lines). The continuous black 
line represents the sum of overlapping peaks, and Ep denotes the Fermi energy. 



coincide: indeed, the resistance of a 2D conductor depends in general on the 
so-called aspect ratio L/W via some factor f(L/W) [5]. However, in the case 
of the transverse Hall resistance it is the length of the conductor itself that 
plays the role of the cross section, such that the Hall resistivity and the Hall 
resistance truely coincide, i.e. / = 1. We will sec in Chap. [3]that this conclusion 
also holds in the case of the quantum Hall effect and not only in the classical 
regime. Moreover, the quantum Hall effect is highly insensitive to the particular 
geometric properties of the sample used in the transport measurement, such that 
the quantisation of the Hall resistance is surprisingly precise (on the order of 
10~ 9 ) and the quantum Hall effect is used nowadays in the definition of the 
resistance standard. 



1.1.3 Shubnikov-de Haas effect 

A first indication for the relevance of quantum phenomena in transport mea- 
surements of 2D electrons in a strong magnetic field was found in 1930 with 
the discovery of the Shubnikov-de Haas effect [5J. Whereas the classical re- 
sult (|1.5[) for the resistivity tensor stipulates that the longitudinal resistivity 
Pl = 1/cto (and thus the longitudinal resistance) is independent of the magnetic 
field, Shubnikov and de Haas found that above some characteristic magnetic 
field the longitudinal resistance oscillates as a function of the magnetic field. 
This is schematically depicted in Fig. Il^f a). In contrast to this oscillation in 
the longitudinal resistance, the Hall resistance remains linear in the B field, in 
agreement with the classical result from the Drude model (|1.7[) . 

The Shubnikov-de Haas effect is a consequence of the energy quantisation of 
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the 2D electron in a strong magnetic field, as it has been shown by Landau at 
roughly the same moment. This so-called Landau quantisation will be presented 
in great detail in Sec. [2j In a nutshell, Landau quantisation consists of the 
quantisation of the cyclotron radius, i.e. the radius of the circular trajectory of 
an electron in a magnetic field. As a consequence this leads to the quantisation 
of its kinetic energy into so-called Landau levels (LLs), e n = Uu>c (n+1/2), where 
n is an integer. In order for this quantisation to be relevant, the magnetic field 
must be so strong that the electron performs at least one complete circular 
period without any collision, i.e. ujc t > 1- This condition defines the critical 
magnetic field B c ~ m^/er = fi^ 1 above which the longitudinal resistance starts 
to oscillate, in terms of the mobility (|1.6jl . Notice that today's samples of highest 
mobility are characterised by fi ~ 10 7 cm 2 /Vs = 10 3 m 2 /Vs such that one may 
obtain Shubnikov-de Haas oscillations at magnetic fields as low as B c ~ 1 mT. 

The effect may be understood within a slightly more accurate theoretical de- 
scription of electronic transport (e.g. with the help of the Boltzmann transport 
equation) than the Drude model. The resulting Einstein relation relates then 
the conductivity to a diffusion equation, and the longitudinal conductivity 



turns out to be proportional to the density of states (DOS) p{Ep) at the Fermi 
energy Ep rather than the electronic density^ Due to Landau quantisation, the 
DOS of a clean system consists of a sequence of delta peaks at the energies 
e„ = huj c (n +1/2), 



where g n is takes into account the degeneracy of the energy levels. These peaks 
are eventually impurity-broadened in real samples and may even overlap [see 
Fig. 11.21T b)] , such that the DOS oscillates in energy with maxima at the positions 
of the energy levels e„. Consider a fixed number of electrons in the sample that 
fixes the zero-field Fermi energy the .B-field dependence of which we omit in the 
argument^ When sweeping the magnetic field, one varies the energy distance 
between the LLs, and the DOS thus becomes maximal when Ep coincides with 
the energy of a LL and minimal if Ep lies between two adjacent LLs. The 
resulting oscillation in the DOS as a function of the magnetic field translates 
via the relation p. lip into an oscillation of the longitudinal conductivity (or 
resistivity), which is the essence of the Shubnikov-de Haas effect. 

1 Notice, however, that the Fermi energy and thus the DOS is a function of the electronic 
density. Furthermore we mention that in a fully consistent treatment also the diffusion con- 
stant D depends on the density of states and eventually the magnetic field. This affects the 
precise form of the oscillation but not its periodicity. 

2 Naturally, this is a crude assumption because if the density of states p(e, B) depends on 
the magnetic field, so does the Fermi energy via the relation 



However, the basic features of the Shubnikov-de Haas oscillation may be understood when 
keeping the Fermi energy constant. 



(7 L = e 2 Dp(Ep) 



(1.11) 




n 
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magnetic field B[T] 



Figure 1.3: Typical signature of the quantum Hall effect (measured by J. Smct, 
MPI-Stuttgart). Each plateau in the Hall resistance is accompanied by a van- 
ishing longitudinal resistance. The classical Hall resistance is indicated by the 
dashed-dotted line. The numbers label the plateaus: integral n denote the IQHE 
and n = p/q, with integral p and q, indicate the FQHE. 

1.1.4 Integer quantum Hall effect 

An even more striking manifestation of quantum mechanics in the transport 
properties of 2D electrons in a strong magnetic field was revealed 50 years later 
with the discovery of the integer quantum Hall effect (IQHE) by v. Klitzing, 
Dorda, and Pepper in 1980 [TU]. The Nobel Prize was attributed in 1985 to v. 
Klitzing for this extremely important discovery. 

Indeed, the discovery of the IQHE was intimitely related to technological 
advances in material science, namely in the fabrication of high-quality field- 
effect transistors for the realisation of 2D electron gases. These technological 
aspects will be briefly reviewed in separate a section (Sec. 11.21) . 

The IQHE occurs at low temperatures, when the energy scale set by the 
temperature fcgT is significantly smaller than the LL spacing tkoc- It consists 
of a quantisation of the Hall resistance, which is no longer linear in B, as one 
would expect from the classical treatment presented above, but reveals plateaus 
at particular values of the magnetic field (see Fig. II. 3p . In the plateaus, the 
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Hall resistance is given in terms of universal constants - it is indeed a fraction 
of the inverse quantum of conductance e 2 /h, and one observes 

H„ _(*)!, (1.12) 

in terms of an integer n. The plateau in the Hall resistance is accompanied 
by a vanishing longitudinal resistance. This is at first sight reminiscent of the 
Shubnikov-de Haas effect, where the longitudinal resistance also reveals minima 
although it never vanished. The vanishing of the longitudinal resistance at the 
Shubnikov-de Haas minima may indeed be used to determine the crossover from 
the Shubnikov-de Haas regime to the IQHE. 

It is noteworth to mention that the quantisation of the Hall resistance (jl.121) 
is a universal phenomenon, i.e. independent of the particular properties of the 
sample, such as its geometry, the host materials used to fabricate the 2D elec- 
tron gas and, even more importantly, its impurity concentration or distribution. 
This universality is the reason for the enormous precision of the Hall-resistance 
quantisation (typically ~ 10~ 9 ), which is nowadays - since 1990 - used as the 
resistance standard!! 

R K -90 = h/e 2 = 25 812.807O, (1.13) 

which is also called the Klitzing constant [TTJ [12]. Furthermore, as already 
mentioned in Sec. 11.1.21 the vanishing of the longitudinal resistance indicates 
that the scattering time tends to infinity [see Eq. (|1.9[) ] in the IQHE. This is 
another indication of the above-mentioned universality of the effect, i.e. that 
IQHE does not depend on a particular impurity (or scatterer) arrangement. 

A detailed presentation of the IQHE, namely the role of impurities, may be 
found in Chap. [3J 

1.1.5 Fractional quantum Hall effect 

Three years after the discovery of the IQHE, an even more unexpected effect 
was observed in a 2D electron system of higher quality, i.e. higher mobility: 
the fractional quantum Hall effect (FQHE). The effect ows its name to the fact 
that contrary to the IQHE, where the number n in Eq. (|1.12j) is an integer, a 
Hall-resistance quantisation was discovered by Tsui, Stormer and Gossard with 
n = 1/3 [13] • From the phenomenological point of view, the effect is extremely 
reminiscent of the IQHE: whereas the Hall resistance is quantised and reveals 
a plateau, the longitudinal resistance vanishes (see Fig. 11.31 where different 
instances of both the IQHE and the FQHE are shown). However, the origins of 
the two effects are completely different: whereas the IQHE may be understood 
from Landau quantisation, i.e. the kinetic-energy quantisation of independent 
electrons in a magnetic field, the FQHE is due to strong electronic correlations, 
when a LL is only partially filled and the Coulomb interaction between the 

3 The subscript K honours v. Klitzing and 90 stands for the date since which the unit of 
resistance is defined by the IQHE. 
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electrons becomes relevant. Indeed, in 1983 Laughlin showed that the origin of 
the observed FQHE with n = 1/3, as well as any n = 1/q with q being an odd 
integer, is due to the formation of a correlated incompressible electron liquid 
with extremely exotic properties [14] . which will be reviewed in Chap. [H As 
for the IQHE, the discovery and the theory of the FQHE was awarded a Nobel 
Prize (1998 for Tsui, Stormer and Laughlin). 

After the discovery of the FQHE with n = 1/30 a plethora of other types of 
FQHE has been dicovered and theoretically described. One should first mention 
the 2/5 and 3/7 states (i.e. with n = 2/5 and n = 3/7), which are part of the 
series p/(2sp± 1), with the integers s and p. This series has found a compelling 
interpretation within the so-called composite-fermion (CF) theory according to 
which the FQHE may be viewed as an IQHE of a novel quasi-particle that 
consists of an electron that "captures" an even number of flux quanta [15| 116] . 
The basis of this theory is presented in Sec. 14.41 Another intriguing FQHE was 
discovered in 1987 by Willet et al., with n = 5/2 and 7/2 [17] - it is in so far 
intriguing as up to this moment only states n = p/q with odd denominators 
had been observed in monolayer systems. From a theoretical point of view, it 
was shown in 1991 by Moore and Read [18] and by Greiter, Wilczek and Wen 
[TP] that this FQHE may be described in terms of a very particular, so-called 
PfafEan, wave function, which involves particle pairing and the excitations of 
which are anyons with non-Abelian statistics. These particles are intensively 
studied in today's research because they may play a relevant role in quantum 
computation. The physics of anyons will be introduced briefly in Sec. 14.31 
Finally, we would mention in this brief (and naturally incomplete) historical 
overview a FQHE with n — 4/11 discovered in 2003 by Pan et al. [2D]: it does 
not fit into the above-mentioned CF scries, but it would correspond to a FQHE 
of CFs rather than an IQHE of CFs. 

1.1.6 Relativistic quantum Hall effect in graphene 

Recently, quantum Hall physics experienced another unexpected boost with the 
discovery of a "relativistic" quantum Hall effect in graphene, a one-atom-thick 
layer of graphite [5J [7] . Electrons in graphene behave as if they were relativistic 
massless particles. Formally, their quantum-mechanical behaviour is no longer 
described in terms of a (non-relativistic) Schrodingcr equation, but rather by 
a relativistic 2D Dirac equation |21j . As a consequence, Landau quantisation 
of the electrons' kinetic energy turns out to be different in graphene than in 
conventional (non-relativistic) 2D electron systems, as we will discuss in Sec. [21 
This yields a "relativistic" quantum Hall effect with an unusual series for the 
Hall plateaus. Indeed rather than having plateaus with a quantised resistance 
according to Rh = h/e 2 n, with integer values of n, one finds plateaus with 
n = ±2(2n' + 1), in terms of an integer n', i.e. with n = ±2, ±6, ±10, .... The 
different signs in the series (±) indicate that there are two different carriers, 



4 The quantity n determines the filling of the LLs, usually described by the Greek letter u, 
as we will discuss in Sec. [5] 
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electrons in the conduction band and holes in the valence band, involved in the 
quantum Hall effect in graphene. As we will briefly discuss in Sec. II .21 one 
may easily change the character of the carriers in graphene with the help of the 
electric field effect. 

Interaction effects may be relevant in the formation of other integer Hall 
plateaus, such as n = and n = ±1 [22], which do not occur naturally in 
the series n = ±2(2n' + 1) characteristic of the rclativistic quantum Hall effect. 
Furthermore, a FQHE with n = 1/3 has very recently been observed, although in 
a simpler geometric (two-terminal) configuration than the standard one depicted 

in Fig. rrrr a) mCH]. 

1.2 Two-Dimensional Electron Systems 

As already mentioned above, the history of the quantum Hall effect is intimitcly 
related to technological advances in the fabrication of 2D electron systems with 
high electronic mobilities. The increasing mobility allows one to probe the fine 
structure of the Hall curve and thus to observe those quantum Hall states which 
are more fragile, such as some exotic FQHE states (e.g. the 5/2, 7/2 or the 4/11 
states). This may be compared to the quest for high resolutions in optics: the 
higher the optical resolution, the better the chance of observing tinier objects. 
In this sense, electronic mobility means resolution and the tiny object is the 
quantum Hall state. As an order of magnitude, today's best 2D electron gases 
(in GaAs/AlGaAs hcterostructures) are characterised by mobilities a, ~ 10 7 
cm 2 /Vs. 

1.2.1 Field-effect transistors 

The samples used in the discovery and in the first studies of the IQHE were so- 
calles metal-oxide-semiconductor field-effect transistors (MOSFET). A metallic 
layer is sepcrated from a semiconductor (typically doped silicon) by an insulating 
oxide (e.g. Si02) layer (see inset Jin Fig. II. 4j) . The chemical potential in the 
metallic layer may be varied with the help of a gate voltage Vq. At Vq = 0, the 
Fermi energy in the semiconductor lies in the band gap below the acceptor levels 
of the dopants [Fig. 11.4( a)]. When lowering the chemical potential in the metal 
with the help of a positive gate voltage Vg > 0, one introduces holes in the metal 
that attract, via the electric field effect, electrons from the semiconductor to the 
semiconductor-insulator interface. These electrons populate the acceptor levels, 
and as a consequence, the semiconductor bands are bent downwards when they 
approach the interface, such that the filled acceptor levels lie now below the 
Fermi energy [Fig. 11.4( b)] . 

Above a certain threshold of the gate voltage, the bending of the semicon- 
ductor bands becomes so strong that not only the acceptor levels are below the 
Fermi energy, but also the conduction band in the vicinity of the interface which 
consequently gets filled with electrons [Fig. 11.4( c)]. One thus obtains a con- 
finement potential of triangular shape for the electrons in the conduction band, 
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Figure 1.4: MOSFET. The inset I shows a sketch of a MOSFET. (a) Level 
structure at Vg = 0. In the metallic part, the band is filled up to the Fermi 
energy Ep whereas the oxide is insulating. In the semiconductor, the Fermi 
energy lies in the band gap (energy gap between the valence and the conduction 
bands). Close to the valence band, albeit above Ep, are the acceptor levels. 
(b) The chemical potential in the metallic part may be controled by the gate 
voltage Vg via the electric field effect. As a consequence of the introduction 
of holes the semiconductor bands are bent downwards, and above a threshold 
voltage (c), the conduction band is filled in the vicinity of the interface with 
the insulator. One thus obtains a 2D electron gas. Its confinement potential 
of which is of triangular shape, the levels (electronic subbands) of which are 
represented in the inset II. 
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Figure 1.5: Semiconductor heterostructure (GaAs/AlGaAs). (a) Dopants are 
introduced in the AlGaAs layer at a certain distance from the interface. The 
Fermi energy lies below the in the band gap and is pinned by the dopant levels. 
The GaAs conduction band has an energy that is lower than that of the dopant 
levels, such that it is energetically favourable for the electrons in the dopant 
layer to populate the GaAs conduction band in the vicinity of the interface, (b) 
This polarisation bends the bands in the vicinity of the interface between the 
two semiconductors, and thus a 2D electron gas is formed there on the GaAs 
side. 



the dynamics of which is quantised into discrete electronic subbands in the per- 
pendicular z-dircction (sec inset II in Fig. I1.4|l . Naturally, the electronic wave 
functions are then extended in the z-direction, but in typical MOSFETs only 
the lowest electronic subband Eq is filled, such that the electrons are purely 
2D from a dynamical point of view, i.e. there is no electronic motion in the 
z-direction. 

The typical 2D electronic densities in these systems are on the order of 
n e i ~ 10 11 cm~ 2 , i.e. much lower than in usual metals. This turns out to be 
important in the study of the IQHE and FQHE, because the effects occur, as we 
will show below, when the 2D electronic density is on the order of the density of 
magnetic flux rig = B/(h/e) threading the system, in units of the flux quantum 
h/e. This needs to be compared to metals where the surface density is on the 
order of 10 14 cm" 2 , which would require inaccessibly high magnetic fields (on 
the order of 1000 T) in order to probe the regime n e i ~ ub- 

1.2.2 Semiconductor heterostructures 

The mobility in MOSFETs, which is typically on the order of \i ~ 10 6 cm 2 /Vs, 
is limited by the quality of the oxide-semiconductor interface (surface rough- 
ness). This technical difficulty is circumvented in semiconductor heterostruc- 
tures - most popular are GaAs/ AlGaAs heterostructures - which are grown 
by molecular-beam epitaxy (MBE), where high-quality interfaces with almost 
atomic precision may be achieved, with mobilities on the order of [i ~ 10 7 
cm 2 /Vs. These mobilities were necessary to observe the FQHE, which was 
indeed first observed in a GaAs/ AlGaAs sample [13] . 
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Figure 1.6: Schematic view of graphene on a Si02 substrate with a doped Si 
(metallic) backgate. The system graphene-Si02-backgate may be viewed as a 
capacitor the charge density of which is controled by a gate voltage Vg- 

In the (generic) case of GaAs/AlGaAs, the two semiconductors do not pos- 
sess the same band gap - indeed that of GaAs is smaller than that of AlGaAs, 
which is chemically doped by donor ions at a certain distance from the interface 
between GaAs and AlGaAs [Fig 11.5( a)]. The Fermi energy is pinned by these 
donor levels in AlGaAs. which may have a higher energy than the originally un- 
occupied conduction band in the GaAs part, such that it becomes energetically 
favourable for the electrons in the donor levels to occupy the GaAs conduction 
band in the vicinity of the interface. As a consequence, the energy bands of Al- 
GaAs are bent upwards, whereas those of GaAs are bent downwards. Similarly 
to the above-mentioned MOSFET, one thus obtains a 2D electron gas at the 
interface on the GaAs side, with a triangular confinement potential. 

1.2.3 Graphene 

Graphene, a one-atom thick layer of graphite, presents a novel 2D electron 
system, which, from the electronic point of view, is either a zero-overlap semi- 
metal or a zero-gap semiconductor, where the conduction and the valence bands 
are no longer separated by an energy gap. Indeed, in the absence of doping, 
the Fermi energy lies exactly at the points where the valence band touches the 
conduction band and where the density of states vanishes linearly. 

In order to vary the Fermi energy in graphene, one usually places a graphene 
flake on a 300 nm thick insulating Si02 layer which is itself placed on top of 
a positively doped metallic silicon substrate (see Fig. II. 6p . This sandwich 
structure, with the metallic silicon layer that serves as a backgate, may thus be 
viewed as a capacitor (Fig. II. 6|) the capacitance of which is 




e eA 



(1.14) 



d 



where Q = eri2i)A is the capacitor charge, in terms of the total surface A, Vg 
is the gate voltage, and d = 300 nm is the thickness of the Si02 layer with the 
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dielectric constant e = 3.7. The field-effect induced 2D carrier density is thus 
given by 

n 2D = aV G with a = % ~ 7.2 x 10 10 ^— . (1.15) 

ea V 

The gate voltage may vary roughly between —100 and 100 V, such that one 
may induce maximal carrier densities on the order of 10 12 cm~ 2 , on top of 
the intrinsic carrier density which turns out to be zero in graphene, as will be 
discussed in the next chapter. At gate voltages above ±100 V, the capacitor 
breaks down (electrical breakdown). 

In contrast to 2D electron gases in semiconductor heterostructures, the mo- 
bilities achieved in graphene are rather low: they are typically on the order of 
/i ~ 10 4 — 10 5 cm 2 /Vs. Notice, however, that these graphene samples are fabri- 
cated in the so-called exfoliation technique, where one "peals" thin graphite crys- 
tals, under ambiant condictions, whereas the highest-mobility GaAs/AlGaAs 
laboratory samples are fabricated with a very high technological effort. The 
mobilities of graphene samples are comparable to those of commercial silicon- 
based electronic elements. 



Chapter 2 

Landau Quantisation 



The basic ingredient for the understanding of both the IQHE and the FQHE is 
Landau quantisation, i.e. the kinetic-energy quantisation of a (free) charged 2D 
particle in a perpendicular magnetic field. In this chapter, we give a detailed 
introduction to the different aspects of Landau quantisation. We have chosen a 
very general presentation of this quantisation in order to account for both a non- 
relativistic and a relativistic 2D particle some properties of which, such as the 
level degeneracy, are identical. In Sec. 12. 1[ we introduce the basic Hamiltonians 
for 2D particles in the absence of a magnetic field and discuss both Schrodinger- 
and Dirac-type particles, and discuss the case of a non-zero -B-field in Sec. 12.21 
Sec. 12.31 is devoted to the discussion of the LL structure of non-rclativistic and 
relativistic particles. 



2.1 Basic One-Particle Hamiltonians for B = 

In this section, we introduce the basic Hamiltonians which we treat in a quantum- 
mechanical manner in the following parts. Quite generally, we consider a Hamil- 
tonian for a 2D particle^] that is translation invariant, i.e. the momentum 
P = (Px,Py) is a constant of motion, in the absence of a magnetic field. In 
quantum mechanics, this means that the momentum operator commutes with 
the Hamiltonian, [p, H] = 0, and that the eigenvalue of the momentum operator 
is a good quantum number. 



2.1.1 Hamiltonian of a free particle 

In the case of a free particle, this is a very natural assumption, and one has for 
the non-relativistic case, 

P 2 

H=%-, (2-1) 
2m 



J A11 vector quantities (also in the quantum-mechanical case of operators) v = (v x , v y ) are 
hence 2D, unless stated explicitly. 
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in terms of the particle mass toU However, we are interested, here, in the motion 
of electrons in some material (in a metal or at the interface of to semiconductors). 
It seems, at first sight, to be a very crude assumption to describe the motion 
of an electron in a crystalline environment in the same manner as a particle in 
free space. Indeed, a particle in a lattice in not described by the Hamiltonian 
(|2.ip but rather by the Hamiltonian 

i 

where the last term represents the electrostatic potential caused by the ions 
situated at the lattice sites r,. Evidently, the Hamiltonian now depends on the 
position r of the particle with respect to that of the ions, and the momentum 
p is therefore no longer a constant of motion or a good quantum number. 

This problem is solved with the help of Bloch's theorem: although an ar- 
bitrary spatial translation is not an allowed symmetry operation as it is the 
case for a free particle (|2.ip . the system is invariant under a translation by an 
arbitrary lattice vector if the lattice is of infinite extension - an assumption 
we make hereH In the same manner as for the free particle, where one defines 
the momentum as the generator of a spatial translation, one may then define a 
generator of a lattice translation. This generator is called the lattice momen- 
tum or also the quasi-momcntum. As a consequence of the discreteness of the 
lattice translations, not all values of this lattice momentum are physical, but 
only those within the first Brillouin zone (BZ) - any vibrational mode, be it a 
lattice vibration or an electronic wave, with a wave vector outside the first BZ 
can be described by a mode with a wave vector within the first BZ. Since these 
lecture notes cannot include a full course on basic solid-state physics, we refer 
the reader to standard textbooks on solid-state physics US] ■ 

The bottom line is that also in a (perfect) crystal, the electrons may be 
described in terms of a Hamiltonian H(p Xl p y ) if one keeps in mind that the 
momentum p in this expression is a lattice momentum restricted to the first 
BZ. Notice, however, that although the resulting Hamiltonian may often be 
written in the form (|2.1[) . the mass is generally not the free electron mass but 
a band mass that takes into account the particular features of the energy 
bandtQ - indeed, the mass may even depend on the direction of propagation, 



2 The statement that p is a constant of motion remains valid also in the case of a relativistic 
particle. However, the Hamiltonian description depends on the frame of reference because the 
energy is not Lorentz-invariant, i.e. invariant under a transformation into another frame of 
reference that moves at constant velocity with respect to the first one. For this reason a 
Lagrangian rather than a Hamiltonian formalism is often prefered in relativistic quantum 
mechanics. 

3 Although this may seem to be a typical "theoretician's assumption", it is a very good 
approximation when the lattice size is much larger than all other relevant length scales, such 
as the lattice spacing or the Fermi wave length. 

4 In GaAs, e.g., the band mass is m b = 0.068mo, in terms of the free electron mass mo. 
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Figure 2.1: (a) Honeycomb lattice. The vectors 6±, 62, and 63 connect nn 
carbon atoms, separated by a distance a = 0.142 nm. The vectors ai and 
a.2 are basis vectors of the triangular Bravais lattice, (b) Reciprocal lattice of 
the triangular lattice. Its primitive lattice vectors are and a^. The shaded 
region represents the first Brillouin zone (BZ), with its centre T and the two 
inequivalent corners K (black squares) and K' (white squares). The thick part 
of the border of the first BZ represents those points which are counted in the 
definition such that no points are doubly counted. The first BZ, defined in a 
strict manner, is, thus, the shaded region plus the thick part of the border. For 
completeness, we have also shown the three inequivalent cristallographic points 
M, M', and M" (white triangles). 



such that one should write the Hamiltonian more generally as 

H= Pi , Pj 
2m x 2m y 



2.1.2 Dirac Hamiltonian in graphene 

The above considerations for electrons in a 2D lattice are only valid in the case 
of a Bravais lattice, i.e. a lattice in which all lattice sites are equivalent from a 
crystallographic point of view. However, some lattices, such as the honeycomb 
lattice that describes the arrangement of carbon atoms in graphene due to the 
sp 2 hybridisation of the valence electrons, are not Bravais lattices. In this case, 
one may describe the lattice as a Bravais lattice plus a particular pattern of 
N s sites, called the basis. This is illustrated in Fig. I2.1f a) for the case of 
the honeycomb lattice. When one compares a site A (full circle) with a site B 
(empty circle), one notices that the environment of these two sites is different: 
whereas a site A has nearest neighbours in the directions north-east, north-west 
and south, a site B has nearest neighbours in the directions north, south-west 
and south-east. This precisely means that the two sites are not equivalent 
from a crystallographic point of view - although they may be equivalent from 
a chemical point of view, i.e. occupied by the same atom or ion type (carbon 
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Figure 2.2: Energy bands of graphene. The valence band touches the conduc- 
tion band in the two inequivalcnt BZ corners K and K' . For undoped graphene, 
the Fermi energy lies precisely in the contact points, and the band dispersion 
in the vicinity of these points is of conical shape. 

in the case of graphene). However, all sites A form a triangular Bravais lattice 
as well as all sites B. Both subsets of lattice sites form the two sublattices, and 
the honeycomb lattice may thus be viewed as a triangular Bravais lattice with 
a two-atom basis, e.g. the pattern of two A and B sites connected by the vector 



In order to calculate the electronic bands in a lattice with N s Bravais sublat- 
tices, i.e. a basis with N s sites, one needs to describe the general electronic wave 
function as a superposition of N s different wave functions, which satisfy each 
Bloch's theorem for all sublattices Formally, this may be described in 

terms of a N s x N s matrix, the eigenvalues of which yield N s different energy 
bands. In a lattice with N s different sublattices, one therefore obtains one en- 
ergy band per sublattice, and for graphene, one obtains two different bands for 
the conducting electrons, the valence band and the conduction band. 

The Hamiltonian for low-energy electrons in reciprocal space reads 



which is obtained within a tight-binding model, where one considers electronic 
hopping between nearest-neighbouring sites with a hopping amplitude t. Be- 
cause the nearest neighbour of a site A is a site B and vice versa [see Fig. 
I2.1f a)]. the Hamiltonian is off-diagonal, and the off-diagonal elements are re- 
lated by complex conjugation due to time-reversal symmetry [H{— k)* = HQs.)]. 
As already mentioned above, the lattice momentum k is restricted to the first 
BZ, which is of hexagonal shape and which we have depicted in Fig. 12.1T b) for 




(2.3) 
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completeness. The precise form of the functions 7k is derived in Appendix [X] 
[Eq. (|A.9[) ]. The band structure is obtained by diagonalising the Hamiltonian, 
and one finds the two bands, labelled by A = ±, e\{k) = \t\jk\, which are 
plotted in Fig. 12.21 The valence band (A = — ) touches the conduction band 
(A = +) in the two inequivalent corners K and K' of the first BZ. Because there 
are as many electrons in the 7r-orbitals, that determine the low-energy conduc- 
tion properties of graphene, as lattice sites, the overall energy band structure 
is half-filled. This is due to the two spin orientations of the electrons, which 
allow for a quantum-mechanical double occupancy of each 7r-orbital. As a con- 
sequence, the Fermi energy lies exactly in the contact points K and K' of the 
two bands unless the graphene sheet is doped, e.g. with the help of the electric 
field effect, as described in Sec. 11.21 of the previous chapter. 

The inset in Fig. 12.21 shows the band dispersion in the vicinity of the con- 
tact points K and K', the linearity of which is sufficient to describe the low- 
energy electronic properties in graphene, i.e. when all relevant energy scales are 
much smaller than the full band widthH The conical form of the two bands 
is reminiscent of that of relativistic particles, the general dispersion of which 
is E = ±y // m 2 c 4 + p 2 c 2 , in terms of the light velocity c and the particle mass 
m. If the latter is zero, one obtains precisely E = ±c|p|, as in the case of low- 
energy electrons in graphene (inset of Fig. I2.2[) . which may thus be treated as 
massless Dirac fermions. Notice that in the continuum description of electrons 
in graphene, we have two electron types - one for the K point and another one 
for the K' point. This doubling is called valley degeneracy which is two-fold 
here. 

The analogy between electrons in graphene and massless relativistic particles 
is corroborated by a low-energy expansion of the Hamiltonian (|2.3[) around the 
contact points K and K', at the momenta K and K' = — K [sec Fig. I2.1f a)]. 
k = ±K+p/7i, where |p/7i| <C |K|. One may then expand the function 7±K+p/fi, 
to first order, and one obtains formallj|^| 

V V7K P J \Px+iPy J 

where cr = (a x , o y ) in terms of the Pauli matrices 

1 \ „ / -i \ , z ( 1 

1 o ) ' G = { i o ) and ° = { -1 

and where we have chosen to expand the Hamiltonian (|2.3p around the K point 
Here, the Fermi velocity v plays the role of the velocity of light c, which is though 
roughly 300 times larger, c ~ 300u. The details of the above derivation may 
be found in Appendix [X] The above Hamiltonian is indeed formally that of 
massless 2D particles, and it is sometimes called Weyl or Dirac Hamiltonian. 



5 Indeed, in graphene, the relevant low-energy scales are in the 10 — 100 meV regime, 
whereas non-linear corrections of the band dispersion become relevant in the eV regime. 
6 Notice that 7±k = by symmetry. 

7 One obtains a similar result at the K' point, see Eq. l)A.15|l in Appendix lAl 
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We will discuss, in the remainder of this chapter, how the two Hamiltonians 



^2 



H s = ^— and H D = vp ■ <x , (2.4) 

for non-relativistic and relativistic particles, respectively, need to be modified 
in order to account for a non-zero magnetic field. 



2.2 Hamiltonians for Non-Zero B Fields 

2.2.1 Minimal coupling and Peierls substitution 

In order to describe free electrons in a magnetic field, one needs to replace the 
momentum by its gauge-invariant form [27j 

p^n = p + eA(r), (2.5) 

where A(r) is the vector potential that generates the magnetic field B = V x 
A(r). This gauge-invariant momentum is proportional the electron velocity 
v, which must naturally be gauge-invariant because it is a physical quantity. 
Notice that because A(r) is not gauge invariant, neither is the momentum p. 
Remember that adding the gradiant of an arbitrary derivable function A(r), 
A(r) — > A(r)+VA(r), does not change the magnetic field because the rotational 
of a gradient is zero. Indeed, the momentum transforms as p — > p — eVA(r) 
under a gauge transformation in order to compensate the transformed vector 
potential, such that II is gauge- invariant. The substitution (|2.5[) is also called 
minimal substitution. 

In the case of electrons on a lattice, this substitution is more tricky because 
of the presence of several bands. Furthermore, the vector potential is unbound, 
even for a finite magnetic field; this becomes clear if one chooses a particular 
gauge, such as e.g. the Landau gauge Aj,(r) = B{— y, 0,0), in which case the 
value of the vector potential may become as large as B x L y , where L y is the 
macroscopic extension of the system in the y-direction. However, it may be 
shown that the substitution (|2.5[) . which is called Pcicrls substitution in the 
context of electrons on a lattice, remains correct as long as the lattice spacing 
a is much smaller than the magnetic length 




l B = \hl5> (2-6) 



which is the fundamental length scale in the presence of a magnetic field. Be- 
cause a is typically an atomic scale (~ 0.1 to 10 nm) and lg — 26 nm/ ^B[T], 
this condition is fulfilled in all atomic lattices for the magnetic fields, which may 
be achieved in today's high-field laboratories (~ 45 T in the continuous regime 
and ~ 80 T in the pulsed regime) H 



8 Higher magnetic fields may be obtained only in semi-destructive experiments, in which 
the sample survives the experiment but not the coil that is used to produce the magnetic field. 
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With the help of the (Peierls) substitution (|2.5p , one may thus immediately 
write down the Hamiltonian for charged particles in a magnetic field if one 
knows the Hamiltonian in the absence of a magnetic field, 

H(p) -> H(U) = H(p + eA) = H B (p, r). 

Notice that because of the spatial dependence of the vector potential, the result- 
ing Hamiltonian is no longer translation invariant, and the (gauge-dependent) 
momentum p is no longer a conserved quantity. We will limit the discussion to 
the B-field Hamiltonians corresponding to the Hamiltonians (|2.4|) 



ttB _ [P + eA(r)] 2 

Hs 2^ (27) 

for non-relativistic and 

HE = v[p + eA(r)] • er (2.8) 
for relativistic 2D charged particles, respectively. 



2.2.2 Quantum mechanical treatment 

In order to analyse the one-particle Hamiltonians (|2.7p and (|2.8p in a quantum- 
mechanical treatment, we use the standard method, the canonical quantisa- 
tion |28j . where one interprets the physical quantities as operators that act on 
state vectors in a Hilbert space. These operators do in general not commute 
with each other, i.e. the order matters in which they act on the state vector 
that describe the physical system. Formally one introduces the commutator 
[0i,02] = 0102 — 0201 between the two operators 0i and 02, which are 
said to commute when [0i, O2] = or else not to commute. The basic physical 
quantities in the argument of the Hamiltonian are the 2D position r = (x,y) and 
its canonical momenta p = (p x ,p y ), which satisfy the commutation relations 

[x,p x ] = ih, [y,p y ] = i% and [x,y] = [p x ,p y ] = [x,p y ] = [y,p x ] = 0, 

(2.9) 

i.e. each component of the position operator does not commute with the mo- 
mentum in the corresponding direction. This non-commutativity between the 
position and its associated momentum is the origin of the Heisenberg inequality 
according to which one cannot know precisely both the position of a quantum- 
mechanical particle and, at the same moment, its momentum, AxAp x > h and 
AyAp y > h. 

As a consequence of the commutation relations (|2.9[) . the components of the 
gauge-invariant momentum no longer commute themselves, 

[n^Hy] = [p x + eA x (r),p v + eA y (r)}=e([p x ,A y }~[p y ,A x }) 

(dAy dAy dA dA A 
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where we have used the relation^ 



df 



[01,02] 



(2.10) 



d0 2 



between two arbitrary operators, the commutator of which is a c-number or an 
operator that commutes itself with both 0\ and O2 [28]. With the help of the 
commutation relations (|2.9[) . one finds that 



and, in terms of the magnetic length (|2.6p , 

[1^,1^] = -^ . (2.11) 

This equation is the basic result of this section and merits some further discus- 
sion. 

• As one would have expected for gauge-invariant quantities (the two com- 
ponents of II), their commutator is itself gauge-invariant. Indeed, it only 
depends on universal constants and the (gauge-invariant) magnetic field 
B, and not on the vector potential A. 

• The components of the gauge-invariant momentum II are mutually conju- 
gate in the same manner as x and p x or y and p y . Remember that p x gen- 
erates the translations in the x-direction (andp y those in the y-dircction). 
This is similar here: IL X generates a "boost" of the gauge-invariant mo- 
mentum in the y-direction, and similarly IL y one in the x-dircction. 

• As a consequence, one may not diagonalise at the same time Tl x and H y , 
in contrast to the zero-field case, where the arguments of the Hamiltonian, 
p x and p y , commute. 

For solving the Hamiltonians (|2.7p and (|2.8p , it is convenient to use the pair 
of conjugate operators IL^ and H y to introduce ladder operators in the same 
manner as in the quantum-mechanical treatment of the one-dimensional har- 
monic oscillator. Remember from your basic quantum-mechanics class that the 
ladder operators may be viewed as the complex position of the one-dimensional 
oscillator in the phase space, which is spanned by the position (x-axis) and the 
momentum (y-axis), 



9 More precisely we have used a gradient generalisation of this relation to operator functions 
that depend on several different operators, 






and 




which is valid if [[Co, Oj],O ] = [[O , Oj],Oj] 



[O ,f(Oi,...,Oj)] 
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where xq = ^JhjrribW and po = s/fvmtfj) are normalisation constants in terms of 
the oscillator frequency u> |28j . The fact that the position x and the momentum p 
are conjugate variables and the particular choice of the normalisation constants 
yields the commutation relation [a, at] = 1 for the ladder operators. 

In the case of the 2D electron in a magnetic field, the ladder operators play 
the role of a complex gauge- invariant momentum (or velocity) , and they read 

R (U x -illy) and a f = (n x + illy) , (2.12) 



V2h " y/2h 

where we have chosen the appropriate normalisation such as to obtain the usual 
commutation relation 

[a,a t ] = l. (2.13) 

It turns out to be helpful for future calculations to invert the expression for the 
ladder operators (|2.12D . 

n, = -A- (at + a) and IL, = -J- ( a t - a) . (2.14) 



2.3 Landau Levels 



The considerations of the preceding section are extremely useful in the calcu- 
lation of the level spectrum associated with the Hamiltonians (|2.7p and (|2.8[) 
of both the non-relativistic and the relativistic particles, respectively. The un- 
derstanding of this level spectrum is the issue of the present section. Because 
electrons do not only possess a charge but also a spin, each level is split into 
two spin branches separated by the energy difference A^e = g^sB, where g is 
the g-factor of the host material and [Ib = eh/2mo the Bohr magneton. In or- 
der to simplify the following presentation of the quantum-mechanical treatment 
and the level structure, we neglect this effect associated with the spin degree of 
freedom. Formally, this amounts to considering spinlcss fcrmions. Notice, how- 
ever, that there exist interesting physical properties related to the spin degree 
of freedom, which will be treated separately in Chap. [5] 



2.3.1 Non-relativistic Landau levels 

In terms of the gauge-invariant momentum, the Hamiltonian (|2.7p for non- 
relativistic electrons reads 

The analogy with the one-dimensional harmonic oscillator is apparent if one 
notices that both conjugate operators Tl x and 11^ occur in this expression in a 
quadratic form. If one replaces these operators with the ladder operators (|2.14[) . 
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one obtains, with the help of the commutation relation (|2.13p , 

h 2 

H§ = 4m ^ 2 [ flt2 + atfl + aqt + q2 ~ ( flt2 ~ flta ~ aqt + fl2 )] 
— — - — (ata J- - — ( 



aa' = — -77- a' a H — 



(2.15) 



where we have used the relation uj c = h/rribl B between the cyclotron frequency 
(|1.2|) and the magnetic length (|2.6|) in the last step. 

As in the case of the one-dimensional harmonic oscillator, the eigenvalues 
and eigenstates of the Hamiltonian (|2 . 1 5[) are therefore those of the number 
operator a^a, with a)a\n) = n\n). The ladder operators act on these states in 
the usual manner 1281 



a T|n) = y/n + l\n + 1) and a\n) = y/n\n - 1), (2.16) 

where the last equation is valid only for n > - the action of a on the ground 
state |0) gives zero, 

a|0)=0. (2.17) 

This last equation turns out to be helpful in the calculation of the eigenstates 
associated with the level of lowest energy, as well as the construction of states 
in higher levels n (see Sec. I2.4.1[) 

■|0). (2.18) 



The energy levels of the 2D charged non-relativistic particle are therefore 
discrete and labelled by the integer n, 

e n = hu> c (n+^\ . (2.19) 

These levels, which are also called Landau levels (LL), are depicted in Fig. 12.3( a) 
as a function of the magnetic field. Because of the linear field-dependence of 
the cyclotron frequency, the LLs disperse linearly themselves with the magnetic 
field. 

2.3.2 Relativistic Landau levels 

The relativistic case (|2.8|) for electrons in graphene may be treated exactly in 
the same manner as the non-relativistic one. In terms of the ladder operators 
(|2.12| . the Hamiltonian reads 



IT- + illy J ^ l B \ ^ 
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magnetic field B 

Figure 2.3: Landau levels as a function of the magnetic field, (a) Non-relativistic 
case with e n = hu>c(n + 1/2) oc B(n + 1/2). (b) Relativistic case with e\, n — 
\(hv/lB)V2n cx AV Bn. 

One notices the occurence of a characteristic frequency u>' = \/2v/Ib, which 
plays the role of the cyclotron frequency in the relativistic case. Notice, however, 
that this frequency may not be written in the form eB/rrib because the band 
mass is strictly zero in graphene, such that the frequency would diverge F°l 

In order to obtain the eigenvalues and the eigenstates of the Hamiltonian 
(|2.20p , one needs to solve the eigenvalue equation Hpip n = e n ip n . Because the 
Hamiltonian is a 2 x 2 matrix, the eigenstates are 2-spinors, 




and we thus need to solve the system of equations 

TiLo'av n — e n u n and huj' a) u n = e n v n , (2.21) 
which yields the equation 




(2.22) 



for the second spinor component. One notices that this component is an eigen- 
state of the number operator n = a^a, which we have already encountered in 
the preceding subsection. We may therefore identify, up to a numerical factor, 
the second spinor component v n with the eigenstate |n) of the non-relativistic 
Hamiltonian (|2.15p . v n ~ \n). Furthermore, one observes that the square of the 

10 Sometimes, a cyclotron mass rriQ is formally introduced via the equality ui' = eB/rriQ. 
However, this mass is a somewhat artificial quantity, which turns out to depend on the carrier 
density. We will therefore not use this quantity in the present lecture notes. 
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energy is proportional to this quantum number, = (huj') 2 n. This equation 
has two solutions, a positive and a negative one, and one needs to introduce 
another quantum number A = ±, which labels the states of positive and neg- 
ative energy, respectively. This quantum number plays the same role as the 
band index (A = + for the conduction and A = — for the valence band) in the 
zero-B-field case discussed in Sec. 12.11 One thus obtains the level spectrum 



the energy levels of which are depicted in Fig. I2.3f b) . These rclativistic Landau 
levels disperse as AV Bn as a function of the magnetic field. 

Once wc know the second spinor component, the first spinor component is 
obtained from Eq. (|2.21[) . which reads u n oc av n ~ a\ri) ~ \n — 1). One then 
needs to distinguish the zero-energy LL (n = 0) from all other levels. Indeed, 
for n = 0, the first component is zero as one may see from Eq. (|2.17[) . In this 
case one obtains the spinor 



In all other cases (n ^ 0), one has positive and negative energy solutions, 
which differ among each other by a relative sign in one of the components. A 
convenient representation of the associated spinors is given by 



Experimental observation of relativistic Landau levels 

Relativistic LLs have been observed experimentally in transmission spectroscopy, 
where one shines light on the sample and measures the intensity of the trans- 
mitted light. Such experiments have been performed on so-called epitaxial 
graphcneO [31] and later on exfoliated graphene [32]. When the monochro- 
matic light is in resonance with a dipole-allowed transition from the (partially) 
filled LL (A, n) to the (partially) unoccupied LL (A',rt ± 1), the light is ab- 
sorbed due to an electronic excitation between the two levels [see Fig. I2.4f a)]. 
Notice that, in a non-relativistic 2D electron gas, the only allowed dipolar tran- 
sition is that from the last occupied LL n to the first unoccupied one n + 1. 
The transition energy is Tiu>c, independently of n, and one therefore observes 
a single absorbtion line (cyclotron resonance). In graphene, however, there are 
many more allowed transitions due to the presence of two electronic bands, the 
conduction and the valence band, and the transitions have the energies 




(2.23) 




(2.24) 




(2.25) 




y/2(n + 1) -£V2n~ , 



11 Epitaxial graphene is obtained from a thermal graphitisation process of an epitaxially 
grown SiC crystal I3U| 
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Figure 2.4: LL spectroscopy in graphene (from Sadowski et al, 2006). (a) 
For a fixed magnetic field (0.4 T), one observes resonances in the transmission 
spectrum as a function of the irradiation energy. The resonances are associated 
with allowed dipolar transitions between rclativistic LLs. (b) These resonances 
are shifted as a function of the magnetic field, (c) If one plots the resonance 
energies as a function of the square root of the magnetic field, \f~B 1 a linear 
dependence is observed as one would expect for relativistic LLs. 
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where £ = + denotes an intraband and £ = — an interband transition. One 
therefore obtains families of resonances the energy of which disperses as A n ^ oc 
\AB, as it has been observed in the experiments [see Fig. I2.4f c). where we show 
the results from Sadowski et al. [31]]. Notice that the dashed lines in Fig. 12.4( c) 
are fits with a single fitting parameter (the Fermi velocity v), which matches 
well all experimental points for different values of n. 



2.3.3 Level degeneracy 

In the preceding subsection, we have learnt that the energy of 2D (non-)relativistic 
charged particles is characterised by a quantum number n, which denotes the 
LLs (in addition to the band index A in for rclativistic particles). However, 
the quantum system is yet underdetermined, as may be seen from the following 
dimensional argument. The original Hamiltonians (|2.7[) and (|2.8[) are functions 
that depend on two pairs of conjugate operators, x and p x , and y and p y , 
respectively, whereas when they are expressed in terms of the gauge-invariant 
momentum II or else the ladder operators a and the Hamiltonians (|2.15j) and 
(|2.20[) depend only on a single pair of conjugate operators. From the original 
models, one would therefore expect the quantum states to be described by two 
quantum numbers (one for each spatial dimension). This is indeed the case in 
the zero-field models (|2.4|) . where the quantum states are characterised by the 
two quantum numbers p x and p y , i.e. the components of the 2D momentum. 
For a complete description of the quantum states, we must therefore search 
for a second pair of conjugate operators, which necessarily commutes with the 
Hamiltonian and which therefore gives rise to the level degeneracy of the LLs 
- in addition to the degeneracy due to internal degrees of freedom such as the 
spiral or, in the case of graphenc, the two-fold valley degeneracy. 

In analogy with the gauge- invariant momentum, II = p + eA(r), we consider 
the same combination with the opposite relative sign, 

II = p-eA(r), (2.26) 

which we call pseudo-momentum to give a name to this operator. One may 
then express the momentum operator p and the vector potential A(r) in terms 
of II and n, 

p=i(n + n) and A(r) = i(II-n). (2.27) 
2 2e 

Notice that, in contrast to the gauge-invariant momentum, the pseudo-momentum 
depends on the gauge and, therefore, docs not represent a physical quantity^ 
However, the commutator between the two components of the pseudo- momentum 
turn out to be gauge- invariant, 

h 2 



n„n. 



*7T . (2.28) 

I 1 



B 



12 The quantum states are naturally only degenerate if one neglects the Zeeman effect. 
13 We will nevertheless try to give a physical interpretation to this operator below, within a 



semi-classical picture. 
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This expression is calculated in the same manner as the commutator (|2.11[) 
between IT^ and U y , as well as the mixed commutators between the gauge- 
invariant momentum and the pseudo-momentum, 





= 2ieh- 




„ . c 


Uy , Uy 


= zien- 




= ieh 1 



dx 



dy 

dA x dA v 



(2.29) 



dy 



0;> 



n K ,iL 



These mixed commutators arc unwanted quantities because they would in- 
duce unphysical dynamics due to the fact that the components of the pseudo- 
momentum would not commute with the Hamiltonian, [Il^/y, H] ^ 0. However, 
this embarrassing situation may be avoided by choosing the appropriate gauge, 
which turns out to be the symmetric gauge 



A s(r) = —{-y,x,0), 



(2.30) 



with the help of which all mixed commutators (|2.29|) vanish such that the com- 
ponents of the pseudo-momentum also commute with the Hamiltonian. 

Notice that there exists a second popular choice for the vector potential, 
namely the Landau gauge, which we have already mentioned above, 



A £ (r)=B(-y,0,0), 



(2.31) 



for which the last of the mixed commutators (|2.29|) would not vanish. This 
gauge choice may even occur simpler: because the vector potential only de- 
pends on the y-componcnt of the position, the system remains then translation 
invariant in the a;-direction. Therefore, the associated momentum p x is a good 
quantum number, which may be used to label the quantum states in addition 
to the LL quantum number n. For the Landau gauge, which is useful in the de- 
scription of geometries with translation invariance in the y-direction, the wave 
functions are calculated in Sec. ()2.4.2|) . However, the symmetric gauge, the 
wave functions of which are presented in Sec. (|2.4.1[) . plays an important role 
in two different aspects; first, it allows for a semi-classical interpretation more 
easily than the Landau gauge, and second, the wave functions obtained from 
the symmetric gauge happen to be the basic ingredient in the construction of 
trial wave functions a la Laughlin for the description of the FQHE, as we will 
see in Chap. [U 

The pseudo- momentum, with its mutually conjugate components and Hi,, 
allows us to introduce, in the same manner as for the gauge-invariant momentum 
II, ladder operators, 



b = 



Ib 

V2h 



n, 



in 



and 



Ib 

V2h 



Tlx - «n 



(2.32) 
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Figure 2.5: Cyclotron motion of an electron in a magnetic field around the 
guiding centre R. The grey region indicates the quantum-mechanical uncer- 
tainty of the guiding-centre position due to the non-commutativity (|2.39p of its 
components. 



which again satisfy the usual commutation relations [b, tf] = 1 and which, in the 
symmetric gauge, commute with the ladder operators a and a 1 , [b, a^j = 0, and 
thus with the Hamiltonian, [6^ , -ffg] = 0. One may then introduce a number 
operator Wb associated with these ladder operators, the eigenstates of which 
satisfy the eigenvalue equation 

b^b\m) = m|m). 

One thus obtains a second quantum number, an integer m > 0, which is nec- 
essary to describe, as expected from the above dimensional argument, the full 
quantum states in addition to the LL quantum number n. The quantum states 
therefore become tensor products of the two Hilbert vectors 

\n,m) = \n) <g> |m) (2.33) 

for non-relativistic particles. In the relativistic case, one has 

1 / \Tl — 1 77 1) \ 
1p\n, m = 1p\n, m ® | m) = — = ' , '/ (2.34) 



for n and 



A|n, to) 



V>™=o,m = V>™=o ® \m) = ( , = ° . ) (2.35) 



for the zero-energy LL. 



2.3.4 Semi-classical interpretation of the level degeneracy 

How can we illustrate this somewhat mysterious pseudo-momentum introduced 
formally above? Remember that, because the pseudo- momentum is a gauge- 
dependent quantity, any physical interpretation needs to be handled with care. 
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However, within a semiclassical treatment, the symmetric gauge allows us to 
make a connection with a classical constant of motion that one obtains from 
solving the classical equations of motion for a massive electron in a magnetic 
field, 

f x = -uj c y 

m b r = -e(r x B) <s> <^ (2.36) 

which is nothing other than the electron's acceleration due to the Lorentz force. 
These equations may be integrated, and one then finds 

= -"o(y-Y) ] ( y = 

^ n (2 - 37) 

where R = (X, Y) is an integration constant, which physically describes a con- 
stant of motion. This quantity may easily be interpreted: it represents the centre 
of the electronic cyclotron motion (see Fig. 12. 5|) . Indeed, further integration of 
the equations (|2.37|) yields the classical cyclotron motion 

x(t) = X — r sin(uJct + (f) and y(t) = Y + r cos(uct + </>), 

where the phase 4> is another constant of motion. The cyclotron motion itself 
is described by the velocities (or else the gauge-invariant momenta) H x /m and 
ILy/m. Whereas the energy depends on these velocities that determine the 
radius r of the cyclotron motion, it is completely independent of the position 
of its centre R, which we call guiding centre from now on, as one would expect 
from the translational invariance of the equations of motion (|2.36p . 

In order to relate the guiding centre R to the pseudo-momentum II, we use 
Eq. (|2.27p for the vector potential in the symmetric gauge, 

•«*)-" {7 )->-**■ 

The postions x and y may then be expressed in terms of the momenta II and 

n, 



fix 




eB 


eB 






eB 


eB 



A comparison of these expresssions with Eq. (|2.37|) allows us to identify 

X = -5* and Y = — . (2.38) 
eB eB V ' 

This means that, in the symmetric gauge, the components of the pseudo- 
momentum arc nothing other, apart from a factor to translate a momentum 
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into a position, than the the components of the guiding centre, which are natu- 
rally constants of motion. In a quantum- mechanical treatment, these operators 
therefore necessarily commute with the Hamiltonian, as we have already seen 
above. Furthermore, the commutation relation (|2.28|) between the components 
of the pseudo-momentum, [n x ,rij,] = ih 2 /l g induces the commutation relation 

[X,Y]=il% (2.39) 

between the components of the guiding-centre operator. This means that there 
is a Heisenberg uncertainty associated with the guiding-centre position of a 
quantum-mechanical state - one cannot know its x- and y-components simulta- 
neously, and the guiding centre is, therefore, smeared out over a surface 

MAY = 2ttZ| (2.40) 

(see grey region in Fig. 12.5ft F^l This minimal surface plays the same role as the 
surface (action) h in phase space and therefore allows us to count the number 
of possible quantum states of a given (macroscopic) surface A, 

NB =AYKY = 2^ = nBXA > 
where we have introduced the flux density 

nB - 2k = w (2 - 41) 

which is nothing other than the magnetic field measured in units of the flux 
quantum h/e. Therefore, the number of quantum states in a LL equals the 
number of Rux quanta threading the sample surface A, and each LL is macro- 
scopically degenerate. We will show in a more quantitative manner than in the 
above argument based on the Heisenberg inequality that the number of states 
per LL is indeed given by Nb when discussing, in the next section, the electronic 
wave functions in the symmetric and the Landau gauges. 

Similarly to the guiding-centre operator, we may introduce the cyclotron 
variable r\ = (rj x ,r]y), which determines the cyclotron motion and which fully 
describes the dynamic properties. The cyclotron variable is perpendicular to 
the electron's velocity and may be expressed in terms of the gauge-invariant 
momentum II, 

Vx = 5| and »?» = —, (2.42) 

as one sees from Eq. (|2.37j) . The position of the electron is thus decomposed into 
its guiding centre and its cyclotron variable, r = R + rj. Also the components 
of the cyclotron variable do not commute, and one finds with the help of Eq. 

arm 

lV*,Vy}= [ j^ = -ilB = -[X,Y]. (2.43) 



14 Mathematicians speak of a non-commutative geometry in this context, and the charged 
2D particle may be viewed as a pardigm of this concept. 
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Until now, we have only discussed a single particle and its possible quantum 
states. Consider now N independent quantum-mechanical electrons at zero- 
temperature. In the absence of a magnetic field, electrons in a metal, due to 
their fermionic nature and the Pauli principle which prohibits double occupancy 
of a single quantum state, fill all quantum states up to the Fermi energy, which 
depends thus on the number of electrons itself. The situation is similar in the 
presence of a magnetic field: the electrons preferentially occupy the lowest LLs, 
i.e. those of lowest energy. But once a LL is filled, the remaining electrons are 
forced to populate higher LLs. In order to describe the LL filling it is therefore 
useful to introduce the dimensionlcss ratio between the number of electrons 
N e i = n e i x A and that of the flux quanta, 

N e i n e i hn e i . , 

;/ = — = — = — -, 2.44) 
N B n B eB ' y ' 

which is called filling factor. Indeed the integer part, [u], of the filling factor 
counts the number of completely filled LLs. Notice that one may vary the filling 
factor either by changing the particle number or by changing the magnetic field. 
At fixed particle number, lowering the magnetic field corresponds to an increase 
of the filling factor. 



2.4 Eigenstates 

2.4.1 Wave functions in the symmetric gauge 

The algebraic tools established above may be used calculate the electronic wave 
functions, which are the space representations of the quantum states |n, to), 
0n,m(^, y) = (x, y\n, to) 111 Notice first that one may obtain all quantum state 
|n, m) from a single state \n = 0, m = 0), with the help of 

(at)" (&t) TO 

\n,m) = ±-^±-=!=\n = 0,to = 0), (2.45) 



which is a generalisation of Eq. (|2.18|) . Naturally, this equation translates into 
a differential equation for the wave functions 4> n ,m{x, y). 

A state in the lowest LL (n = 0) is characterised by the condition (|2.f 71) 

a\n = 0,m)=0, (2.46) 

which needs to be translated into a differential equation. Remember from Eq. 
(|2~T2"|) that a = (l B /y/2h)(Il x - iU y ) and, by definition, II = -iKV + eA(r) 
where we have already represented the momentum as a differential operator in 
position representation, p = — zftV. One then finds 



iV2 



Ib ,r, x-iy 



15 We limit the discussion to the non-relativistic case. The spinor wave functions for rela- 
tivistic electrons are then easily obtained with the help of Eqs. H2.34I I and 1 12.351 . 
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where d x and d y are the components of the gradient V = (d x , d y ), and one sees 
from this expression that it is convenient to introduce complex coordinates to 
describe the 2D plane. We define z = x — iy, z* = x + iy, d = (d x + id y )/2 and 
d = (d x — id y )/2. The lowest LL condition (|2.46j) then becomes a differential 
equation, 

4^+^) <t>n=o(z,z*) = 0, (2.47) 
which may easily be solved by the complex function 

cj> n=0 (z,z*)=f(z)e-W il °, (2.48) 

where f(z) is an analytic function, i.e. df(z) = 0, and \z\ 2 — zz* . This means 
that there is an additional degree of freedom because f(z) may be any analytic 
function. It is not unexpected that this degree of freedom is associated with the 
second quantum number to, as we will now discuss. 

The ladder operators b and may be expressed in position representation in 
a similar manner as a, and one obtains the space representation of the different 
ladder operators, 

V2~(^- + lBd), ^ = iy/2 ( 4r~ - iBd 



V2(|- + / B c)j, =iV2~(^--l B 5) . (2.49) 



In the same manner as for a state in the lowest LL, the condition for the reference 
state with m = is b\n, m = 0) = 0, which yields the differential equation 

{z* + 4l 2 B d) ti n=o (z,z*) = 

with the solution 

<t>' m=0 {z,z*)=g{z*)e-W 2 t il °, 

in terms of an anti-analytic function g(z*) with dg(z*) = 0. The wave function 
4>n=o,m=o(z, z*) must therefore be the Gaussian with a prefactor that is both 
analytic and anti-analytic, i.e. a constant that is fixed by the normalisation. 
One finds 

^=o,m=o(2,2*) = (z,z> = 0,to = 0) = — i=e-l z l 2 / 4 ^, (2.50) 

V 27t1 b 

and a lowest-LL state with arbitrary to may then be obtained with the help of 

Eq. msD , 
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The states within the lowest LL are therefore, apart from the Gaussian, given 
by the usual polynomial basis states z m of analytic functions. In an arbitrary 
LL, the states may be obtained in a similar manner, but they happen to be 
more complicated because the differential operators (|2.49p no longer act on the 
Gaussian only but also on the polynomial functions. They may be expressed in 
terms on Laguerre polynomials. 

To conclude the discussion about the wave functions in the symmetric gauge, 
we calculate the average value of the guiding-centre operator in the state \n = 
0, to). With the help of Eqs. (|2.32p and ()2.38j) . one may express the components 
of the guiding-centre operator in terms of the ladder operators b and w , 

X = j2=(bt-b) and F=i|( & t + 6)) (2.52) 

and the ladder operators act, in analogy with Eq. (|2.16p , on the states \n, m) 
as 

£>' \n, to) = \/m + l\n, to + 1) and 6|n, to) = y/m\n, m — 1). 

The average value of the guiding-centre operator is therefore zero in the states 
|n, to), 

(R) = (n = 0, m\K\n = 0, to.) = 0, 

but we have 

(|R|) = (yx 2 + Y^ = l B (y/2tfb+l} = W2to+1. (2.53) 

This means that the guiding centre is situated, in a quantum state \n, to), some- 
where on a circle of radius I b \j2m + 1 whereas its angle (or phase) is completely 
undetermined. 

The symmetric gauge is the natural gauge to describe a sample in the form 
of a disc. Consider the disc to have a radius R ma x (and a surface A = TrR^ nax ). 
How many quantum states may be accomodated within the circle? The quan- 
tum state with maximal to quantum number, which we call M, has a radius 
IbV^M + 1, which must naturally coincide with the radius Rmax of the disc. 
One therefore obtains A = irl B (2Al +1), and the number of states within the 
disc is then, in the thermodynamic limit M 1, 

M = -4r =n B xA = N B , (2.54) 

in agreement with the result (|2.41[) obtained from the argument based on the 
Heisenberg uncertainty relation. 

2.4.2 Wave functions in the Landau gauge 

If the sample geometry is rectangular, the Landau gauge (|2.3ip . Al(t) = 
B{—y 1 0, 0), is more appropriate than the symmetric gauge to describe the phys- 
ical system. As already mentioned above, the momentum p x = fik is a good 
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quantum number due to translational invariance in the x-direction. One may 
therefore use a plane-wave ansatz 

i> n .k{x,y) = —7=Xn,k(v)i 

v L 



for the wave functions. In this case, the Hamiltonian (|2.7jl becomes 

Hs = 2^ + 2^ = 2^ + 2 m " c(y Vt)) ' (2 ' 55) 



2/o = kl%. (2.56) 



where we have defined 

The Hamiltonian (|2 . 55[) is nothing other than the Hamiltonian of a one-dimensional 
oscillator centred around the position yo, and the eigenstates are 



in terms of Hermite polynomials H n (x) [28| . The coordinate yo plays the role 
of the guiding centre component Y, the component X being smeared over the 
whole sample length L, as it is dictated by the Heisenberg uncertainty relation 
resulting from the commutation relation ()2.39p [X, Y] = il B . 

Using periodic boundary conditions k = m x 2tt/L for the wave vector in 
the IE-direction, one may count the number of states in a rectangular surface of 
length L and width W (in the y-direction) , similarly to the above arguments in 
the symmetric gauge. Consider the sample to range from y m { n = to y max = 
W, the first corresponding via the above-mentioned condition (|2.56p to the 
wave vector k = and the latter to a wave vector k max = M x lixjL. Two 
neighbouring quantum states are separated by the distance Ay = Akl B = 
Am(2ir / L)l% = 2irl B /L, and each state therefore occupies a surface a = Ay x 
L = 2nl%, which agrees with the result (|2.40j) obtained above with the help of the 
consideration based on the Heisenberg uncertainty relation. The total number 
of states is, as in the symmetric gauge and the general argument leading to Eq. 

(EUD, 

M = N B = n B x LW = n B x A, 
i.e. the number of flux quanta threading the (rectangular) surface A = LW. 
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The quantum-mechanical treatment of the 2D electron in a perpendicular mag- 
netic field is the backbone for the understanding of the basic properties of the 
quantum Hall effect. However, we need to relate the kinetic-energy quantisa- 
tion to the resistance quantisation, which is the essential feature of the IQHE. 
In the present chapter, we discuss the transport properties of electrons in the 
IQHE, namely the somewhat mysterious role that disorder plays in this type of 
transport. Remember from the introduction that the Hall resistance is quan- 
tised with an astonishingly high precision (10 -9 ), such that it is now used as the 
standard of resistance [see Eq. (| 1 . 1 3[) ] . The resistance quantisation in the IQHE 
therefore does reflect neither a particular disorder distribution nor a particular 
sample geometry Nevertheless, disorder turns out to play an essential role in 
the occurence of the IQHE, as we will see in this chapter. 

We will first consider, in Sec. 13. 1[ the motion of a 2D electron in a perpen- 
dicular magnetic field when also an external electrostatic potential is present, 
such as the one generated by disorder or the confinement potential that defines 
the sample boundaries. In Sec. 13.21 we then calculate the conductance of a 
single LL within a mesoscopic picture and discuss the difference between a two- 
terminal and a six-terminal transport measurement in Sec. 13.31 Furthermore, 
we discuss, in Sec. 13.41 the IQHE within a percolation picture and present some 
scaling properties that characterise the plateau transitions. We terminate this 
chapter with a short discussion of the pecularities of the relativistic quantum 
Hall effect in graphene the understanding of which requires essentially the same 
ingredients as the IQHE in non-relativistic quantum Hall systems. 
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Figure 3.1: Potential landscape of an electrostatic potential in a sample. The 
metallic contacts are described by the chemical potentials [i^ and for the 
left and right contacts, respectively. We consider L 3> W » £ S> Ib, where £ 
is the typical length scale for the variation of the electrostatic potential. The 
sample is confined in the y-direction between y max and y m in- The thin lines 
indicate the equipotential lines. When approaching one of the sample edges, 
they become parallel to the edge. The grey lines indicate the electronic motion 
with the guiding centre moving along the equipotential lines. The electron turns 
clockwise around a summit of the potential landscape, which is caused e.g. by 
a negatively charged impurity (— ), and counter-clockwise around a valley (+). 
At the sample edges, the equipotential lines due to the confinement potential 
connect the two contacts on the left and on the right hand side. 
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3.1 Electronic Motion in an External Electro- 
static Potential 

We consider a system the length L of which is much larger than the width W 
(see Fig. 13. ip . This may be modeled by a confinement potential V con f(y) that 
only depends on the y-direction, i.e. the system remains translation-invariant in 
the x-direction with respect to this potential!]] In addition to the confinement, 
we consider a smoothly varying electrostatic potential V{ m p{x,y) that is caused 
by the impurities in the sample. This impurity potential breaks the translation 
invariancc in the x-dircction as well as that in the y-direction, which is already 
broken by the confinement potential. The Hamiltonian of a 2D particle in a 
perpendicular magnetic field then needs to be completed by a potential term 

V(r) = V con{ (y)+V imp (x,y), (3.1) 

which creates a potential landscape that is schematically depicted in Fig. 13.11 

3.1.1 Semi-classical treatment 

In a first step, we consider a potential V(r) that varies smoothly on the scale 
set by the magnetic length, i.e. £ 3> Ib, where £ describes the characteristic 
length scale for the variation of V(r). Notice first that the external electrostatic 
potential lifts the LL degeneracy because the Hamiltonian H = Hb + V(r = 
R.+T7) no longer commutes with the guiding-centre operator R, in contrast to the 
"free" Hamiltonian Hb, [H, R] = [V, R] 7^ 0. Physically, this is not unexpected: 
the guiding centre is a constant of motion due to translation invariancc, i.e. it 
does not matter whether the electron performs its cyclotron motion around a 
point Ri or R2 in the 2D plane as long as the cyclotron radius is the same. 
However, the electrostatic potential V(r) breaks this translation invariance and 
thus lifts the degeneracy associated with the guiding centre. 

In the case where the electrostatic potential varies smoothly on a length 
scale set by the magnetic length and does not generate LL mixing, i.e. when 
|VV| <C Hu>c/Ib, we may approximate the argument r in the potential (|3.ip by 
the guiding-centre variable RH 

V(v) ~ V(R). (3.2) 

Notice that this approximation may seem unappropriate when we consider the 
confinement potential in the y-direction which may vary abruptly when ap- 

1 Naturally, the system is also confined in the x-direction, but since we consider a sample 
with L 3> W, the system appears as translation-invariant in the x-direction when one considers 
intermediate length scales. The latter may be taken into account with the help of periodic 
boundary conditions that discretise the wave vector in the ^-direction, as we have seen in the 
preceding chapter within the quantum-mechanical treatment of the 2D electron in the Landau 
gauge (see Sec. 12. 4. 21 . 

2 This approximation may be viewed as the first term of an expansion of the electrostatic 
potential in the coherent (or vortex) state basis, where the states are maximally localised 
around the guiding-centre position R. 1331 . 
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proaching the sample edges. The confinement potential with translation invari- 
ance in the a;-direction will be discussed separately in the following subsection. 

As a consequence of the non-commutativity of the potential term V(R) with 
the guiding-centre operator, the latter quantity acquires dynamics, as may be 
seen from the Heisenberg equations of motion 



dV lYV] , ]2 dV 

3V 

ihY = [Y,V(R)]=-il% 1 r V , (3.3) 



ihX = [X,H] = [X,V(R)) = ^[X,Y] =V 2 

dV 



■dX 1 

i.e. R _L VV . Here, we have used the commutation relation (|2.39[) for the 
guiding-centre components and Eq. (|2.10p . The Heisenberg equations of motion 
are particularly useful in the discussion of the semi-classical limit because the 
averaged equations satisfy the classical equations of motion, 

(R) _L W, (3.4) 

which means that, within the semi-classical picture, the guiding centres move 
along the equipotcntial lines of the smoothly varying external electrostatic po- 
tential. This feature, which is also called the Hall drift, 

ExB ■ -WxB 

VD = = ( R > = eB 2 ■ ( 3 - 5 ) 

in terms of the (local) electric field E — — VV/e, is depicted in Fig. 13. II by the 
grey lines. 

In the bulk, the potential landscape is created by the charged impurities in 
the sample, and the electrons turn clockwise on an equipotcntial line around a 
summit that is caused by a negatively charged impurity and counter-clockwise 
around a valley created by a positively charged impurity. If the equipotential 
lines are closed, as it is the case for most of the equipotential lines in a poten- 
tial landscapcjj an electron cannot move from one point to another one over a 
macroscopic distance, e.g. from one contact to the other one. An electron mov- 
ing on a closed equipotential line can therefore not contribute to the electronic 
transport, and the electron is thus localised. Notice that this type of localisation 
it different from other popular types. Anderson localisation in 2D, e.g., is due 
to quantum interferences of the electronic wave functions |34j . Here, however, 
the localisation is a purely classical effect. The high-field localisation is also 



3 In order to illustrate this point, consider a hiking tour in the mountains, e.g. around Les 
Houchcs in the French Alps. To go from one point to another one at the same height, one 
usual needs go downhill as well as uphill. It is very rare to be able to stay on the same height 
unless one wants to turn in circles that are just the closed countour lines which correspond to 
closed equipotcntial lines in our potential landscape. For those who participated at the Les 
Houches session which was outsourced to Singapore, where there are no mountains and where 
the whether is anyway too hot for hiking, just look at a hiking map of some mountainous 
region. Then search for countour lines that connect one border of the map to the opposite 
border. It turns out to be very hard to find such lines as compared to a large number of closed 
countour lines. 



Electronic Motion in an External Electrostatic Potential 



47 



different from the interaction-driven Mott insulator, where the electrons freeze 
out in order to minimise the mutual Coulomb repulsion between the electrons. 

At the edge, the equipotential lines reflect the confinement potential, which 
is zero in the bulk but rapidly increases when approaching the sample edge at 
y m in and y max (see Fig. 13. ip . In this case, the equipotential lines are open and 
therefore connect the two different electronic contacts. The electrons occupying 
quantum states at these equipotential lines then contribute to the electronic 
transport, in contrast to those on closed equipotential lines in the bulk. These 
quantum states are called extended sfcatesp as opposed to the localised states 
discussed above. The difference between localised and extended states turns 
out to be essential in the understanding of the IQHE, as we will see below (Sec. 



3.1.2 Electrostatic potential with translation invariance in 
the .x-direction 

Although the above semi-classical considerations yield the correct physical pic- 
ture of localised and extended states, it is based on the assumption that the 
electrostatic potential varies smoothly on the scale set by the magnetic length, 
such that we may replace the electron's position by that of its guiding centre [Eq. 
(|3.2[) ]. This assumption is, however, problematic in view of the confinement po- 
tential which varies strongly at the sample edges, i.e. in the vicinity of y m in and 
Umax- We will therefore treat the y-dependent confinement potential in a quan- 
tum treatment. Naturally, the appropriate gauge for the quantum- mechanical 
treatment is the Landau gauge (|2.3ip , which respects the translation invariance 
in the ^-direction, and the Hamiltonian (|2.55|) becomes 

H = ^ + \™u c {y - yo) 2 + V coni (y). (3.6) 
2m 2 

Remember that for a fixed wave vector k in the x-direction, the position around 
which the one-dimensional harmonic oscillator is centred is fixed by Eq. (|2.56l) , 
Ho = kl 2 B . We may therefore expand the confinement potential, even in the case 
of a strong variation, around this position, 



V{y) ~ V(y = kl B ) - eE(y )(y - y ) + O 



f d 2 V 

V d v 2 



where the local electric field is given in terms of the first derivative of the 
potential at yo, 

eE(yo) = - 



dy 

This expansion yields the Hamiltonian 

p 2 1 1 
H = - ]L + -mu) C (y - y' ) 2 + V con{ {y a ) - -mv D (y ), 
Am I 2 



*In the semi-classical picture the extended states are also called skipping orbits. 
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where the local drift velocity reads Vd = E(y )/B and the position of the 
harmonic oscillator is shifted, yo — ► y' = yo + eE(y ) jrauj 2 c . Notice that the 
last term is quadratic in the electric field E(y ) and therefore a second-order 
term in the expansion of the confinement potential. We neglect this term in the 
following calculations. The final Hamiltonian then reads 

H =^ + 2 muJc( y - Vo) + y co„f (2/o)> (3-7) 

where we have replaced the argument yo by the shifted harmonic-oscillator posi- 
tion y' , which is valid at first order in the expansion of the confinement potential. 
One therefore obtains the energy spectrum 

tn, V0 = huj c [n + + V{y ), (3.8) 

where we have omitted the prime at the shifted harmonic-oscillator position to 
simplify the notation. One therefore obtains the same LL spectrum as in the 
absence of a confinement potential, apart from an energy shift that is determined 
by the value of the confinement potential at the harmonic-oscillator position, 
which may indeed vary strongly. This position y$ plays the role of the guiding- 
centre position, as we have already mentioned in the last chapter, where we have 
calculated the electronic wave functions in the Landau gauge (|2.4.2p . One thus 
obtains a result that is consistent with the semi-classical treatment presented 
above. 



3.2 Conductance of a Single Landau Level 

We now calculate the conductance of a completely filled LL for the geometry 
depicted in Fig. 13.11 i.e. when all quantum states (described within the Landau 
gauge) of the n-th LL are occupied. In a first step, we calculate the current of 
the n-th LL, which flows from the left to the right contact, with the help of the 
formula [35] 

7 » =-^5Z<»,*l«*Kfc>. (3-9) 

k 

i.e. as the sum over all Nb quantum channels labelled by the wave vector 
k = 2Trm/L, with the velocity 

^ x \n,k) = -— = ——, 

in terms of the dispersion relation (|3.8p Fl Notice that the velocity in the in- 
direction is zero because the energy does not disperse as a function of the y- 

5 This relation may be obtained from the Heisenberg equations of motion, ihx = [x, H] = 
(dH/dp x )[x,p x ] = idH/dk, where we have used Eq. 112.101 1 and p x = hk. One therefore 
obtains the operator equation 

, _ 1 dH 
X ~ h~dk' 

which we evaluate in the state \n,k). In the last step we have used the periodic boundary 
conditions. 
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Figure 3.2: Edge states, (a) The LLs are bent upwards when approaching the 
sample edge, which may be modeled by an increasing confinement potential. 
One may associate with each LL n a maximal value y^ax °f the y-component 
where the LL crosses the chemical potential fJ, max . (b) At each position y 7 ^^, 
the filling factor decreases by a jump of 1. The n-th edge state is associated 
with the jump at y^ax an d the gradient of the confinement potential imposes a 
direction to the Hall drift of this state (chirality). This chirality is the same for 
all edge states at the same edge. 



component of the wave vector. The above expression is readily evaluated with 
Am = 1, and one obtains 

(n, k\v x \n, k) = ~ (e n m+1 - . 
h 

With the help of this expression, the current (|3.9p of the n-th LL becomes 

In ~J~ ^ ^ ~T (^n,m+l ^n,m) j 

m 

and one notices that all terms in the sum cancel apart from the boundary terms 
e n _m min arid ^n,m min i which correspond to the chemical potentials [i m in and 
fJ-max, respectively. The difference between these two chemical potentials may 
be described in terms of the (Hall) voltage V between the upper and the lower 
edge, Umax — Umin = —eV. One thus obtains the final result 

e e 2 

In = ~T (t^max- Vmin) = ~rV. (3.10) 

h h 

This means that each LL contributes one quantum of conductance G n = 
e 2 /h to the electronic transport and n completely filled LLs contribute a con- 
ductance^] 

n— 1 2 

G=J2 G n' =n-. (3.11) 

n'=Q 



6 Notice that, because the lowest LL is labelled by n' = 0, the last one has the index n — 1 
in the case of n completely filled levels. 
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Notice furthermore that this is a particular example of the Landauer-Biittiker 
formula of quantum transport 

e 2 

r — —T 

' n , J- n 

h 

through a conduction channel n, where T n is the transmission coefficient of the 
channel [36, 35, 37J. Because T n + R n = 1, in terms of the reflexion coefficient, 
the above result (|3 . 10[) indicates that each filled LL may be viewed as a conduc- 
tion channel with perfect transmission T n = 1, i.e. where an injected electron 
is not reflected or backscattered. 



3.2.1 Edge states 

The astonishing feature of perfect transmission, which is independent of the 
length L (or more precisely of the aspect ratio L/W, see the discussion in Sec. 
11.1.21 of the introduction) or the particular geometry of the sample, may be 
understood from the edge-state picture which we have introduced above (see 
Fig. 13. 2[) . Consider the upper edge, without loss of generality. The current- 
transporting edge state of the n-th LL is the one situated at y^ nax , where the n-th 
LL crosses the Fermi energy and where the filling factor jumps from v = n + 1 to 
7i0 Due to the upward bent of the confinement potential a particular direction is 
imposed on the electronic motion, which is nothing other than the Hall drift (sec 
Fig. 13. ip . This uni-dircctional motion is also called chirality of the edge state. 
Notice that this is the same chirality which one expects from the semi-classical 
expression (|3.5j) for the drift velocity. The chirality is the same for all edge states 
n at the same sample edge where the gradient of the confinement potential docs 
not change its direction. Therefore, even if an electron is scattered from one 
channel n to another one n' at the same edge it does not change its direction 
of motion, and the electron cannot be backscattered unless it is scattered to 
the opposite edge with inverse chirality. However, in a usual quantum Hall 
system, the opposite edges are separated by a macroscopic distance ~ W, and 
backscattering processes are therefore strongly (exponentially) suppressed in the 
ratio Ib/W between the magnetic length, which determines the spatial extension 
of quantum-mechanical state, and the macroscopic sample width W. Notice 
that the quantum Hall system at integer filling factors v = n is therefore a very 
unusual electron liquid: it is indeed a bulk insulator with perfectly conducting 
(non-dissipative) edges. 



7 Strictly speaking the filling factor does not jump not abruptly when one takes interactions 
between the electrons into account. In this case, two incompressible strips, of v = n + 1 and 
v = n are separated by a compressible strip of finite width. The picture of chiral electron 
transport remains, however, essentially the same when considering such compressible regions. 
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Figure 3.3: Two-terminal measurement. The current is driven through the 
sample via the left and the right contacts, where one also measures the voltage 
drop and thus a resistance. The upper edge is in thermodynamic equilibrium 
with the left contact (blue), whereas the lower one is in equilibrium with the 
right contact (red) . The chemical potential drops abruptly when the upper edge 
reaches the right contact, and when the lower edge reaches the left contact. Dis- 
sipation occurs in these hot spots (red dots). The measured resistance between 
the two contacts thus equals the Hall resistance. 



3.3 Two-terminal versus Six- Terminal Measure- 
ment 

3.3.1 Two-terminal measurement 

In the preceding section Sec. (|3.2p . we have calculated the conductance of 
a single LL (and n filled LLs) within a so-called two-terminal measurement, 
where we inject a current in the left contact with chemical potential (1l and 
collect the outcoming current at the right contact with u,r. As a consequence of 
Eq. (|3.10p , this current builds up a voltage V between the upper and the lower 
sample edge. This voltage drop is therefore associated with a Hall resistance, 
which is the inverse of the conductance G = ne 2 /h, 

i? ff = G- 1 = 4-' ( 3 - 12 ) 

e n 

and which coincides with the contact (or interface) resistance of a mcsoscopic 
system [37]. However, the voltage drop Vl between the left and the right contact 
is given by the difference of the chemical potentials in the contacts, /i/? — /it/, = 
—cVl, and the associated longitudinal resistance Vl/I is non-zero, in contrast 
to what we have seen in the introduction. This is due to the fact that the 
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difference between the longitudinal and the Hall resistance is not clearly defined 
in such a two-terminal measurement. 

This feature may be understood from Fig. 13.31 Indeed, due to the above- 
mentioned absence of backscattering, the chemical potential is constant along 
a sample edge, but there is a potential difference between the two edges. This 
means that the chemical potential must change somewhere along the edge. Con- 
sider the upper edge that is fed with electrons by the left contact, i.e. the upper 
edge is in thermodynamic equilibrium with the left contact and the chemical 
potentials therefore coincide, /il = [i max (see Fig l3.3j) . Now, when the upper 
edge touches the right contact which is at a different chemical potential [Ir, the 
chemical potential of the upper edge must rapidly relax to be in equilibrium 
with the right contact. In the same manner, the lower edge is in equilibrium 
with the right contact, /i m in = ^r, arid abruptly changes when touching the left 
contact. The rapid change in the chemical potential is associated with a dissi- 
pation of energy (at so-called hot spots) that has been observed experimentally 
|38| . In this experiment, the sample was put in liquid helium and the heating 
at the hot spots caused a local vaporisation of the helium observable in form of 
a fountain of gas bubbles. 

Due to the equivalence of the chemical potentials (1l = Umax and fJ, m i n = j-iR, 
the voltage drops V, between the upper and the lower edge, and Vl between 
the current contacts are equal, V — Vl- An unexpected consequence of this 
equation is that in a resistance measurement between the two contacts, in the 
two-terminal configuration, the two-terminal resistance equals the Hall resis- 
tance, 

Rr-l = Rh = -,— , (3.13) 

e" 2 n 

and not the (vanishing) longitudinal resistance, when the bulk is insulating (at 
v = n). 

3.3.2 Six-terminal measurement 

A more sophisticated geometry that allows for the simulaneous measurement 
of a well-defined longitudinal and Hall resistance is the six-terminal geometry, 
with two additional contacts at the upper and two at the lower edge [see Fig. 
I3.4f a)]. These additional contacts (2 and 3 at the upper and 5 and 6 at the lower 
edge, the left and the right contacts being labelled by 1 and 4, respectively) are 
used to measure a voltage, i.e. they have ideally an infinitely high internal 
resistance to prevent electrons to leak out of or into the sample. The chemical 
potential therefore remains constant at the upper edge /i^ = /i2 = as well 
as that at the lower edge fiR = /15 = fie, and one measures a zero- resistance, 
Rl = (fi2 — ^3) /si = (fJ-5 — fJ-e)/ e I = 0, as one expects from the calculation 
of the conductance through n LLs (see Sec. I3.11[) , which is entirely transverse. 
The conductance matrix is thus off-diagonal, as well as the resistance matrix, 




(3.14) 
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R H ~ h 5 -m^ r -^ l 

Figure 3.4: (a) Six-terminal measurement. The current I is driven through 
the sample via the contacts 1 and 4. Between these two contacts the chemical 
potential on the upper edge /i^ (blue) does not vary because the electrons do 
not leak out or in at the contacts 2 and 3, where one measures the longitudi- 
nal resistance. In the same manner, the chemical potential fj,R (red) remains 
constant between the contacts 5 and 6 on the lower edge. The longitudinal 
resistance measured between 2 and 3 as well as between 5 and 6 is therefore 
Rl = (fJ,2 — M3)/ e ^ = (^5 ~~ l l s)/ e I — 0. The Hall resistance is determined by 
the potential difference between the two edges and thus measured, e.g. between 
the contacts 5 and 3, where /i3 = /j_r — /iL, and thus Rh = (jU3 — (b) 
Four-terminal measurement in the van-der-Pauw geometry. In a Hall-resistance 
measurement, one drives a current through the sample via the contacts 1 and 
3 (connected by the continuous blue line) and measures the Hall resistance via 
the contacts 2 and 4 (dashed blue line). In a measurement of the longitudinal 
resistance, the current is driven through the sample via the contacts 1 and 4 
(continuous red line) and one measures a resistance between the contacts 2 and 
3 (connected by the dashed red line). 
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and one precisely measures the diagonal elements of the resistance matrix, the 
longitudinal resistance, between the contacs 3 and 2 (or 6 and 5). The off- 
diagonal elements, i.e. the Hall resistance, may e.g. be measured between the 
contacts 5 and 3 [as shown in Fig. 13.4( a)], and one measures then the result 
Rh = G^ 1 = h/e 2 n obtained from the calculation presented in Sec. 13.111 
because of the voltage drop V = ((II — Mi?)/e = (fi3 — ^z)/e between the upper 
and the lower edge. 

Finally wc mention that there exists an intermediate geometry that consists 
of four terminals (van-der-Pauw geometry) , where the resistance measurements 
are equally well defined [Fig. I3.4f b)]. If one labels the contacts from 1 to 4 in 
a clockwise manner, one may measure a Hall resistance between the contacts 
2 and 3 while driving a current through the sample by the contacts 1 and 3 
[blue lines in Fig. 13.41 Tb')] . In this case, one may use the clear topological 
definition mentioned at the beginning of the introduction. If one connects the 
contacts 2 and 3 by an imaginary line through the sample (dashed blue line) 
it necessarily crosses the imaginary line (continuous blue) which connects the 
current contacts 1 and 3 through the sample. This is precisely the topological 
definition of a Hall-resistance measurement. 

Similarly, one may measure the longitudinal resistance between the contacts 

2 and 3 if one drives a current through the sample via the contacts 1 and 4. 
In this case, the imaginary line (dashed red) which connects the contacts 2 and 

3 where one measures a resistance does not need to cross the line (continuous 
red) between the contacts 1 and 4 at which one injects and collects the current, 
respectively. As we have already mentioned at the beginning of the introduction, 
this defines topologically a measurement of the longitudinal resistance. 

These considerations show that a resistance measurement, although it docs 
not depend on the microscopic details of the sample, depends nevertheless on 
the geometry in which the contacts are placed at the sample [351 ISS]. This 
aspect is often not sufficiently appreciated in the literature, namely the fact 
that one measures, in a two-terminal geometry, a Hall resistance between the 
contacts that are used to inject and collect the current and not a longitudinal 
resistance, as one may have expected naively, when the system is in the IQHE 
condition. 

3.4 The Integer Quantum Hall Effect and Per- 
colation 

Until now we have shown that the Hall resistance is quantised [Eq. (|3. 12|) ] when 
n LLs are completely filled, i.e. when the filling factor is exactly v = n. However, 
we have not yet explained the occurence of plateaus in the Hall resistance, i.e. 
a Hall resistance that remains constant even if one varies the filling factor, 
e.g. by sweeping the magnetic field, around v = In order to explain the 



8 Strictly speaking, we have not gained anything because the quantum treatment allows us 
only to determine the Hall resistance at certain points of the Hall curve, those at the magnetic 
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Figure 3.5: Quantum Hall effect. The (impurity-broadened) density of states 
is shown in the first line for increasing fillings (a) - (c) described by the Fermi 
energy Ef- The second line represents the impurity-potential landscape the 
valleys of which become successively filled with electrons when increasing the 
filling factor, i.e. when lowering the magnetic field at fixed particle number. 
The third line shows the corresponding Hall (blue) and the longitudinal (red) 
resistance measured in a six-terminal geometry, as a function of the magnetic 
field. The first figure in column (c) indicates that the bulk extended states are 
in the centre of the DOS peaks, whereas the localised states are in the tails. 
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Constance of the Hall resistance over a rather large interval of magnetic field 
around v = n, we need to take into account the semi-classical localisation of 
additional electrons (or holes) described in Sec. 13.11 This is shown in Fig. 13. 5i 
where we represent the filling of the LLs (first line), the potential landscape of 
the last partially-filled level (second line) and the resistances as a function of 
the magnetic field, measured in a six-terminal geometry (third line). We start 
with the situation of n completely filled LLs [column (a) of Fig. 13. 5j . which 
we have extensively discussed above: the LL n (and its potential landscape) is 
unoccupiedH In a six-terminal measurement, one thcrfore measures the Hall 
resistance Rh = h/e 2 n and a zero longitudinal resistance, as we have seen in 
Eq. IpOijl . 

In column (b) of Fig. 13.51 we represent the situation where the LL n gets 
moderately filled by electrons when the magnetic field B is decreased. These 
electrons in n populate preferentially the valleys of the potential landscape, or 
more precisely the closed equipotential lines that enclose these valleys. The 
electrons in the LL n are thus (classically) localised somewhere in the bulk and 
do not affect the global transport characteristics, measured by the resistances, 
because they are not probed by the sample contacts. Therefore, the Hall re- 
sistance remains unaltered and the longitudinal resistance remains zero despite 
the change of the magnetic field. This is the origin of the plateau in the Hall 
resistance. 

If one continues to lower the magnetic field, the regions of the potential land- 
scape in the LL n occupied by electrons become larger, and they arc eventually 
enclosed by equipotential lines that pass through the bulk and that connect the 
opposite edges. In this case, an electron injected at the left contact and trav- 
elling a certain distance at the upper edge may jump into the state associated 
with this equipotential line and thus reach the lower edge. Due to its chirality, 
the electron is then backscattered to the left contact, which causes an increase 
in the longitudinal resistance. Indeed, if one measures the resistance between 
the two contacts at the lower edge, a potential drop is caused by the electron 
that leaks in from this equipotential connecting the upper and the lower edge. 
It is this potential drop that causes a non-zero longitudinal resistance. At the 
same moment the Hall resistance is no longer quantised and jumps to the next 
(lower) plateau, a situation that is called plateau transition. This situation of 
electron-filled equipotential lines connecting opposite edges, which are thus ex- 
tended states [see first line of Fig. 13.5( c)] as opposed to the bulk localised states, 
arises when the LL n is approximately half-filled. Notice that these extended 
states, which are found in the centre of the DOS peaks [see upper part of Fig. 
13.5( c)]. are bulk states in contrast to the above-mentioned edge states, which 
are naturally also extended 

The clean jump in the Hall resistance at the plateau transition accompanied 

fields corresponding to v = hn e i/eB = n. If we substitute the filling factor in Eq. 113.121 1. we 
see immediately that Rh = h/e 2 u = B/en e i, i.e. one retrieves the classical result for the Hall 
resistance. 

^Remember that due to the label for the lowest LL, all LLs with v! = 0, n— 1 are the 
completely filled and the LL n is then the lowest unoccupied level. 
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Figure 3.6: STS measurements by Hashimoto ct al., 2008, on a 2D electron 
system on a n-InSb surface. The figures (a) - (g) show the local DOS at various 
sample voltages, around the peak obtained from a dl/dV measurement (h). 
Figure (i) shows a calculated characteristic LDOS, and figure (j) an STS result 
on a larger scale. 

by a peak in the longitudinal one is only visible in the six- (or four-)terminal 
measurement. As we have argued in Sec. 13.3.21 there is no clear cut between the 
longitudinal and the Hall resistivity in the two-terminal configuration, where 
the resistance measured between the current contacts is indeed quantised in 
the IQHE. At the plateau transition, however, the chemical potential at the 
edges is no longer constant because of backscattered electrons and the resistance 
is no longer quantised. One observes indeed the resistance peak associated 
with the longitudinal resistance in the six- or four-terminal configuration. As a 
consequence, one measures, at the plateau transition, the superposition of the 
Hall and the longitudinal resistances. 

If one increases even more the filling of the LL n, the same arguments apply 
but now in terms of hole states. The Hall resistance is quantised as Rh = 
h/e 2 (n + 1), and the holes (i.e. the lacking electrons with respect to n + 1 
completely filled LLs) get localised in states at closed equipotential lines around 
the potential summits. As a consequence, the longitudinal resistance drops to 
zero again. 

3.4.1 Extended and localised bulk states in an optical mea- 
surement 

The physical picture presented above, in terms of localised and extended bulk 
states, has recently been confirmed in scanning-tunneling spectroscopy (STS) 
of a 2D electron system that was prepared on an n-InSb surface instead of the 
more common GaAs/AlGaAs heterostructure [50]. Its advantage consists of its 
accessibility by an "optical" (surface) measurement that cannot be performed 
if the 2D electron gas is buried deep in a semiconductor heterostructure. In an 
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STS measurement one scans the sample and thus measures the local density of 
states at a certain energy that can be tuned via the voltage between the tip 
of the electron microscope and the sample. When measuring the differential 
conductance dl/dV, which is proportional to the DOS, one observes a peak 
that corresponds to the centre of a LL [Fig. 13.6( h)] where the extended states 
are capable of transporting a current between the different electric contacts, as 
mentioned above. Whereas the quantum states at energies corresponding to 
closed cquipotential lines of the impurity landscape are clearly visible as closed 
orbits in Fig. I3.6( a).(b) and (f),(g), the states in the vicinity of the peak are 
more and more extended, as shown by the spaghetti-like lines in Figs. I3.6f c).(d) 
and (e), as one would expect from the arguments presented above. 

3.4.2 Plateau transitions and scaling laws 

The physical picture presented above suggests that the plateau transition in the 
Hall resistance is related to a percolation transition, where initially separated 
electron-filled valleys start to percolate between the opposite sample edges be- 
yond a certain threshold of the filling. Because of the second-order character 
of a percolation transition, this scenario suggests that the plateau transition is 
a second-order quantum phase transition described by universal scaling laws, 
where the control parameter is just the magnetic field B [?TJ[12]. We finish this 
chapter on the IQHE with a brief overview over these scaling laws, and refer the 
interested reader to the literature [41] [42] an d the class given by G. Batrouni at 
the same Singapore session of Les Houches Summer School 20090 

The phase transition occurs at the critical magnetic field B c and is charac- 
terised by an algebraically diverging correlation length 

t~\B-B c \-\ (3.15) 

where v is called the critical exponent^ In the same manner, the temporal 
fluctuations are described by a correlation "length" £ r that is related to the 
spatial correlation length £, 

£ T ~£*~ \B-B c \-*\ (3.16) 

where z is called dynamical critical exponent. It is roughly a measure of the 
anisotropy between the spatial and temporal fluctuations, which is often en- 
countered in non-relativistic condensed-matter systems P^l 

At the phase transition B c , the longitudinal and transverse resistivities Pl/h 
arc described in terms of universal functions that are functions of the ratio t/£ t 

10 The lecture notes for this class are availabcl on the School's program webpage: 
| http://www.ntuTe du.sg/ias/upcorningcvcnts/LHSOPS09/Pagcs/progr amme.aspxl 

11 Although we use the same Greek letter v for the critical exponent, it must not be con- 
funded with the filling factor, which plays no role in this subsection. 

12 Notice that in relativity, time is considered as the "fourth" dimension, and Lorentz in- 
variancc would require that spatial and temporal fluctuations be equivalent, i.e. z = 1. 
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Figure 3.7: Experiment by Wei et a/., 1988. The width of the transition AB 
and of the derivative of the Hall resisitivity dp xy /dB, measured as a function 
of temperature, reveals a scaling law with an exponent 1/zv = 0.42 ± 0.04, for 
the transition between the filling factors 1 — ► 2 (TV = |), 2^3 (TV = 1 f) and 
3^4 (TV = 11). 



between the (imaginary) time r, which is proportional to the inverse tempera- 
ture, Ti/t = ksT [1TJ 52] and the temporal correlation length £ r , 

Pl/h = Jl/h (yj-^j = Jl/h , (3-17) 

where we have defined AB = \B — B c \. In the case of an AC (alternating 
current) measurement at frequency u, another dimensionlcss quantity, namely 
the ratio between the frequency and the temperature, tko/ksT, needs to be 
taken into account such that the universal function reads 

AC _ f ( r hu 

However, we do not consider an alternating current here. Equation (|3.17[) then 
yields the scaling of the width of the peak in the longitudinal resistance (or else 
the plateau transition) 

AB-T 1 '™. (3.18) 

A measurement of this width by Wei et al. [13] has confirmed such critical 
behaviour with an exponent 1/zv = 0.42 ± 0.04 (see Fig. 13. 7[) . 

Furthermore, one may distinguish between the two exponents v and z within 
a measurement of the plateau-transition width as a function of the electric field 
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E via current fluctuations. One may identify the energy fluctuation eE£ at 
the correlation length £ with the energy scale 7i/£ T ~ ft/t; z set by the temporal 
fluctuation £ T) which yields E ~ £-( 1+z ) ~ AB^ 1+Z \ and thus 

AB~B 1 /Ki+«). (3.19) 

Other measurements by Wei et al. [44] have shown that these types of fluctu- 
ations yield z — 1, i.e. ^ ~ 2.3. The precision of the measured critical expo- 
nent has since been improved - more recent experiments [45 1 146j have revealed 
v = 2.38 ±0.06. 

Theoretically one knows that the critical exponent for classical 2D perco- 
lation is fciass = 4/3 and thus much smaller than the measured one. This 
discrepancy is due to the quantum nature of the percolation in quantum Hall 
systems. Indeed, quantum- mechanical tunneling and the typical extension ~ Is 
of the wave functions associated with the equipotcntial lines enhance percola- 
tion, i.e. the electron puddles in the potential valleys may percolate before they 
touch each other in the classical sense. A model that takes into account this 
effect has been proposed by Chalkcr and Coddington [47] , though with simplify- 
ing assumptions for the puddle geometry^ and one obtains a critical exponent 
v = 2.5 ± 0.5 from numerical studies of this model [471 EH], in quite a good 
agreement with the experimental data [131 EI] • 

In spite of the good agreement with experimental findings, these theoreti- 
cal results need to be handled with care - indeed, analytical calculations have 
shown that the dynamical exponent should be exactly z = 2 for non-interacting 
electrons, whereas the measured value z ~ 1 is obtained when interactions are 
taken into account on the level of the Hartrec-Fock approximation [J5]. Fur- 
thermore, very recent numerical calculations within the Chalker-Coddington 
model have shown that the accurate value of the critical exponent is slightly 
larger (y ~ 2.59) than the measured one when interactions are not taken into 
account [50] . 

3.5 Relativistic Quantum Hall Effect in Graphene 

We finish this chapter on the IQHE with a short presentation of the relativis- 
tic quantum Hall effect (RQHE) in graphene, which is understandable in the 
same framework of LL quantisation and (semi-classical) one-particle localisa- 
tion as the IQHE in a non-relativistic 2D electron system. Indeed, the above 
arguments also apply to relativistic electrons in graphene, but we need to take 
into account the two different carrier types, electrons and holes, which carry a 
different charge. This is not so much a problem in the case of the impurity po- 
tential with its valleys and summits: in a particle-hole transformation, a valley 
becomes a summit and vice versa Furthermore, the direction of the Hall drift 

13 Notice, however, that due to the universality of the scaling laws and the fluctuations at all 
length scales, the results are expected to be independent on these microscopic assumptions. 

14 The particle-hole transformed landscape corresponds to an impurity distribution in which 
one interchanges negatively and positively charged impurities. 
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Figure 3.8: (a) Mass confinement for relativistic Landau levels. Whereas the 
electron-like LLs (A = +) arc bent upwards when approaching the sample edge 
(Umax), the hole-like LLs (A = — ) are bent downwards. The fate of the n = 
LL depends on the valley (parity anomaly) - in one valley (K), the level energy 
decreases, whereas it increases in the other valley (K'). (b) Filling of the bulk 
Landau levels at v = 0. All electron-like LLs (A = +) are unoccupied whereas 
all hole-like LLs (A = — ) are completely filled. The n = LL is altogether 



changes in this transformation. Because of the universality of the quantum Hall 
effect, both types of impurity distributions related by particle-hole symmetry 
yield the same quantisation of the Hall resistance. The picture of semi-classical 
localisation therefore applies also in the case of relativistic electrons in graphene. 

The situation is different for the confinement potential. An ansatz of the 
form V(y)t - remember that the Hamiltonian of electrons in relativistic graphene 
is a 2 x 2 matrix that reflects the two different sublattices A and B - has 
the problem that an increase V{y — y m ax/min) — * oo at the sample edge con- 
fines electrons but not the holes of the valence band for which we would need 
V(y — Dmax/min) ~ ¥ ~ 00 f° r an efficient confinement. A possible confinement 
potential may be formed with the Pauli matrix a z , 



which, together with the Hamiltonian (|2.8|) . yields the Hamiltonian which cor- 
responds to the non-relativistic model (|3.6j) . For a constant term M = V(y) 
the contribution (|3.20|) plays the role of a mass of a relativistic particle (see also 
Appendix |B|) . Therefore, the confinement (|3.20|) is sometimes also called mass 
confinement. The corresponding energy spectrum, which one obtains within the 
same approximation as in Sec. 13.11 via the replacement y — > y = kl 2 B in the 
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Landau gauge, reads [c.f. Eq. (|B.8[) in Appendix [B] 



exn, yo = AWM 2 (y ) + 2 — n > (3-21) 

and is schematically represented in Fig. I3.8f a). Notice that Eq. (|3.21j) is only 
valid for n ^= indeed, the n = acquires a non-zero energy M(j/o), which is 
negative for our particular choice (see Appendix [Bj . This feature is sometimes 
called parity anomaly in high-energy physics. Remember that in the case of 
graphene, one has two inequivalent low-energy points in the first BZ which give 
rise to a relativistic energy spectrum. The Dirac Hamiltonians (|2.4p and (|2.8[) 
for the zero-B and magnetic-field case, respectively, applies principally only to 
one of the two valleys (say K), whereas that for the other valley is given by 
— Hd (or —Hp) if one interchanges the A and B components [c.f. Eq. (|A.16|) 
in Appendix |A] . The confinement term <|3.20[) therefore reads — V CO ni(y) in the 
other valley, i.e. with a negative mass. The n = LL thus shifts to positive 
energies in the second valley, and the two-fold valley degeneracy is lifted in this 
level. A more detailed discussion of the mass confinement (|3.20[) in graphene 
may be found in the Appendix [Bl 

This type of confinement may seem to be somewhat artificial, whereas the 
confinement in the non-relativistic case is easier to accept. Notice, however, that 
the whole model of massless Dirac fermions [second Hamiltonian in Eq. (|2 .4|) ] 
only describes the physical properties at length scales that are large compared 
to the lattice spacing (in graphene). In the true lattice model, the electrons are 
naturally confined because one does not allow for hopping from a lattice site 
at the edge into free space. The expression ()3.20[) is therefore only an effective 
model to describe confinement. We notice that, although the effective model 
yields a qualitatively correct picture, the fine structure of the dispersion at the 
edge depends on the edge geometry |51j . For further reading, we refer the 
interested reader to the literature [21]. 

With the help of these preliminary considerations, we are now prepared to 
understand the RQHE in graphene - the semi-classical localisation is the same 
as in the non-relativistic case, and the confinement, which needed to be adopted 
to account for the simultaneous presence of electron- and hole-like LLs, yields 
the edge states which are responsible for the quantum transport and, thus, the 
resistance quantisation. The RQHE was indeed discovered in 2005 by two differ- 
ent groups OB|, and the results are shown in Fig. 13.91 [7], The phenomenology 
of the RQHE is the same as that of the IQHE in non-relativistic LLs: one ob- 
serves plateaus in the Hall resistance while the longitudinal resistance vanishes. 
Notice that one may vary the filling factor either by changing the -B-field at 
fixed carrier density [Fig. 13.9( a)] or one keeps the i?-field fixed while changing 
the carrier density with the help of a gate voltage [Fig. I3.9f c)]. The latter 
measurement is much easier to perform in graphene than in non-relativistic 2D 
electron gases in semiconductor heterostructures. 

In spite of the similarity with the non-relativistic IQHE, one notices, in Fig. 
13.91 an essential difference: the quantum Hall effect is observed at the filling 
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Figure 3.9: Measurement of the relativistic quantum Hall effect (Zhang efc al., 
2005). (a) RQHE at fixed carrier density (V G = 15 V) at T = 30 mK. The 
filling factor is varied by sweeping the magnetic field, (b) Sketch of the DOS 
with the Fermi energy between the LLs n = and +, n = 1. (c) RQHE at fixed 
magnetic field (B = 9 T) at higher temperatures, T = 1.6 K. The filling factor 
is now varied by changing the gate voltage. 
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factors 

i/ = ±2(2n+l), (3.22) 

in terms of the LL quantum number n, whereas the IQHE is observed at v = n 
(or v = 2n if the LLs are spin-degenerate). The step in units of 4 is easy to 
understand: each relativistic LL in graphene is four-fold degenerate (in addition 
to the guiding-centre degeneracy), due to the two-fold spin and the additional 
two-fold valley degeneracy. However, there is an "offset" of 2. This is due to the 
fact that the filling factor v = corresponds to no carriers in the system, i.e. 
to a situation where the Fermi energy is exactly at the Dirac point (undoped 
graphene). In this case, one has a perfect electron- hole symmetry, and the 
n = LL must therefore be half-filled [see Fig. I3.8f b)]. or else: there are as 
many electrons as holes in n = 0. According to the considerations presented in 
Sec. 13.41 this does not correspond to a situation where one observes a quantum 
Hall effect due to percolating extended states. Indeed, the system turns out 
to be metallic at v = with a finite non-zero longitudinal resistance [6J [7] . A 
situation, where one would expect a quantum Hall effect, arises when the central 
LL n = is completely filled (or completely empty). As a consequence of the 
four-fold level degeneracy, one obtains the quantum Hall effect at v = 2 (or 
;/ = —2) observed in the experiments (see Fig. I3.9|) . This is the origin of the 
particular filling-factor sequence (|3.22|) of the RQHE in graphene. 



Chapter 4 



Strong Correlations and the 
Fractional Quantum Hall 
Effect 

In the preceding chapter, we have seen that one may understand the essential 
featues of the IQHE within a one-particle picture, i.e. in terms of Landau 
quantisation; at integer filling factors v = n, which correspond to n completely 
filled LLsQ an additional electron is forced, as a result of the Pauli principle, 
to populate the next higher (unoccupied) LL [see Fig. 14.1( a)]. It therefore, 
needs to "pay" a finite amount of energy hue [or y/2(hv /lB){y/n — \Jn — 1) in 
the case of the RQHE in graphene] and is localised by the impurities in the 
sample, due to the classical Hall drift which forces the electron to move on 
closed cquipotential lines. The system is said to be incompressible because one 
may not vary the filling factor and pay only an infinitesimal amount of energy - 
indeed in the case of a fixed particle number, consider an infinitesimal decrease 
of the magnetic field which amounts to an infinitesimal change of the surface 
2ttI b occupied by each quantum state. Since the total surface of the system 
remains constant, the infinitesimal increase of 2irl B may not be accomodated 
by an infinitesimal change in energy, due to the gap between the LL n — 1 and 
n where at least one electron must be promoted to. This gives rise to a zero 
compressibility. 

In view of this picture of the quantum Hall effect, it was therefore a big 
surprise to observe a FQHE at a filling factor v = 1/3, with the corresponding 
Hall quantisation Rh = h/e 2 v = 3/i/e 2 [13], and, later, at a large set of other 
fractional filling factors. Indeed, if only the kinetic energy is taken into account, 
the ground state at v = 1/3 is highly degenerate and there is no evident gap 
present in the system: the Pauli principle no longer prevents an additional 
electron to populate the next higher LL, but it finds enough place in the lowest 

As before, we neglect the electron spin to render the discussion as simple as possible. The 
role of spin will be discussed briefly in the last chapter on multi-component systems. 
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(a) 



(b) 
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Figure 4.1: (a) Sketch of a completely occupied LL. An additional electron (grey 
circle) is forced to populate the next higher LL because of the Pauli principle, (b ) 
Sketch of a partially filled LL. Because of the presence of unoccupied states in the 
LL (crosses), the Pauli principle does not prevent an additional electron (grey 
circle) to populate the next higher LL. The low-energy dynamical properties of 
the electrons are described by excitations within the same LL (no cost in kinetic 
energy), and inter-LL excitations are now part of the high-energy degrees of 
freedom. 

LL which is only one-third filled. 

Notice that we have neglected so far the mutual Coulomb repulsion between 
the electrons, which happens to be responsible for the occurence of the FHQE. 
The relevance of electronic interactions is discussed in the next section (Sec. 
14. ip . In Sec. 14. 2\ we present the basic results of Laughlin's theory of the 
FQHE, such as the ground-state wave functions, fractionally charged quasi- 
particles and the interpretation of Laughlin's wave function in terms of a 2D 
one-component plasma. The related issue of fractional statistics is introduced 
in a section apart (Sec. 14. 3p . and we close this chapter with a short discussion 
of different generalisations of Laughlin's wave function, such as CF theory or 
the Moore- Read wave function in half-filled LLs. 

4.1 The Role of Coulomb Interactions 

As already mentioned above, the situation of a partially filled LL is somewhat 
opposite to that of n completely occupied levels, where one observes the IQHE. 
This difference is summarised in Fig. 14. II and it is also the origin of the different 
role played by the Coulomb repulsion between the electrons. In the case of n 
completely filled LLs, one has a non-degenerate (Fermi-liquid-like) ground state, 
where the interactions may be treated within a perturbative approach. Indeed, 
any type of excitation involves a transition between two adjacent LLs that are 
separated by an energy gap of hue [see Fig. l4.lT a)]Fl and we need to compare 
the Coulomb energy at the characteristic length scale Rc = IsV^n + 1 to this 



2 In order to simplify the discussion, we consider only the IQHE in non-relativistic quantum 
Hall systems, but the arguments apply also to the RQHE in graphene. 
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which turns out to be nothing other than the usual dimensionless coupling 
constant 

me 2 _i 

for the 2D Coulomb gas in Fermi liquid theory [SH [53j ■ The last expression is 
obtained by identifying the Fermi energy Ep = K k F /2m, in terms of the Fermi 
wave vector fcp, with the energy of the last occupied LL fiucn. The perturba- 
tive approach allows one, e.g., to describe collective electronic excitations in the 
IQHE, such as magncto-plasmon modes (the 2D plasmon in the presence of a 
magnetic field) or magncto-excitons (intcr-LL excitations that acquire a disper- 
sion due to the Coulomb interaction) [M], or else the corresponding modes of 
the RQHE in graphene [531 [56]. 

In the case of a partially filled LL n the situation is inverted: for an electronic 
excitation, there are enough unoccupied states in the LL n for an electron of the 
same level to hop to. From the point of view of the kinetic energy, there is no 
energy cost associated with such an excitation (low-energy degrees of freedom) 
whereas an excitation to the next higher (unoccupied) LL costs an energy Tiujc- 
Inter-LL excitations may then be neglected as belonging to high-energy degrees 
of freedom [Fig. 14.1T b)]. Notice that all possible distributions of N electrons 
within the same partially filled LL n therefore have the same kinetic energy, 
which effectively drops out of the problem. The macroscopic degeneracy may 
be lifted by phenomena due to other energy scales, such as those associated with 
the impurities in the sample or else the electron-electron interactions. The first 
hypothesis (impurities) may be immediately discarded as the driving mechanism 
of the FQHE because, in contrast to the IQHE, the FQHE only occurs in high- 
quality samples with low impurity concentrations. Indeed, the hierarchy of 
energy scales in the FQHE may be characterised by the succession 

hue <; v c > v imp , (4.1) 

and we therefore need to consider seriously the Coulomb repulsion, which govern 
the low-energy electronic properties in a partially filled LL|f| Notice that we thus 
obtain a system of strongly-correlated electrons for the description of which all 
perturbative approaches starting from the Fermi liquid are doomed to fail. The 
only hope one may have to describe the FQHE is then a well-educated guess of 
the ground state. 

The most natural guess would be that the electrons in a partially filled LL 
behave as classical charged particles that form a crystalline state in order to 
minimise their mutual Coulomb repulsion. Such a state is also called Wigncr 
crystal (WC) because it was first proposed by Wigner in 1934 [57] . A WC has 
indeed been thought - before the discovery of the FQHE - to be the ground state 
of electrons in a partially filled LL [55]. Even if the WC is the ground state at 
very low filling factors, as it has been shown experimentally [59] , this state may 

3 As for the IQHE, impurities play nevertheless an important role in the localisation of 
quasi-particlcs, which we need to invoke later in this chapter in order to explain the transport 
properties of the FQHE. 
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not allow for an explanation of the FQHE. Indeed, the WC is a state that breaks 
a continuous spatial symmetry (translation invariance) and any such state has 
gapless long- wave- length excitations (Goldstone modes). The Goldstone mode 
of the WC (as of any other crystal) is the acoustic phonon the energy of which 
tends to zero at zero wave vector. One may thus compress the WC by changing 
the occupied surface in an infinitesimal manner or else by adding an electron 
without changing the macroscopic surface and pay only an infinitesimal amount 
of energy. The ground state is therefore compressible, i.e. it is not separated 
by an energy gap from its single-particle excitations, a situation that is at odds 
with the FQHE. 

4.2 Laughlin's Theory 

As a consequence of the above-mentioned considerations on the WC, one thus 
needs to search for a candidate for the ground state that does not break any 
continuous spatial symmetry and that has an energy gap. Such a state is the 
incompressible quantum liquid which was proposed by Laughlin in 1983 [14j 
the basic features of which we present in the present section. We consider, 
here, only the FQHE in the lowest LL (LLL), for simplicity. There are different 
prescriptions to generalise the associated wave functions to higher LLs, e.g. with 
the help of Eq. (|2~T8j) (see MacDonald, 1984). Experimentally, several FQHE 
states have been observed in the next higher LL n = 1 although the majority 
of FQHE states is found in the LLL0 

4.2.1 Laughlin's guess from two-particle wave functions 

In order to illustrate - one cannot speak of a derivation - Laughlin's wave func- 
tion, we first need to remember the one-particle wave function of the LLL and 
then consider the corresponding two-particle wave function. We have already 
seen in Sec. 12.4.11 that a one-particle wave function in the LLL is described in 
terms of an analytic function times a Gaussian[f| 

V> ~ z m 'e-l*l 2 / 4 , 

in terms of the integer m' = 0, ...,Nb — 1, where we have absorbed now (and 
in the remainder of these lecture notes) the magnetic length in the definition of 
the complex position, z = (x — iy)/lB- 

Consider, in a second step, an arbitrary two-particle wave function. This 
wave function must also be an analytic function of both postions z\ and Zi 
of the first and second particle, respectively, and may be a superposition of 
polynomials, such as e.g. of the basis states 

^ 2 \z, Z) ~ Z M z m e~^ 2+ ^/\ (4.2) 

4 There is even some slight indication for a 1/5 FQHE state in the next excited LL n = 2 

5 We neglect the numerical prefactors here that account for the normalisation of the wave 
functions. 
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where we have defined the centre of mass coordinate Z = (zi + Z2)/2 and the 
relative coordinate z = (z\ — z^). The quantum number m plays the role of 
the relative angular momentum between the two particles, and M is associated 
with the total angular momentum of the pair. Because of the analyticity of the 
LLL wave functions, m must be an integer, and the exchange of the positions 
z\ and Z2 imposes on m to be odd because of the electrons' fermionic nature. 

Laughlin's wave function [T3] is a straight-forward TV-particle generalisation 
of the two-particle wave function (|4.2p . 

^({ Zj ,z*}) =H(z k -z i re-Z i \^, (4.3) 

k<l 

where we have omitted the normalisation constants in order to simplify the 
notation and where all indices run from 1 to the total number of particles N . 
Notice that there is no dependence on the centre of mass, but only on the relative 
coordinates between the particle pairs. Had there been such a dependence, 
described by a non-zero value of the total angular momentum quantum number 
M 7^ 0, one would have broken a continuous spatial symmetry, in which case the 
state would describe a compressible rather than an incompressible state required 
for the FQHE, as we have mentioned above. We emphasize once again that 
Laughlin's wave function is not based on a mathematical derivation, although we 
will see below that there exist some mathematical models for which it describes 
the exact ground state, but it is more appropriately characterised as a variational 
wave function. 

Variational parameter 

The variational parameter in Laughlin's wave function (|4.3|) is nothing other 
than the exponent m, with respect to which we would, in principle, need to 
optimised the wave function in order to approximate the true ground state of 
the system. Notice, however, that due to the LLL analyticity condition and 
fermionic statistics, the exponent is restricted to odd integers, m = 2s + 1, in 
terms of the integer s. Furthermore, this variational parameter turns out to 
be fully determined by the filling factor v, as we will show with the following 
argument @ 

Consider Laughlin's wave function as a function of the position Zk of some 
arbitrary but fixed electron k. There are N — 1 factors of the type (z k — zi) m , 
one for each of the remaining N — 1 electrons, I, occuring in the ansatz (|4.3|) . 
such that the highest power of Zk is m(N — 1), 

[z h -zi) ~ z k y 

k<l 

Now, remember from Sec. 12.4.11 [see Eq. (|2.54j) ] that the highest power of the 
complex particle position is fixed by the number of states Nb in each LL. This 



6 We are therefore confronted with the somewhat bizarre situation where we dispose of a 
variational wave function with no possible variation. 
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yields the relation 

mN -5 = N B (4.4) 

between the number of particles N and the number of flux quanta Nb threading 
the system. Here, S is some shift that is on the order of unity and that plays 
no role in the thermodynamic limit N, Nb — * oo0 Because the ratio between 
the number of particles and that of flux quanta is nothing other than the LL 
filling factor (|2.44|1 . v = N/Nb, one notices that, in the thermodynamic limit, 
the "variational parameter" is entirely fixed by the filling factor, i.e. 



m = 2s + 1 = 



1 1 
v = — = . 

m 2s +1 



(4.5) 



and Laughlin's wave function is therefore a candidate wave function for the 
ground state at the filling factors 

v = 1, 1/3, 1/5,... 

Remember that the odd value m = 2s + 1 is required by the fermionic nature 
of the electrons. Formally, one may though lift this restriction and generalise 
Laughlin's wave function to bosonic particles by choosing an even exponent 2s. 
Such bosonic Laughlin wave functions have been studied theoretically in the 
context of rotating cold Bose gases in an optical trap [61] . 



Laughlin's wave function at v = 1 

It may seem, at first sight, astonishing that also the case of a completely filled 
LL for v = 1 is described in terms of a Laughlin wave function with m = 1 (or 
s = 0). Indeed, the state 



^({z 3 })=f N ({z 3 })e 



7-i 



would be non-degenerate and could thus be described in terms of a Slater de- 
terminant, 



fN({zj}) = det 



V 



N-1 \ 



N-1 
z 2 



.N-1 



(4.6) 



-N 



'AT 



-N 



where we have omitted the ubiquitous Gaussian factor exp(— \zj\ 2 /4). No- 
tice that the j-th line in this determinant corresponds to all LLL states of the 
j-th particle described in terms of the polynomials zj 1 . The determinant takes 
into account all permutations of the N particles over the TV particle positions, 
z\, ...,zn, and may be rewritten in a compact manner with the help of the 
co-called Vandcrmondc determinant, 



lN{{Zj}) = Y[(Zi~ Zj), 



(4.7) 



i<j 



7 Notice, however, that this shift plays an important role in numerical calculations, such as 
exact diagonalisation, when performed on special geometries, such as on a sphere. 
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which is indeed nothing other than the polynimial prefactor in Laughlin's wave 
function (|4.3[) with m = 1. 

Until now we have obtained an iV-particle wave function from some very 
general symmetry considerations (LLL analyticity condition, fermionic statis- 
tics, no broken continuous spatial symmetries), but we have not at all shown 
that it describes indeed the ground state responsible of the FQHE. In the fol- 
lowing parts, we will therefore discuss the basic physical properties of this, for 
the moment rather abstract, mathematical entity. In a first step, we will discuss 
some energy properties of the ground state and show that Laughlin's wave func- 
tion is the exact ground state of a certain class of models that are qualitatively 
compared to the physical one (Coulomb interaction). We will then discuss the 
fractionally charged quasi-particle excitations of this wave function. 



4.2.2 Haldane's pseudopotentials 

In order to describe the energetic properties of Laughlin's wave function 
we consider again the two-particle wave function (|4.2p . Notice that this wave 
function is an exact eigenstate for any central interaction potential that depends 
only on the relative coordinate z between particle pairs, such as it is the case 
for the Coulomb interaction, V = V(|z|). One may therefore decompose the 
interaction potential in the relative angular momentum quantum numbers m, 

„ _ (m,M\V\m,M) 



(to, M\m, M) 



where the denominator takes into account the fact that we have not properly 
normalised the two-particle wave functions O, ^ 2 \z,Z) = {z,Z\m,M)ft 
The fact that there is no dependence on M is a direct consequence of the as- 
sumption that we deal with a central interaction potential, i.e. (z, Z\V\z', Z') = 
V(\z\)5 z , z >Sz.z' ■ Furthermore, there are no off-diagonal terms of the form (to, M|T^|to', M), 
with to' 7^ to, as one may show explicitly in the polar representation z = 
pexp(^), 

(m,M\V\m',M) cx / d(j) dp p m+m ' +1 V{p)e- t{m - m '^ cx 6 m , m >, 



due to the integration over the polar angle. The potentials v m obtained from the 
decomposition into relative angular momentum states are also called Haldane's 
pseudopotentials [62] . They fully characterise the two-particle energy spectrum 
because the kinetic energy is the same for all two-particle states \m,M), as de- 
scribed above. Notice that this is a very special case: normally any repulsive 
interaction potential yields unbound states with a continuous energy spectrum, 
such as the plane-wave states in scattering theory. Here, the energy spectrum is 
discrete even if the interaction is repulsive, due to the presence of a quantising 
magnetic field. Notice further that Haldane's pseudopotentials are an image of 



8 In order to simplify the notations, we have omitted the LL quantum number n = 0, which 
is the same for both particles in this wave function. 
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Relative Angular Momentum 



10 



Figure 4.2: Haldanc's pseudopotcntials for the Coulomb interaction in the LLs 
n = and n = 1. Notice that we have plotted the pseudopotentials for both 
odd and even values of the relative angular momentum m even though only odd 
values matter in the case of fcrmions. 



the real-space form of the interaction potential. Indeed, if a pair of electrons is 
in a quantum state with relative angular momentum m, the average distance be- 
tween the electrons is \z\ ~ Zs\/2mO Haldanc's pseudopotcntial v m is therefore 
roughly the value of the original interaction potential at the relative distance 

v m ^V (\z\ = l B V2^) , (4.9) 

and the small-m components of Haldanc's pseudopotcntials correspond to the 
short-range components of the underlying interaction potential. Figure 14.21 
shows the pseudopotcntial expansion for the Coulomb interaction in the lowest 
(n = 0) and the first excited (n = 1) LL. 

Haldane's pseudopotentials are extremely useful in the description of the 
iV-particle state as well. Indeed, the iV-particle interaction Hamiltonian V may 
be rewritten in terms of pseudopotentials as 

oo 

V = ]T>(N - z i\) =EE v m ,V m '(ij), (4.10) 

i<j i<j m'—O 

where the operator V m ' {ij) projects the electron pair ij onto the relative angular 
momentum state m' . Notice that due to the factor Yik<i ( z k — z i)™ m Laughlin's 
wave function (|4.3[) . no particle pair is in a relative angular momentum state 

9 This is similar to the average value of the radius at which the electron's guiding centre is 
placed in the symmetric gauge (see Sec. 12.4.11 . Remember (e.g. from classical mechanics) that 
the decomposition of a two-particle wave function in relative and centre-of-mass coordinates 
maps the two-body problem to an effective one-body problem. 
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m' < m. If one then chooses, though somewhat artificially, all pseudopotentials 
with a to < to' to be positive (say 1) and all others zero, 

m I lor to > to K ' 

one obtains Vip^ = 0, i.e. Laughlin's wave function is the zero-energy eigenstate 
of the model (|4.11[) . Since the model describes an entirely repulsive interaction, 
all possible states must have an energy E > 0. Therefore, Laughlin's wave 
function is even the exact ground state of the model (|4.1ip . Furthermore, it 
is the only zero-energy state because if one keeps the total number of particles 
and flux fixed, any other state different from that described by Laughlin's wave 
function involves a particle pair in a state with an angular momentum quantum 
number different from to. If it is smaller than to, this particle pair is affected 
by the associated non-zero pseudopotcntial to' and thus costs an energy on the 
order of v m i > 0. If the particle pair is in a momentum state with to' > m. 
there is at least another pair with to" < to in order to keep the filling factor 
fixed, and this pair raises the energy. These general arguments show that any 
excited state involves a finite (positive) energy given by a pseudopotential v m > , 
with to' < to, which plays the role of an energy gap. In this sense, the liquid 
state described by Laughlin's wave function is indeed an incompressible state 
that already hints at the possibility of a quantum Hall effect if we can identify 
the correct quasi-particle of this TV-particle state that becomes localised by the 
sample impurities. 

Notice that the above considerations are based on an extremely artificial 
model interaction (|4. 1 1 [1 that has, at first sight, very little to do with the physical 
Coulomb repulsion. However, the model is often used to generate numerically 
(in exact-diagonali-sation calculations) the Laughlin state, which may then be 
compared to the Coulomb potential decomposed in Haldane's pseudopotentials. 
This procedure has shown that thcL aughlin state generated in this manner 
has an overlap of more than 99% with the state obtained from the Coulomb 
potential 63, 64], which is amazingly high for a wave function obtained from 
a well-educated guess. This high accuracy of Laughlin's wave function may 
be understood in the following manner: when one decomposes the Coulomb 
interaction potential in Haldane's pseudopotentials, one obtains a monotonically 
decreasing function when plotted as a function of to (see Fig. 14. 2p . Furthermore, 
the component v\ is much larger than v$ and all other pseudopotentials v m 
with higher values of to These higher terms may be treated in a perturbative 
manner and do not change the ground state which is protected by the above- 
mentioned gap on the order of v\ > w m , with to > 1. 

Furthermore, we mention that, apart from its successful verification by exact- 
diagonalisation calculations [63l [64] , Laughlin, in his original paper [14] , showed 
within a variational calculation that the quantum liquid described by his wave 



10 One has vi/v$ ~ 1.6 in the LLL. Notice that pseudopotentials with even angular momen- 
tum quantum number m do not play any physical role because of the fermionic nature of the 
electrons. 
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function (|4.3[) has indeed a lower energy than the previously proposed WC. 
Again the reason for this unexpected feature is the capacity of Laughlin's wave 
function, which varies as r 2m when two particles i and j approach each other 
with r = \zi — Zj\, to screen the short-range components of the interaction 
potential. Notice that for a WC of fermions, the corresponding iV-particle wave 
function decreases as r 2 , as dictated by the Pauli principle. 

4.2.3 Quasi-particles and quasi-holes with fractional charge 

Until now, we have discussed some ground-state properties of Laughlin's wave 
function. We have seen that the Laughlin state at v = 1/to is insensitive to 
the short-range components of the interaction potential described by Haldane's 
pseudopotentials v m i with m! < to, whereas excited states must be separated 
from the ground state by a gap characterised by these short-range pseudopoten- 
tials. However, we have not characterised so far the nature of the excitations. 

There are two different sorts of excitations: (i) elementary excitations (quasi- 
particles or quasi-holes) that one obtains by adding or removing charge from the 
system, and (ii) collective excitations at fixed charge. The latter are simply a 
charge-density-wave excitation which consist of a superposition of particle-hole 
excitations at a fixed wave vector q (the momentum of the pair) and which may 
be shown to be gapped at all values of q. Its dispersion reveals a minimum 
(called magneto-roton minimum) at a non-zero value of the wave vector that 
indicates a certain tendency to form a ground state with modulated charge 
density, such as a WC. The characteristic dispersion relation of these collective 
excitations is shown in Fig. I4.3f a). However, we do not discuss collective 
excitations here and refer the interested reader to the literature for a more 
detailed discussion [55j E EL an( i concentrate here on a presentation of the 
elementary excitations. 

Quasi-holes 

Elementary excitations are obtained when sweeping the filling factor slightly 
away from v = 1/to. Remember that there are two possibilities for varying the 
filling factor: adding charge to the system by changing the electronic density or 
adding (or removing) flux by varying the magnetic field. Remember further [see 
Eq. ()4.4|) ] that the number of flux is intimitely related to the number of zeros 
in Laughlin's wave function. We therefore consider the ansatz 

N 

% h (Z , {zj, Z*}) = Y[( Zj - Zo) (fa, ( 4A2 ) 
3=1 

for an excited state. Each electron at the positions Zj thus "sees" an additional 
zero at zq. In order to verify that this wave function adds indeed another flux 
quantum to the system, we may expand Laughlin's wave function (|4.3p formally 
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Figure 4.3: (a) Dispersion relation for collective charge-density-wave excita- 
tions (Girvin, MacDonald and Platzman, 1986; Girvin, 1999). The continuous 
lines have been obtained in the so-called single-mode approximation (Girvin, 
MacDonald and Platzman. 1986) for the Laughlin states at v = 1/3, 1/5 and 
1/7, whereas the points are exact-diagonalisation results (Haldane and Rezayi, 
1985; Fano, Ortolani and Colombo, 1986). The arrows indicated the character- 
istic wave vector of the WC state at the corresponding densities, (b) Quasi- hole 
excitation. Each electron jumps from the state m to the next-higher angular 
momentum state m + 1 . 



in a polynomial, 

t l'm({ z jj z j}) = Oim 1 ,...,in N Z 1 1 ...Z N N e ^ ' ' I , 

{mi} 

where the a mi ,...,m N describe the expansion coefficients. We now choose the 
position zo at the centre of the disc, in which case the wave function of the 
excited state (|4.12[) simply reads 

*Pqh({zj , Zj }) = ^ ] Oi mi ,,,, tmN Z 1 ... Z N N+ e 3 ' 3 ' I , 

{mi} 

i.e. each exponent is increased by one, rrn — ► rtii + 1. This may be illustrated 
in the following manner: each clcctronjumps from the angular momentum state 
m to a state in which the angular momentum is increased by one (see Fig. 14. 3[) . 
leaving behind an empty state at m = 0. The excitation is therefore called a 
quasi-hole as we have already suggested by the subscript in Eq. (|4.12p . This 
also affects the quantum state with highest angular momentum M, i.e. we 
have increased the sample size by the surface occupied by one flux quantum, 
while keeping the number of electrons fixed Furthermore, this quasi-hole is 



11 Naturally, the total surface of the quantum Hall system remains constant, but physically 
we have slightly increased the B-field. Each quantum state occupies then an infinitcsimally 
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associated with a vorticity if one considers the phase of the additional factor in 
Eq. (|4~T2l) . 

N 

iP qh (z = 0, {z 3 , z*}) <x J] e""' x ^ ({*,-, «;}) , 

3 

i.e. each particle that circles around the origin z = experiences an additional 
phase shift of 2ir as compared to the original situation described in terms of 
Laughlin's wave function (|4.3|) . This is reminiscent of the vortex excitation in 
a type-II superconductor [66] . 

We have seen above that one can create a quasi-hole excitation at the postion 
Zq by introducing one additional flux quantum, Nb — > Nb + 1 , which lowers the 
filling factor by a tiny amount. However, we have not yet determined the charge 
associated with this elementary excitation. This charge may be calculated by 
considering the filling factor fixed, i.e. we need to add some (negative) charge 
to compensate the extra flux quantum in the system. From Eq. (|4.4p . we notice 
that the relation between the extra flux ANb and the compensating extra charge 
AN is simply given by 

AiV R 

mAN = AN B & AN = . (4.13) 

m 

This very important result is somewhat unexpected: in order to compensate one 
additional flux quantum (ANb = 1), one would need to add the m-th fraction of 
an electron. The charge deficit caused by the quasi-hole excitation is therefore 

e* = -, (4.14) 
m 

i.e. the quasi-hole carries fractional charge. 
Quasi-particles 

In the preceding paragraph, we have considered a quasi-hole excitation that is 
obtained by introducing an additional flux quantum in the system [or, math- 
ematically, an additional zero in the Laughlin wave function, see Eq. (|4.12|) ]. 
Naturally, one may also lower the number of flux quanta by one in which case 
one obtains a quasi-particle excitation with opposite vorticity as compared to 
that of the quasi-hole excitation. This opposite vorticity suggests that we use a 
prefactor n^=i( 2 j' — z o)> instead of n^Li( z i — z o) as m the expression (|4.12p . 
in order to create a quasi-particle excitation at the position z . Remember, 
however, that the resulting wave function would have unwanted components in 
higher LLs because the analyticity condition of the LLL is no longer satisfied. 
In order to heal the quasi-particle expression, one formally projects it into the 
LLL, 

N 

i> qp (z , {z 3 , z*}) = V LLL Y[(z* - z*) ^ ({*,-,*;}) • (4.15) 

3=1 

smaller surface 2tt1 2 b , such that the system may accomodate for one more quantum state, 
M = N B -> N B + 1. 
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Figure 4.4: Experimental setup for the observation of fractionally charged 
quasi-particles. In addition to the usual geometry, one adds, at the upper and 
the lower edges, side gates that are used to deplete the region around the gates 
by the application of a voltage V sg . The filling factor is chosen to be v = 1/3. 
As a result, the edge states at the opposite edges are brought into close vicinity. 
(a) Weak-backscattering limit. The incompressible liquid has a bottleneck at 
the side gates, i.e. the edges are so close to each other that a tunneling event 
between them has a finite probability. A particle injected at the left contact 
may thus be backscattered (grey arrow) in a region filled by the incompressible 
Laughlin liquid, although the majority of the particles reaches the right contact 
(black arrows), (b) Strong-backscattering limit. If one increases the side-gate 
voltage V sg , the incompressible v = 1/3 liquid is eventually cut into two parts 
separated by a fully depleted region (v = 0). In this case, backscattering is the 
majority process (black arrow), and a tunneling may occur over the depleted 
region such that a particle injected at the left contact may still reach the right 
one (grey arrows). 

There are several manners of taking into account this projection Vlll- A com- 
mon one consists of replacing each occurence of the non-analytic variables z* 
(and powers of them) in the polynomial part of the wave function by a deriva- 
tive with respect to Zj in the same polynomial [67] . By partial integration, this 
amount to deriving the Gaussian factor by (d Zj ) m which, up to a numerical 
prefactor, yields exactly the non-analytic polynomial factor z* m . We will en- 
counter this projection scheme again in the discussion of the CF generalisation 
of Laughlin's wave function (Sec. 14.4. 1|) . 

4.2.4 Experimental observation of fractionally charged quasi- 
particles 

That the fractional charge of Laughlin quasi-particle^jfl is not only a mathe- 
matical concept but a physical reality has been proven in a spectacular manner 

12 From now on, we use the term "(Laughlin) quasi-particles" generically in order to denote 
quasi-particles and quasi-holcs. 



78 



Strong Correlations and the Fractional Quantum Hall Effect 



in so-called shot-noise experiments on the v = 1/3 FQHE state [55J IM] El In 
these experiments, one constrains the quantum Hall system with the help of 
side gates (see Fig. I4.4[) that are used to deplete the region in their vicinity via 
the application of a gate voltage V sg . As a consequence of this depletion the 
quantum Hall system has a bottleneck where the corresponding edge states are 
brought into spatial vicinity [Fig. 14. 4f a)] or where the incompressible quantum 
liquid may even be cut into two parts separated by a completely depicted bar- 
rier [Fig. I4.4f b)] . In the first case, an injected charge may be backscattered in 
a tunneling event at the bottleneck over a region filled by the v = 1/3 liquid 
(weak-backscattcring limit). If one increases the side-gate voltage, the incom- 
pressible liquid becomes eventually cut into two parts separated by a completely 
depleted barrier, and one obtains the strong-backscattering limit. 

In a shot-noise measurement, one does not only measure the average current 
I (over a certain time interval) but simultaneously the (square of the) current 
fluctuation (A/) 2 which is proportional to the carrier charge. If the elementary 
charged excitations arc e* = e/3 quasi-particles and not electrons, one may ex- 
pect to measure this particular charge. The experiments j68j [69] have indeed 
shown that the charge measured in the shot noise is e* — e/3 if the tunneling 
process takes place at a bottleneck filled with the incompressible quantum liq- 
uid [Fig. 14.4( a)]. whereas it is the usual elementary charge e in the case of a 
tunneling process over a depleted region [Fig. 14.4( b)] . 

4.2.5 Laughlin's plasma analogy 

A compelling physical picture of Laughlin's wave function (|4.3p and the prop- 
erties of its elementary excitations (|4. 12[) and (|4.15[) with fractional charge has 
been provided by Laughlin himself |14j , in terms of an analogy with a classical 
2D one- component plasma. In the present subsection, we present the basic ideas 
and results of this plasma analogy, for completeness and pedagogical reasons. 
However, no new results will come out of this analogy here, as compared to 
those derived above. 

Remember from basic quantum mechanics that the modulus square of a 
quantum-mechanical wave function may be interpreted as a statistical probabil- 
ity distribution. For Laughlin's wave function ()4.3[) . one obtains the probability 
distribution 

\^({^})\ 2 = U^-^ 2m ^ M2/2 - 

Now, remember from classical statistical mechanics that a probability distribu- 
tion in the canonical ensemble is the Boltzmann weight, exp(— (3Tt), of some 
Hamiltonian TL and that the classical partition function, which encodes all rele- 
vant statistical information, is obtained from a sum over the Boltzmann weights 
of all possible configurations C, Z = ^ c exp[—(3H(C)]. Laughlin's plasma anal- 
ogy consists precisely of the identification of the modulus square of his wave 



Later this kind of experiment has been repeated for other FQHE states. 
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function with the Boltzmann weight of some mock Hamiltonian L/^o The 
mock Hamiltonian may be obtained exactly from this identification, 

-/3U cl = ]n\^({z j })\ 2 , (4.16) 
and one obtains, by choosing somewhat artificially (3 = 2/qF^l 

U cl = -a 2 J2 In \ Zi - Zj \ +q J2 l -^f- (4.17) 

i<3 3 

This is nothing other than the classical Hamiltonian of a 2D one-component 
plasma, in terms of the charge 

q = m = 2s + l (4.18) 

of the plasma particles. The first term of Eq. (|4.17[) reflects the interactions 
between the charged plasma particles, whereas the second term describes their 
interaction with a neutralising background of positive charge, as in the case of 
the jellium model of the Coulomb gas [SH [55]. This may be seen best with 
the help of Poisson's equation, — Acf) = 2irqn q (r), for an electrostatic 2D po- 
tential due to the charge density qn q . The first term describes then indeed 
particles with charge q interacting via the 2D Coulomb interaction potential 
4>(r) = — ln(|r|//s), and the second term is the interactions with the neutralis- 
ing background because A|r| 2 /4/^ = 1/1% = 2ttub, where the flux density ns 
may thus be viewed as the charge density of the positively charged background. 

In order to minimise the energy of the mock Hamiltonian U c i, which cor- 
responds to a distribution of highest weight, the 2D plasma thus needs to be 
charge-neutral, i.e. the charge density of the plasma particles qn e i must be 
compensated by that of the background n b , 

n B -qn e i=Q, (4.19) 

which, together with Eq. (|4.18[) , yields nothing other than the relation between 
the filling factor v and the exponent in Laughlin's wave function (14. 5p . v = 
n e i/riB = l/m. 

The plasma analogy does not only apply to the ground-state wave func- 
tion (|4.3p but also to the quasi-hole excitation (|4.12|) . The additional factor 
rij=i( z j — z o) m the quasi-hole wave function (|4.12j) yields, within the plasma 
analogy (|4.16p . an additional term 

JV 

V = -q^tofa - zo\ (4.20) 

j=i 

to the mock Hamiltonian (|4. 17|) . U c i — > U c i + V. This additional term may 
be interpreted as the interaction of the plasma particles with an "impurity" 



1 mock: singlish for fake; mainly used the description of Singaporean catering food. 
15 Notice that Laughlin's wave function describes a system at T = 0, such that temperature 
does not intervene in the expressions. The choice is purely formal. 
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Figure 4.5: (a) Process in which a particle A moves on a path C around another 
particle B. In three space dimensions, one may profit from the third direction 
(z-direction) to lift the path over particle B and thus to shrink the path into 
a single point, (b) Process equivalent to moving A on a closed path around B 
which consists, apart from a topologically irrelevant translation, of two succes- 
sive exchanges of A and B. 



of unit charge at the position zq. In order to maintain charge neutrality, the 
impurity needs to be screened by the plasma particles. Since the charge of each 
plasma particle is q = m = 2s + 1 and thus greater than unity, one needs 1/q 
plasma particles to screen the impurity of charge one. Remember that each 
plasma particle represents one electron of unit charge in the original Laughlin 
liquid. One therefore obtains the same charge fractionalisation of the Laughlin 
quasi-particle (|4.14p . e* = e/m, as in the original quantum model. 

4.3 Fractional Statistics 

4.3.1 Bosons, fermions and anyons — an introduction 

One of the most exotic consequences of charge fractionalisation in 2D quantum 
mechanics, exemplified by Laughlin quasi-particlcs, is fractional statistics. Re- 
member that, in three space dimensions, the quantum-mechanical treatment of 
two and more particles yields a supcrsclection rule according to which quantum 
particles are, from a statistical point of view, either bosons or fermions. This 
supcrsclection rule is no longer valid in 2D (two space dimensions), and one 
may find intermediate statistics between bosons and fermions. The correspond- 
ing particles are called anyons, because the statistics may be any. The present 
section is meant to illustrate these amazing aspects of 2D quantum mechanics, 
and we try to avoid a too formal or mathematical treatment. We refer, again, 
the interested reader to the more detailed literature [70] . 

In order to illustrate the different statistical (i.e. exchange) properties of 
two quantum particles in three and two space dimensions, consider a particle A 
that moves adiabatically on a closed path C in the xy-planc around another one 
B of the same species (see Fig. 14. 5|) . We choose the path to be sufficiently far 
away from particle B and the two particles to be sufficiently localised such that 
we can neglect corrections due to the overlap between the two corresponding 
wave functions. Notice first that such a process T is equivalent, apart from a 
topologically unimportant translation, to two successive exchange processes £ , 
in which one exchanges the positions of A and B. Algebraically, this may be 
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expressed in terms of the corresponding operators as 



£ 2 =T 



or 



£ = ±VT, 



(4.21) 



modulo a translation. 

Let us discuss first the three-dimensional case. Because of the presence of the 
third direction (z-dircction), one may elevate the closed path in this direction 
while keeping the position of particle A fixed in the xy plane. We call the 
elevated path C . Furthermore, one may now shrink the closed loop C into a 
single point at the position A without passing by the position of particle B which 
remains in the xy-plane. This final (point-like) path is called C" . Although this 
procedure may seem somewhat formal, a quantum-mechanical exchange process 
does principally not specify the exchange path in order to define whether a 
particle is a boson or a fermion, but only its topological properties. From a 
topological point of view, all paths that may be continuously deformed into each 
other define a homotopy class [71] . Equation (|4.21[) must therefore be viewed as 
an equation for homotopy classes in which a simple translation and an allowed 
deformation are irrelavant. As a consequence of these considerations, the simple 
point-like path C" at the position of particle A, which may be formally described 
by C" = 1, is in the same homotopy class as the original path C. Therefore, the 
associated processes are the same, and one has 



where the last equation is symbolic in terms of the one operator. It indicates 
that the quantum- mechanical operator £, corresponding to particle exchange, 
has two eigenvalues that are the two square roots of unity, es = exp(2i7r) = 1 
and &f = exp(z7r) = — 1. This is precisely the above-mentioned superselection 
rule, according to which all quantum particles in three space dimensions are 
either bosons (eg = 1) or fcrmions [ep = —1). 

In two space dimensions, this topological argument yields a completely dif- 
ferent result. It is not possible to shrink a path C enclosing the second particle B 
into a single point at the position of A, without passing by B itself. This means 
that the position of B must be an element of the path at a certain moment of 
the shrinking process, which cannot profit from a third dimension in order to 
elevate the loop on which it moves above the xj/-planc. The single point still 
represents a homotopy class of paths, but these paths do not enclose another 
particle, and C is therefore an element of another homotopy class, i.e. the one 
of all paths starting from A and enclosing only the particle B. If there are more 
than two particles present, the homotopy classes are described by the integer 
number of particles enclosed by the paths in this class. From an algebraic point 
of view, the exchange processes are no longer described by the two roots of unity 
1 and —1, but by the so-called braiding group, and the classification in bosons 
and fermions is no longer valid. In the simplest case of Abelian statisticsOf] one 

16 There are more complicated cases of non- Abelian statistics, in which the exchange pro- 
cesses of more than two different particles no longer commute, but we do not discuss this case 
here and refer the reader to the review by Nayak et aJ. 170) . 



T = T(C) = T(C") = 1 and thus £ = VI, 



(4.22) 
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needs to generalise the commutation relation 



^(rOV'fo) = ±V>(r 2 )^(ri), 



(4.23) 



for bosons and fermions, respectively, to 



V>(ri)?/>(r 2 ) 



(4.24) 



where a is also called the statistical angle. One has a = for bosons and 
a = 1 for fermions, and all other values of a in the interval between and 2 for 
anyons. Sometimes anyonic statistics is also called fractional statistics - indeed 
all physical quasi-particles, such as those relevant for the FQHE, have an angle 
that is a fractional (or rational) number, but there is no fundamental objection 
that irrational values of the statistical angle should be excluded. 

Before discussing the anyonic nature of Laughlin quasi-particles, we need to 
mention an important issue in these statistical considerations. We know that 
fermions are forced to satisfy Pauli's principle which excludes double occupancy 
of a single quantum state, whereas the number of bosons per quantum state is 
unrestricted. What about anyons then? In the context of quantum fields the 
Pauli principle yields, via Eq. (|4.23|) for r = ri = r 2 , 

ip(r)ip{r) = 0. 

For an arbitrary statistical angle, one obtains in the same manner, from Eq. 



which may be viewed as a generalised Pauli principle for 2D anyons |72j . Only 
if a = modulo 2, we may have ip(r)ip(r) ^ in order to satisfy Eq. (|4.25|) . 
Otherwise, when a / modulo 2, we necessarily have ip(r)ip(r) = 0. Anyons 
are, thus, from an exclusion-principle point of view more similar to fermions 
than to bosons. 

4.3.2 Statistical properties of Laughlin quasi-particles 

We may now apply the above general statistical considerations to the case of 
Laughlin quasi-particles. The basic idea is to describe the statistical angle as 
an Aharonov-Bohm phase due to some gauge field that is generated by the flux 
bound to the charges included in a closed loop <9E. This closed loop, around 
which a quasi-particle moves adiabatically, encloses a surface E. The gauge field 
is not to be confunded with the one which generates the true magnetic field B 
- it is rather a mock (or fake) field Am (with B^/ = V x Am) that generates 
the flux bound, e.g., by the electrons in the Laughlin liquid via the relation 
(|4.4[) . We consider the case where the area E is filled with A e ;(E) electrons 
condensed in an incompressible quantum liquid described by Laughlin's wave 
function (|4. 3|) and A^^(E) quasi-hole excitations (|4.12p . such that there are two 
contributions to Bm — |Bm|> 




(4.25) 



B M E = AW- = [mA e; (E) + N qh (T,)] - 



(4.26) 
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The corresponding Aharonov-Bohm phase, which the quasi-particle picks up 
when turning around the area E on the boundary path <9E, is given by 

r A _ B = 27r^-/ dr ■ A M (r) = 27r^- / d 2 rB M (r), 

where e* = e/m is the charge of the quasi-particle and where we have used 
Stoke's theorem to convert the line integral of Am on <9E into a surface integral 
of Bm over the area E. The Aharonov-Bohm phase has therefore two contri- 
butions, one T e i that stems from the electrons condensed in the Laughlin liquid 
and the other one T qh that is due to the enclosed quasi- holes. One obtains from 
Eq. fOffl l 

T el = 2ir—mN el = 27rN el , (4.27) 
e 

for the enclosed electrons, i.e. an integer times 2ir. Notice that this contribution 
to the Aharonov-Bohm phase may not be interpreted in terms of a statistical 
angle because it does not describe a true exchange process: the involved particles 
are not of the same type - we have chosen a quasi-particle to move on a path 
enclosing condensed electrons. However, had we chosen an electron rather than 
a quasi-hole to move along the path 9E, the Aharonov-Bohm phase, 

e 

r e j_ ei = 27T-TOA e/ (E), 

e 

would give rise to a statistical angle a = mN„r (E)Fl If we have only one electron 
enclosed by the path, 7V e /(E) = 1, the statistical angle is simply the odd integer 
m, which is equal to 1 (modulo 2), as it should be for fermions. 

A more interesting situation arises when the path encloses Laughlin quasi- 
holcs, in which case the Aharonov-Bohm phase reads 

T el = 2n-N qh = 2n^. (4.28) 

e m 

Consider a single quasi-hole in the area E, N q h = 1: one encounters the rather 
unusual situation in which the Aharonov-Bohm phase is a fraction of 2-7T, and 
the associated statistical angle is a = \ jm. This illustrates that Laughlin quasi- 
holcs arc indeed anyons with fractional statistics, as we have argued above. 

4.4 Generalisations of Laughlin's Wave Function 

Although Laughlin's wave function (|4.3| has been extremely successful in the de- 
scription of the FQHE at v = 1/3 and 1/5, it is not capable of describing all ob- 
served FQHE states. Indeed, there are e.g. FQHE states at v — 2/5, 3/7, 4/9, ... 
corresponding to the series p/(2p + 1), or more generally to p/(2sp + 1), in 

17 Remember that the statistical angle is defined with respect to an exchange process £ 
which is the square root of the process T considered here [Eq. I|4,21j l] . The relation between 
the statistical angle and the Aharononv-Bohm phase is therefore T = 2ira and not net. 
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Figure 4.6: Schematic view of composite fermions. The electronic state at 
v = 1/3 may be interpreted as a CF state at an integer CF filling factor v* = 1, 
where each vortex bound to an electron carries 2s (here s = 1) flux quanta. In 
the same manner a CF filling factor v* = 2 gives rise to an (electronic) FQHE 
state at v = 2/5. 



terms of the integers s and p, which may be accounted for within composite- 
fermion (CF) theory, which we present below. Furthermore, even-denominator 
FQHE states have been observed at v = 5/2 and 7/2 [17], in the first excited 
LL (n = 1), and, in wide quantum wells or bilayer quantum Hall systems, at 
v = 1/2 and v = 1/4 [73l[74j. Whereas the latter may be understood within a 
multi-component picture, which we will briefly introduce in Chap. [5j the states 
at v = 5/2 and 7/2 may find their explanation in terms of a so-called PfafEan 
wave function. Both the CF and the Pfaffian wave functions are sophisticated 
generalisations of Laughlin's original idea. 



4.4.1 Composite Fermions 

Soon after the discovery of the most prominent FQHE state at v = 1/3, a 
lot of other states have been observed at the filling factors v — p/(2sp + 1). 
In a first theoretical approach, these states were interpreted in the framework 
of a hierarchy scheme [62l [75] according to which the quasi-particles of the 
Laughlin (parent) state, such as v — 1/3, condense themselves into a Laughlin- 
type (daughter) state, due to their residual Coulomb repulsion - remember that 
the Laughlin quasi-particles are charged with charge e* = e/m. In this picture, 
the 2/5 state would be the daughter state formed of Laughlin quasi-particle 
excitations of the 1/3 state. 

An alternative picture, though related to the above-mentioned hierarchy 
scheme, was proposed by Jain in 1989 [TSJ, [15]. The basic idea consists of a 
reinterpretation of Laughlin's wave function (|4.3[) : consider only the polynomial 
part, the Gaussian cxp(— |zj| 2 /4) being an ubiquitous factor which finally 
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needs to be multiplied with the polynomial wave function, 

k<l k<l k<l 

In the last step of this equation, we have split the product into two parts, one 
with the exponent 2s, which we call the vortex part, and another one with the 
exponent 1. 

Before introducing Jain's generalisation, let us interpret the above wave 
function in terms of the statistical properties introduced in the last section. 
Quite generally, one may express any LLL TV-particle wave function i^lll as a 
product of such a vortex factor and another (residual) wave function ip res , 

TpLLL ({Zj}) = ]"[(<2fc ~ Zl) m ^res {{Zj}) ■ 
k<l 

If the original wave function is fermionic, i.e. anti-symmetric with respect to an 
exchange process of an arbitrary particle pair, the symmetry properties of ipres 
depend on the parity of m'. If it is odd, i/j rea must be a symmetric (bosonic) 
wave function, and if it even, both the original and the residual wave functions 
are anti-symmetric (fermionic). In terms of the above-mentioned gauge field 
Aju(r), the statistical angle associated with the vortex factor is just given by 
the parity of m', which may be viewed as the number of flux quanta attached 
to each particle at the positions Zj. Flux attachment may thus be used, in 2D 
quantum mechanics, to transform fermions into bosons and vice versa. 

In the case of the above decomposition (|4.29|) of Laughlin's wave function, 
the vortex part attaches s pairs of flux quanta to each particle position and 
therefore does not affect the statistical properties of the wave function. The 
second factor 

X>=i({zj}) = Y[(z k - z{) 

k<l 

is indeed fermionic and corresponds, as we have mentioned in Sec. 14.2. l\ to a 
completely filled LL at a virtual (CF) filling factor of v* = 1, the true filling 
factor being still v = l/(2s + I). This is schematically represented in Fig. 14.61 
Jain's generalisation consists of replacing the term Ilfc<i( Zfe — Zl ) by any 
other Slater determinant Xw=p{{ z ji z j}) °f P completely filled LLs, with a CF 
filling factor v* = p, 

t J ({ Zj , z*}) = V LLL J] (z k - Zl ) 2s xw=p{{zj, z*}), (4.30) 

k<l 

where we need to take into account the same projection Vlll to the LLL as in 
the case of quasi-particlc excitations (|4.I5|) because, contrary to the v* = 1 case, 
the wave function Xv*=p{{zj, z*}) has by construction non-analytic components, 
i.e. components in higher (CF) LLs. 

Jain's wave function (I4.30jl may be illustrated in the following manner. Via 
the first factor nfc<z( z fc — z () 2s i we have effectively bound 2s flux quanta to 



8(> 



Strong Correlations and the Fractional Quantum Hall Effect 



each of the electrons, as we have already mentioned above. This novel type of 
particle is what we call the composite fermion (CF). The residual (free) flux 
quanta effectively determine the effective number of states per (CF) LL, 

N B ^N* B =N B -2sN el , 

which correspond to a renormalised magnetic field 

B^B* = B-2s(^)n e i. (4.31) 

Similarly the CF filling factor is defined with respect to the renormalised number 
of flux quanta, 

*/■ = — v*- 1 = v- 1 -2s, (4.32) 



which leads to the relation 



v 



v = - (4.33) 

2sv* + 1 v ; 

between the CF filling factor and the usual one v [Eq. (|2.44p ]. For completely 
filled LLs, v* = p, this yields the above-mentioned series 

v = — ^ — (4.34) 
2sp +1 y ' 

for the FQHE states which may thus be interpreted as IQHE states of CFs. To 
be explicit, the physical picture of CF theory is the following: the ground state 
is described by the wave function (|4.30p . which describes an incompressible 
quantum liquid in the same manner as Laughlin's wave function does. The 
elementary excitation in the CF theory consists of a CF promoted to the next 
higher CF LL, which is separated from the ground state by an energy gap, in 
analogy with the electron as compared to n completely filled (electronic) LLs 
in the IQHE A gain, these elementary CF excitations become localised by the 
sample impurities, and one therefore obtains a plateau in the Hall resistance 
which is thus quantised. 

Numerically, Jain's CF wave function (|4.30[) has been successful in the de- 
scription of the series (|4.34p of FQHE states: even if the overlap with the exact 
ground states decreases when the quantum number p, which describes the num- 
ber of completely filled CF LLs, increases, the overlap is still reasonably high 
(above 95%) for the number of particles accessible in state-of-the-art exact di- 
agonalisation calculations. Notice, however, that the physical interpretation is 
more involved as compared to Laughlin's wave function, because of the LLL 
projection Vlll, which is rather complicated to implement in analytical as well 
as numerical calculations. For a further review of CF theory, we refer the inter- 
ested reader to the literature. The above-mentioned wave-function approaches 
are thoroughly reviewed in Jain's recent book [76] . Furthermore, there have 



18 Remember, however, that the energy scale of this gap is not given in terms of a kinetic 
energy heB/m, but in terms of the Coulomb interaction e 2 /elg. 
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been field-theoretical approaches beyond the numerical wave-function descrip- 
tion presented above, such as in terms of Chern-Simons theories [771 178] or in 
terms of a Hamiltonian theory [5] . For a review of these complementary the- 
ories we refer the reader to the book edited by Heinonen [79] or the excellent 
pedagogical review by Murthy and Shankar [5] . 

4.4.2 Half-filled LLs and Pfaffian states 

Within the CF picture, we have seen that the effective magnetic field becomes 
rcnormalised due to flux attachment [Eq. (|4.31[) ]. An interesting situation arises 
when the filling factor is v = 1/2, which corresponds to the limit p — > oo in Eq. 
(|4.34p . In this limit the effective magnetic field (|4.31|) vanishes, B* = 0, and 
one may then expect the corresponding phase to be described in terms of a 
metallic state, such as a Fermi liquid that one would obtain for electrons when 
the magnetic field vanishes. A natural ansatz for the TV-particle wave function 
of such a Fermi-liquid state is given by the Slater determinant 

^ = det (e lk ^), 

where the N electrons occupy the states described by the wave vectors k, 
i = 1,...,N, the modulus of which is delimited by the Fermi wave vector 
|k;| < kp, and is the position of the j-th particle. Notice that this state 
is nevertheless unappropriate in the description of a state in the LLL. Indeed, 
if the scalar product in the exponent is rewritten in terms of complex variables, 
ki • Yj = (k[Z* + k*Zj)/2, one realises that the Fermi-liquid state violates the 
LLL condition of analyticity. Formally, one may again avoid this problem by 
projecting the Fermi-liquid state into the LLL, and one obtains indeed a state, 

Vvl 1/2 = Tlll I] (Zk ~ zif det (e ik * r >) , (4.35) 
k<i 

that was proposed by Rezayi and Read for the description of a compressible 
metallic state at v = 1/2 [80]. The first term is the same factor as in CF theory, 
which attaches 2 flux quanta to each particle and which cancels thus the external 
magnetic field, B* = B - 2(h/e)n e i = 0. 

Because the wave function (|4.35|) describes a compressible state, one should 
not observe a quantised Hall resistance, in agreement with most experimental 
data. A FQHE at v = 1/2 (and 1/4) has only been observed in very wide 
quantum wells [731 174] , which are likely to be described by two-component wave 
functions [HI] that we will briefly introduce in Chap. 

In contrast to the LLL, the half-filled LL n = 1 reveals, in both spin branches, 
a FQHE (5/2 and 7/2 states). The difference between the half-filled LL n = 
and n = 1 is due to a different effective interaction potential that takes into 
account the wave function overlap between two (interacting) particles, which 
we do not discuss in detail here. Indeed, the Fermi-liquid-like state ()4.35[) turns 
out to be quite unstable with respect to particle pairing. This is reminiscent of 
the BCS (Bardeen-Cooper-Schrieffer) instability of a conventional Fermi liquid 



8S 



Strong Correlations and the Fractional Quantum Hall Effect 



that gives rise to superconductivity [5H [55], although the glue between the 
particles is no longer a phonon-mediated attractive interaction, but only the 
repulsive Coulomb interaction in a strong magnetic field. As we have already 
mentioned in Sec. I4.2.2( such an interaction may yield a discrete two-particle 
spectrum, in contrast to a repulsive interaction in the absence of a magnetic 
field. As a consequence, pairing may occur at certain relative angular momenta 
for particular pseudopotential sequences and for sufficiently high filling factors F^l 
In the present case, one may exclude s-wave pairing, i.e. in the relative angular 
momentum state with m = due to the Pauli principle, and the most natural 
candidate would therefore be p-wave pairing in the relative angular momentum 
state m = 1 [19] . 

A wave function that accounts for p-wave pairing was proposed by Moore 
and Read in 1991 [18], 



where we have again omitted the ubiquitous Gaussian factor. As for the CF 
wave functions (|4.30j) and the Rczayi-Read wave function (|4.35j) . the factor 
n fc<; (zfc — zi) 2 attaches two flux quanta to each electron and therefore does 
not change the statistical properties of the wave function. If the wave function 
consisted only of this factor (times the Gaussian), one would have a bosonic 
Laughlin wave function that describes an incompressible quantum liquid at the 
desired filling factor v = 1/2. However, it does not have the correct statistical 
properties. This problem is healed by the first factor Pf[l/(zj — Zj)] which 
represents the Pfaffian of the N x N matrix M.^ — l/(zi — Zj). The Pfaffian 
may be viewed as the square root of the more familiar determinant, Pf(.M) = 
^Jdet(M.), and has the same anti-symmetric properties as the determinant in 
an exchange of two particles i and j, such that it generates a fcrmionic wave 
function. Notice, furthermore, that this Pfaffian seems, at first sight, to take 
away some of the zeros such that one could expect the filling factor to increase. 
However, the function Ilk^C 2 * — z i) 2 1S a product of N(N — 1) ~ TV 2 terms, 
whereas the Pfaffian is a sum of products of N/2 ~ N terms. Therefore, the 
number of zeros, and thus the filling factor, is unchanged in the thermodynamic 
limit, ff-tco. 

A particularly interesting feature of the Pfaffian state are the quasi-particle 
excitations of charge e/4 which satisfy non-Abelian anyonic statistics |18j . in 
contrast to the corresponding excitations of Laughlin's (|4.3p or Jain's (|4.30[) 
wave functions. These non-Abelian quasi-particles are currently investigated 
in detail within the proposal of topologically-protected quantum computation 
[83] . A more detailed discussion of this issue is beyond the scope of these lecture 
notes, and we refer the reader to the review article by Nayak et al. [TO] . 



There have been attempts in the literature to formalise this point [63 , 82 . 




(4.36) 



Chapter 5 



Brief Overview of 
Multicomponent 
Quantum-Hall Systems 

5.1 The Different Multi-Component Systems 
5.1.1 The role of the electronic spin 

In the preceding chapters, we have completely neglected the physical conse- 
quences of possible internal degrees of freedom, apart from an occasional degen- 
eracy factor that has been smuggled in to account for experimental data. This 
choice has been made simply for pedagogical reasons, but it is clear that one 
prominent internal degree of freedom - the electronic spin - may not be put 
under the carpet so easily. Naively, one may expect that each LL is split into 
two distinct spin-branches separated by the energy gap Az due to the Zeeman 
effect. If this gap is large, one may use the same one-particle arguments as in 
the case of the IQHE, but now for each spin branch separately: once the low- 
est spin branch of a paticular LL is completely filled, additional electrons must 
overcome an energy gap that is no longer given by the LL separation but by 
Az- This would indeed not change the presented explanation of the IQHE 
instead of a localised electron in the next higher LL, one simply needs to invoke 
localisation in the upper spin branch. 

Also in the case of the FQHE, the explanation would need to be modified 
only in the fine structure if the Zeeman gap is sufficiently large. If the electrons 
fill partially the lower spin branch of the lowest (or any) LL, one may omit all 
transitions to the upper spin branch and argue that they constitute the high- 
energy degrees of freedom, in the same manner as inter-LL excitations in the 
case of the "spinless" fermions which we have discussed in Sec. 14.11 

However, the situation is not so easy as the above picture might suggest. 
Indeed, already in 1983 Halpcrin pointed out [53] that the Zeeman gap in GaAs, 
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with a ^-factor of g = -0.4, is A z = gn B B = g(he/2m )B ~ 0.33B[T] K and 
therefore much smaller than both the LL separation hue = (he/m)B ~ 24B[T] 
K, due to the rather small band mass (m = 0.068too, in terms of the bare 
electron mass mo, in GaAs), and the Coulomb energy scale Vc = e 2 /els — 
50y/S[T] K with a dielectric constant of e ~ 13. For a characteristic field of 6 
T, for which one typically reaches the LLL condition v = 1, one therefore has 
the energy scales 

e 2 

A Z ~2K < — -120K < ^wc^140K, (5.1) 

in GaAs. The situation is qualitatively the same in graphene, where one finds 
for a fielcfl of 6 T 

A Z ~7K < -^-~ 620K < v^^^lOOOK, (5.2) 
eIb Ib 

for g ~ 2 and e ~ 2.5, which are the appropriate values for graphene on a Si02 
substrate @ 

The inevitable consequence of these considerations is that, even if one may 
neglect the kinetic energy scale in a low-energy description of a partially filled 
LL, one cannot do so with the Zeeman energy scale. One must therefore take 
into account the electron spin within a two-component picture in which each 
quantum state \n,m) is doubled, |n,m;er) with a =\ and |. 

5.1.2 Graphene as a four-component quantum Hall sys- 
tem 

Another multi-component system that we have already discussed is precisely 
graphene, not only because of the tiny Zeeman gap which requires to take into 
account the electronic spin, but also because of its double valley degeneracy due 
to the two inequivalcnt Dirac points situated at the corners K and K' in the 
first BZ. Each quantum state |n,m) therefore occurs in four copies, \n,m;a) 
with a = (K, t), {K, j), (K', f) and (K', |). Formally this four- fold degeneracy 
may be described with the help of an SU(4) spin, whereas the two-fold spin 
degeneracy in GaAs, e.g., is represented by the usual SU(2) spin. Notice that 
it is very difficult in graphene to lift the valley degeneracy, and the associated 
energy scale is expected to be on the same order of magnitude as the Zeeman 
gap, i.e. it is tiny with respect to the one set by the Coulomb interactions. 

5.1.3 Bilayer quantum Hall systems 

A third multi-component system that we would like to mention consists of a dou- 
ble quantum well [see Fig. I5.1f a)]. These bilayer systems, which arc fabricated 

1 Remember that this field is somewhat arbitrary because the situation v = 1 may also be 
obtained easily for other fields by varying the gate voltage Va- 

2 Naturally, the dielectric constant depends on the dielectric environment around the 
graphene sheet and thus also on the substrate. 
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(a) 



-w- 



(b) 



w- 




Figure 5.1: (a) Profile of a double quantum well. The two wells are separated 
by a distance d that is typically on the same order of magnitude as the well 
width W, d ~ W ~ 10 nm. In the presence of a tunneling term t between 
the two wells, the electronic subband is split into a symmetric and an anti- 
symmetric combination, separated by the energy scale As^s = 2i. (b) Wide 
quantum well. In a wide quantum well the energy gap between the occupied 
lowest electronic subband and the unoccupied first excited subband, A S 6, is 
decreaased as compared to a narrow quantum well. 



by molecular-beam epitaxy, consist of two quantum wells spatially separated by 
an insulating barrier that is on the same order of magnitude as the width of each 
of the wells. Formally, each of the wells (layers) may be described in terms of an 
SU(2) pseudo-spin, a =t for an electron in the left well and a =[ for one in the 
right well. In contrast to the true electron spin, the Coulomb interaction does 
not respect this SU(2) symmetry - indeed, the repulsion is stronger between 
particles within the same layer (i.e. with the same pseudo-spin orientation) 
than between particles in different layers (with opposite pseudo-spin orienta- 
tion) because, in the second case, electrons may not be brought together closer 
than the distance d between the layers. In order to minimise the interaction 
energy, it is therefore favourable to charge both layers equally. Alternatively, 
this may be viewed as some capacitive energy, if one interprets the two-layer sys- 
tem in terms of a capacitor, that favours an equal charge distribution between 
the two layers as compared to a charging of only one layer. Notice, further- 
more, that tunneling, with the tunneling energy t, between the two quantum 
wells lifts the pseudo-spin degeneracy: whereas the symmetric superposition 
|+) = (| T) + I i))/\/2 of the layer pseudo-spin lowers the energy, the anti- 
symmetric superposition | — ) = (| f) — | l))/\/2 describes anti-binding. The 
energy separation between the associated subbands is given by As^s = 2i [sec 
Fig. I5.1f a)]. but it may be strongly reduced experimentally with the help of a 
high potential barrier separating the two wells. The term As^s, which plays 
the role of a Zeeman gap (though in the x-quantisation axis), may become the 
lowest energy scale in the system, such that the SU(2) pseudo-spin symmetry 
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breaking only stems from the difference in the Coulomb interaction between 
particles in the same and in different layers. 

5.1.4 Wide quantum wells 

Another quantum Hall system that may be characterised as a multi-component 
system is a wide quantum well [Fig. 15.1T b)]. Indeed, the samples which reveal 
the highest mobilities are those fabricated in wide quantum wells, where the 
well width w is often much larger than the magnetic length Ib- As compared 
to a narrow quantum well, the energy difference between the lowest and the 
first excited electronic subbands, which are the energy levels of the confinement 
potential in the z-direction, is strongly decreased. Although the Fermi level still 
resides in the lowest electronic subband (pseudo-spin a =T), the energy gap to 
the next unoccupied one (pseudo-spin a =\) may then become smaller than 
the relevant Coulomb energy scale. In the same manner as for the electronic 
spin, one must therefore no longer discard higher electronic subbands. In a first 
approximation one may restrict the calculations to these two lowest subbands 
[85l |8T] although the next higher subbands also shift to lower energies and need 
eventually be taken into account. Similarly to the quantum-Hall bilayer, which is 
sometimes also used in the description of the large quantum well, the Coulomb 
interaction decomposed in these electronic subband states is not pseudo-spin 
SU(2)-symmetric. 

In the remainder of this chapter, we discuss some aspects of correlated states 
that one encounters in multi-component quantum Hall systems in general, start- 
ing (Sec. I5.2[) with the completely spin-polarised state at v = 1 (quantum Hall 
ferromagnet) and its various manifestations in the different quantum Hall sys- 
tems described above. We will not discuss, for reasons of space limitation, the 
amazing physical properties of the elementary excitations of the quantum Hall 
ferromagnet, which is a topological spin-texture state (skyrmion), and refer the 
interested reader to the literature [86l [87J El E] • In the line of the preceding 
chapter, we have chosen to discuss a generalisation of Laughlin's wave function, 
which we owe to Halperin [84] . in order to account for the electronic spin (Sec. 
15.31) . These wave functions are further generalised to even more components 
than two, and we close this section with a discussion of their possible use in the 
description of multi-component FQHE states. 

5.2 The State at v = 1 

If one takes into account internal degrees of freedom, the state at v = 1 is no 
longer simply a Slater determinant of all occupied quantum states in the lowest 
LL, but one must take into account the macroscopic degeneracy due to the fact 
that each state \n,m) may now be occupied by 0, 1 or 2 particles. In this 
sense the situation at v = 1 is much more similar to the FQHE in a partially 
filled LL than to the IQHE which one obtains for completely filled LL [86] , and 
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the macroscopic degeneracy is again lifted by the mutual Coulomb interactions 
between the electrons. 

5.2.1 Quantum Hall ferromagnetism 

We first consider the generic case of electrons at v = 1 in the conventional 
monolayer quantum Hall system while taking into account their physical spin. 
In view of the above-mentioned energy arguments, we completely neglect the 
Zeeman effect, which would otherwise trivially lift the macroscopic degeneracy 
at v = 1 by polarising all electron spins. Because of the fact that two electrons, 
with opposite spin, may now occupy the same quantum state \n, m), the electron 
pair may in principle be in a relative angular momentum state with m = - 
the Pauli principle, which only applies to fermions of the same species, does no 
longer prevent this quantum number to be odd. Indeed, such an electron pair is 
described by a two-particle wave function with the rather unspectacular poly- 
nomial factor (z,,! — = 1, where z^f is the position of an arbitrarily chosen 
spin-f electron and Zj^i that of a spin- 1 electron. Such an electron pair therefore 
interacts via the Haldane pseudopotential vq , which is the largest pseudopoten- 
tial in the case of a repulsive Coulomb interaction because it characterises the 
interaction at the shortest possible length scale (see Fig. !4.2p FI Since vq ~ 2vi, 
the system thus tends to avoid double occupancy, and the ground state is de- 
scribed by the fully anti-symmetric (orbital) wave function (14. 7p regardless of 
whether the electron at the position Zj is spin-f or spin-|. 

Notice that, although both spinless and spin- 1/2 electrons are described by 
the same wave function, the physical origin of these ground states is different: 
in the case of spinless fermions, it is simply the non-degenerate wave function 
described by a Slater determinant, whereas in the case of electrons with spin, 
the state is formed in order to minimise the mutual Coulomb repulsion. 

Because the orbital wave function (|4.7p for electrons with spin at v = 1 is 
fully anti-symmetric, the spin wave function describing the internal degrees of 
freedom must be fully symmetric, e.g. 

Xfm = | h,h, ---An), (5.3) 

in order to form an overall wave function that is anti-symmetric. The subscript 
indicates the index of the particle that the spin is associated with. The global 
wave function, therefore, reads 

Vv=i,fm = (zk ~ zi) ® | Ti, T2, Tiv)- (5.4) 

k<l 

This is nothing other than a (spin) wave function of a quantum fcrromagnct, 
similar to ferromagnetism in a usual Fermi liquid. Indeed, the spontaneous spin 
polarisation in a Fermi liquid is also due to a minimisation of the Coulomb 

3 This pseudopotential, as well as any other with an even value of m, does not play any 
physical role due to the Pauli principle if one considers only spinless electrons, as we have 
mentioned in Sec. 14.2.21 
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repulsion by the formation of an anti-symmetric orbital wave function. Notice, 
however, that the spin polarisation in a Fermi liquid comes along with an energy 
cost as a consequence of the mismatch between the Fermi energies of spin-f and 
spin- J, electrons. The competition between the gain in interaction energy and 
the cost in kinetic energy determines the final polarisation of the system, which 
is never perfect. In the case of the quantum Hall fcrromagnct, there is no cost 
in kinetic energy when the system is fully polarised because all quantum states 
have the same kinetic energy, and the system is therefore fully polarised. 

Collective excitations 

Because the spontaneous spin polarisation in the quantum Hall fcrromagnet 
chooses, in the absence of a Zeeman effect, an arbitrary direction in the three- 
dimensional spin space, one is confronted with a spontaneous SU(2) symmetry 
breaking. As a consequence of this broken continuous symmetry, there exists a 
gapless collective excitation (Goldstone mode) the energy of which tends to zero 
in the long wave-length limit. Indeed, even if we have chosen the ferromagnct 
in Eq. (|5.3[) to be oriented in the z-direction, any other orientation, such as the 
one described by the wave function 

JY 

I ii.la, -An) or ®|+,-) = |+x,+2, -,+iv), 

i=i 

where the +j sign indicates the symmetric superposition \+j) = (| Tj) + I lj 
))/v / 2ofboth spin orientations of the j-th electron, would also describe a ground 
state. The Goldstone mode in the large wave-length limit may then be viewed 
as a global rotation of all spins into another ground-state configuration, which 
naturally does not imply an energy cost. 

In the case of a ferromagnet, the Goldstone mode is nothing other than the 
spin-density wavcQ that disperses as lu oc q 2 in the small wave vector limit, 
qls <C 1- At first sight, this mode seems in contradiction with the observation 
of a quantum Hall effect at v = 1, even in the absence of a Zeeman effect, which 
requires a gap as we have seen above. Notice, however, that this gap needs to 
be a transport gap in which a quasi-particle moves independently from a quasi- 
hole in order to transport a current. This is not the case in a spin wave with 
qls <C 1, but one obtains freely moving quasi-particles and quasi-holes in the 
limit qls 3> 1. In this limit, the spin- wave dispersion tends to a finite value that 
is given by the exchange energy between particles of different spin orientation 
and that is proportional to the interaction energy scale e 2 /els, as in the case of 
the FQHE [87]. 

There are more exotic spin-texture excitations (skyrmions), which are de- 
scribed by a topological quantum number associated with the winding of the 
spin-texture. These are gapped excitation which carry an electric charge related 



4 Remember that for a crystaline ground state (WC), the Goldstone mode is the acoustic 
phonon, as we have briefly discussed in the previous chapter in Sec. 14.11 
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to this topological quantum number. As mentioned above, a detailed discussion 
of these amazing excitations is beyond the scope of the present lecture notes. 

5.2.2 Exciton condensate in bilayer systems 

The v = 1 in a bilayer system is remarkably different from the quantum Hall 
fcrromagnct described in the preceding subsection. Although the electronic 
interactions still favour a fully anti-symmetric orbital wave function (|4.7p and 
thus a symmetric, i.e. ferromagnetic, pseudo-spin wave function, the interaction 
potential is no longer SU(2) symmetric in the pseudo-spin degree of freedom^ 
As we have already mentioned above, a charge imbalance Q between the two 
quantum wells (layers) is penalised by a charging energy, Ec = Q 2 /2C, in terms 
of the capacitance C = eA/d, where A is the area of the 2D system. Because 
Q = —evn e iA = —evnsA = —evA/2irlg when all electrons reside in a single 
layer and Q = if they are equally distributed between the two layers, one 
obtains an energy cost 

E c f_d_ 

N e i d B Ib ' 

per particle in the charge-imbalcnccd state, in agreement with a more sophisti- 
cated microscopic calculation |87j . In terms of the pseudo-spin magnetisation, 
this means that in the ground-state configuration, with a homogeneous charge 
distribution over both layers, all pseudo-spins are oriented in the xy-plane. Re- 
member that a pseudo-spin | corresponds to an electron in the upper layer and 
I to one in the lower layer, and a configuration as the one described in Eq. (|5.3|) 
is therefore excluded, whereas the symmetric and anti-symmetric combinations 

N N 

X+ = (g)l+ J ) and X- = ®|-j>, 

j'=i j=i 

with \±j) = (| tj) ± | lj))/y/2 is not. These two states, which correspond to a 
fcrromagnet in the x- and the y-direction, respectively, may be generalised by 
choosing any other direction described by the angle (j) in the xy-plane, 

N 

X* = (8)|&>, (5-5) 

j'=i 

where \<pj) = [| T) + ex P(^)l D1/V2- The states \+ an d X- are obtained for 
= and <j) — 7T (modulo 2n), respectively 

Contrary to the case of the spin ferromagnet with full SU(2) symmetry, 
where a general state would be described in terms of two angles 9 and (j), the 
different possible easy-plane pseudo-spin ferromagnetic are characterised by the 
angle which may vary between and 27r. The low-energy degrees of freedom 

5 Naturally, such an anti-symmetric orbital wave function is only physical if the layer 
separation d is not too large (as compared to the magnetic length) - otherwise one would 
simply have completely decoupled layers. 
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Figure 5.2: Hall resistance measurement used to detect excitonic condensation, 
adopted from (Eisenstein and MacDonald, 2004). (a) Counterflow configuration, 
in which one drives a current I] = I through the upper layer that is flowing in 
the opposite direction as that, Ij_ = —I in the lower layer. The hole component 
of the excitonic quantum state in one layer thus moves in the same direction as 
the electron component in the other one. (b) The two curves schematically rep- 
resent, when taking into account only excitonic superfluidity, the Hall resistance 
in both layers within the counterflow configuration. Because of the relative sign 
between the currents in the two layers, the measured Hall resistances are of 
opposite sign. Electrons with no intcrlayer correlations yield the usual linear 
-B-field dependence of the Hall resistance in order to compensated the Lorentz 
force acting on them individually. In the case of exciton condensation (around 
B = 5 T), charge tranport is due to a uniform current of charge-neutral excitons, 
which are not affected by the Lorentz force, and the Hall resistance vanishes, 
as it has been observed in the experiments (Kellogg ct ai, 2004; Tutuc et al., 
2004). 
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are therefore described by a different universality class that turns out to be the 
same as the one that describes superfluidity or superconductivity. The relation 
between superfluidity and the easy-plane pseudo-spin ferromagnet in bilayer 
systems at v = 1 may indeed be understood in the following manner: on the 
average, the average filling factor per layer is = = 1/2 in order to minimise 
the charging energy due to the capacitive term, i.e. there are as many electrons 
as holes in the LLL of each layer. Naturally, because of the Coulomb interaction 
between the particles in the two different layers, an electron in one layer wants 
to be bound to a hole in the other one. Since the number of electrons in each 
layer equals, on the average, that of holes in the other one, all particles find 
their appropriate partner in the opposite layer. The electron-hole pair in the 
two layers may be viewed as a charge-neutral interlayer exciton that satisfies 
bosonic statistics [Fig. I5.2f a)] . Below a certain temperature, these bosons 
condense into a collective state that is nothing other than the exciton superRuid 
[88] 189] 190] [87] . The phase coherence between the different excitons is precisely 
described by the angle 0. 

The first experimental indication of excitonic superfluidity in bilayer quan- 
tum Hall systems was a zero-bias anomaly in tunneling experiments |91j . Indeed, 
if one injects a charge in a tunneling experiment into one of the layers and collects 
it in a contact at the other layer, the tunneling conductance dI z /dV is expected 
to be weak in the case of uncorrelated electrons because of the Coulomb repul- 
sion between electrons in the opposite layers. However, below a critical value 
of d/ls, where one expects the interlayer correlations to be sufficiently strong 
to form a phase-coherent excitonic condensate, the injected electron systemat- 
ically finds a hole in the other layer, such that tunneling between the layers 
is strongly enhanced. This strong enhancement, which due to its reminiscence 
with the Joscphson effect in superconductors [66] is also called quasi-Josephson 
effect^ has indeed been observed experimentally [9T] , 

Another strong indication for excitons in bilayer quantum Hall systems stems 
from transport measurements in the counterflow configuration, where the cur- 
rent in the upper layer 1-^=1 flows in the opposite direction as compared to 
that in the lower layer Ij = —I [see Fig. I5.2f a)]. From a technical point of 
view, it is indeed possible to contact the two layers separately such that one 
may measure the Hall resistance (and also the longitudinal resistance) in both 
layers independently In the case of exciton condensation, the charges involved 
in transport are zero because the excitons are charge-neutral objects, which are 
not coupled to the magnetic field and thus not affected by the Lorentz force. 
In addition to a vanishing longitudinal resistance, one would therefore expect a 
vanishing Hall resistance because no density gradient between opposite edges is 
built up to compensate the Lorentz force [89] [90]. This is schematically shown in 
Fig. 15.2( b). The simultaneous vanishing of the Hall and longitudinal resistances 
was indeed observed in 2004 by two different experimental groups [92] [93] . 



e Contrary to the Josephson effect, only the tunneling conductance dI z /dV is strongly 
enhanced whereas the tunneling current remains zero in the quasi-Josephson effect in bilayer 
systems. 
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5.2.3 SU(4) ferromagnetism in graphene 

The arguments in favour of a quantum Hall ferromagnetism may easily be gener- 
alised to the case of graphene, where the Coulomb interaction respects to great 
accuracy the four-fold spin- valley degeneracy, as we have described above. In or- 
der to avoid confusion about the filling factor, one first needs to remember that 
the filling factor vq in graphene is defined with respect to the charge-neutral 
point, which happens to be in the centre of the central n = LL (see Sec. 13. 5|) . 
Two of the four (degenerate) spin- valley branches are therefore completely filled 
at vq = 0, which in non-relativistic quantum Hall systems would correspond 
rather to a filling factor v = 2. Similarly the filling factor v = 1 would cor- 
respond to a graphene filling factor vq = — 1, whereas vq = 1 implies three 
completely filled spin- valley branches (v = 3). 

Let us first consider the filling factor vq = — 1 and see how the above consid- 
erations apply to graphene with its SU(4) symmetry^ In the same manner as 
for the spin quantum Hall fcrromagnct at v — 1 , the short-range component vq 
of the Coulomb potential is screened in the completely anti-symmetric orbital 
wave function (|4.7j) . and the spin part of the wave function must therefore be 
completely symmetric. Notice, however, that one may now distribute the elec- 
tron over the four internal states \m; K, f), \m; K, j), \m; K' , f) and \m; K\ J,}. 
The general spin wave function is therefore a superposition of all these states 

N 

Xsu(4) = (^) (u m ,i\m;K,1) + u m ^\m;K,\) + u„ lt3 \m; K', j) + u m , 4 |m; K', J.)) , 

rn—1 

t (5-6) _ 

where the complex coefficients u m< i satisfy the normalisation condition X)i=i l u m,i| 2 = 
1. In the case of global coherence, all coefficients are independent of the guiding- 
centre quantum number to, u m ^ = Ui, and one thus obtains the spin wave func- 
tion of an SU(4) ferromagnetism [94, 95, 96, 97 . These arguments may also be 
generalised to the case of vq =0, where two branches are completely filled [97] , 
but the ground state does not reveal the same degeneracy as the SU(4) ferro- 
magnct at vq = ±1. Indeed, a general argument on l^-componcnt quantum 
Hall system shows that one has generalised ferromagnetic states at all integer 
values of the filling factor v = 1, K — 1 [98] . 

As a consequence of the SU(4) quantum Hall fcrromagnct, one may expect 
a quantum Hall effect in graphene at the unusual filling factors vq = 0, ±1. 
Remember that these states do not belong to the series p.22[) . vq = ±2, ±4, ... 
of the RQHE which may be explained by LL quantisation within the picture 
of non-interacting relativistic particles. In the same manner as for the spin 
quantum Hall fcrromagnct, the gapless spin-density-wave modes, which reveal 
a higher degeneracy due to the larger SU(4) symmetry, do not imply that the 
charged modes are also gapless. Indeed, the elementary charged excitations of 
the SU(4) quantum Hall fcrromagnct are generalised skyrmions [97l [99] which 
are separated by a gap from the ground state, which therefore describes an 

7 The filling factor vq = 1 is related to vq = — 1 by particle- hole symmetry and therefore 
does not require a separate discussion. 



Multi-Component Wave Functions 



99 



incompressible quantum liquid that displays the quantum Hall effect. A quan- 
tum Hall effect has indeed been observed at these unusual filling factors [22] , 
in agreement with the formation of an SU(4) quantum Hall ferromagnet. How- 
ever, there exist alternative scenarios to describe the appearance of a quantum 
Hall effect at these filling factors [10011101] 1102] and a clear indication of SU(4) 
quantum Hall fcrromagnctism is yet lacking. 

Wc finally emphasise that an SU(4) description is not restricted to graphcnc. 
Indeed, if one takes into account the electron spin, the bilayer quantum Hall 
system and its excitations may also be treated within the SU(4) framework 
[9"8" ]ll03| l3"]l9"9] although the interaction does not respect the full SU(4) symmetry 
because of the asymmetry in the layer pseudo-spin described above. 

5.3 Multi-Component Wave Functions 

Until now, we have considered a multi-component quantum Hall effect at the 
integer filling factor v = 1 (or other integer fillings in the case of graphcnc) 
that is described in terms of the Vandermonde determinant (|4.7[) Ofe<z( z fc ~ z i) 
regardless of whether the particle at the position is in a state | or J.. The 
spin orientation has only been taken into account within a spin wave function 
that is multiplied to the Vandermonde determinant. One may naturally ask the 
question whether one may also describe other filling factors than v = 1. 

A simple generalisation of the quantum Hall fcrromagnctism to other filling 
factors consists of replacing the Vandermonde determinant by, e.g., the Laugh- 
lin ([473]) at v = l/(2s + 1) or the Jain wave function (|4~30]) at v = p/(2sp + 1) 
and to multiply it again with a spin wave function that is naturally ferromag- 
netic because the orbital wave function remains anti-symmetric. There are, 
however, more general states for which the orbital wave function is not fully 
anti-symmetric, but only in the intra-component parts as it is required by the 
Pauli principle. These states are described in terms of wave functions proposed 
by Halperin in 1983 [84] that we present in this section, as well as a natural 
generalisation to systems with more components than K = 2. 

5.3.1 Halperin's wave function 

Halperin's wave function for spin- 1/2 electrons is a straight-forward generali- 
sation of Laughlin's proposal (|4.3|) . We consider the particle positions to be 
separated into two sets {z\, z\, Zpf } for spin-| particles and z\, z^ N } 
for spin-j particles. If the particles with different spin orientation could be 
treated as independent from one another, i.e. in the absence of an interaction 
between spin-| and spin-j particles, one would simply write down a product 
ansatz 

^, mi ({z}}) x # )roa ({*j}) = n ( z i - z i) mi n (** - **T 2 (57) 

k<l k<l 
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of two independent Laughlin wave functions that need not necessarily be de- 
scribed by the same exponent m. The total filling factor would then be simply 
the sum u = i/f + of the filling factors = l/m\ and = \jmi for spin-f 
and spin- J. particles, respectively. 

Apart from the fact that this situation is not particularly interesting, it is 
also unphysical because the Coulomb interaction does not depend on the spin 
orientation of the particle pairs. In the wave function (|5.7p . two particles of 
opposite spin orientation may be at the same position, i.e. the wave function 
does not vanish in general for z\ — z\ . Remember that such a double occupancy 
of the same position would be penalised by an energy cost on the order of the 
short-range component vq in a pseudopotcntial expansion. 

In order to account for these inter-component correlations, Halperin pro- 
posed to add a factor Ilfcli Ilila^fc ~ z \) n to the wav e function (|5.7p the 
exponent of which does not necessarily need to be odd because particles of op- 
posite spin orientation are not constrained by the Pauli principle. Halpcrin's 
wave function 



JV T N: Nr N: 

• mi 



<,*,,„({*;, 4» = n {4 - 4) mi n H - 4) 2 n n ( 



J _ I 
z k ~Z 



k<l k<l k=l 1=1 



(5.8) 

is therefore characterised by the set (mi,m,2,n) of three exponents. 

In analogy with Laughlin's wave function, for which we have v = 1/m, the 
exponents fix the (component) filling factors, as one may see from the power- 
counting argument (see Sec. 14. 2p . According to this argument, the maximal 
exponent for a particular particle position cannot exceed the number of flux 
quanta Nb threading the area A of the 2D electron system. Apart from the 
shift that vanishes anyway in the thermodynamic limit, one obtains the two 
equations 

Nb = m 1 N- ] + nN l and N B = m 2 N i + nN v (5.9) 

This means that, contrary to the simpler case of Laughlin's wave function, the 
number of zeros in one component is not simply given by the corresponding 
exponent times the number of particles in this component (first term in the 
above expressions). Instead, it is also affected by the particles in the other 
component that contribute each a zero of order n (second term) due to the 
mixed term in Halperin's wave function (|5.8|) . In terms of the component filling 
factors, 

u. = £, (5.10) 
Eq. (|5.9p may be rewritten in matrix form 

1 \ _ / mi n \ / 
1 J \ n m% J \ vi 

from which one obtains the component filling factors by matrix inversion 

v-\ \ _ 1 / m 2 -n \ / 1 
fj. / miTO2 - n 2 \ —n m\ J \ 1 



(5.11) 



(5.12) 
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and one finds 

m 1 +m 2 -2n 
v = vt + v I — ^— (5. 13) 

m\m 2 — n 

for the total filling factor. 

One first notices that, in Eq. (|5.12|) . not only the filling factors are fixed 
by the exponents but also, for a given magnetic field (i.e. a given number of 
flux quanta), the number of particles per component. Contrary to what one 
could have expected from the expression of Halperin's wave function (|5.8[) . the 
numbers N a , namely the ratio between them, cannot be chosen arbitrarily. 

Furthermore, the above expressions (|5.12[) and (|5 . 13|) for the filling factors 
are ill-defined if the exponent matrix in Eq. (|5 . is not invertible, i.e. when 
its determinant is zero, mim 2 — n 2 = 0. The only physically relevant situation 
arises when all exponents are equal odd integers mi = m 2 = n. However, this 
result should not surprise us: we are then confronted again with a completely 
anti-symmetric wave function, actually a Laughlin wave function, which requires 
a ferromagnetic spin wave function. As we have seen above, in the discussion of 
the quantum Hall ferromagnetism, the ground-state manifold comprises states 
with different polarisation along the z-axis: the state with iVj = N and N± = 
is an equally valid ground state as a state with TV-j- = iVj = N/2 or iVf = and 
iVj = TV, where N = TVj + iVj is the total number of particles. The component 
filling factor is therefore not well-defined and depends on the polarisation 

TVt - N i — v i 

Vz = AT = ~, 5.14 

N v 

whereas the total filling factor is simply given by v = 1/m, in terms of the com- 
mon odd exponent m. Notice that contrary to the quantum Hall fcrromagnet, 
a state with an invertible exponent matrix has a polarisation that is completely 
fixed, 

Pz = ■ (5.15) 

TOl?7l2 — n 

Wc finally mention that not all states that can be written down in terms of 
Halperin's wave function are good candidates for the description of the ground 
state chosen by the system. One may show, e.g. within a generalisation of 
Laughlin's plasma analogy (presented in Sec. I4.2.5[) to two or more compo- 
nents, that several of Halperin's wave functions do not describe a homogeneous 
liquid but a liquid in which the different components phase-separate [104) . For 
two components, the condition for a homogeneous state is simply that both 
the exponents m\ and m 2 , which describe the intra-component correlations, 
must be larger than n for the inter-component correlations. As an example, 
we may study the states (3,3, 1) and (1, 1,3), which would both be candidates 
for a possible two-component FQHE at v — 1/2 and which have indeed been 
investigated in the literature [105] . However, only the first one describes a ho- 
mogeneous liquid, such that the second one may be discarded right from the 
beginning. 

Furthermore, some of Halperin's wave functions, even if they satisfy the 
above-mentioned condition, turn out to be problematic if the interaction is 
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SU(2) symmetric, such as for the true electron spin. In this case, one may 
show that (m, m, n) states are only eigenstates of the total-spin operator, which 
commutes with the interaction Hamiltonian, if n = m (i.e. in the ferromag- 
netic state) or if n = m — 1 1 . This restriction may be omitted though in 
bilayer quantum Hall systems or in wide quantum wells where the interaction 
Hamiltonian is not pseudo-spin SU(2)-symmetric. 

Physical relevance of Halperin states 

A physically relevant Halperin state is e.g. the unpolarised (3, 3, 2) state which 
would occur at a filling factor v = 2/5. Remember from the discussion of CF 
theory in Sec. 14.41 that there is also a (naturally polarised) CF candidate, with 
p = 2 completely filled CF LLs, to describe the ground state at this filling fac- 
tor. Which of them is now the better one? This question could be answered 
within exact-diagonalisation calculations, which showed that, in the absence of 
a Zceman effect, the true ground state is described in terms of the unpolarised 
Halperin wave function (3,3,2) 106 . Notice, however, that the energy differ- 
ence between the two states is extremely small, as may be seen from variational 
calculations [76], such that the polarised CF state becomes the ground state 
above a critical value of the energy Az associated with the Zeeman effect. This 
critical value would therefore describe a phase transition between an unpolarised 
and a fully polarised FQHE state. Such transitions have indeed been observed 
in polarisation experiments, where the strength of the Zceman effect was varied 
by a simultaneous change in the magnetic field and in the electronic density 
p7l[T08] . 

5.3.2 Generalised Halperin wave functions 

We would finally mention that Halperin's wave function may easily be gener- 
alised to describe possible FQHE states in systems with a larger number of 
components, such as the four spin-valley components in graphene. This gener- 
alised wave function for l^-componcnt quantum Hall systems may be written 
as a product 



({*M , .->*£ ) }) =^,...,™ K x cr (5.i6) 



of a product of Laughlin wave functions 




for each of the components and a term 
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that takes into account the correlations between particles in different compo- 
nents [109 1 . Here, the indices i and j denote the component, i,j = 1, K, and 
z^} is the complex position of the fc^-th particle in the component i. 

Although the wave function (|5. 16|) may seem scary at the first sight, it is as 
easily manipulated as Halperin's original wave function (|5.8|) . The component 
filling factors fj = Nj/Ns may be determined, in the same manner as in the 
two-component case (|5.11[) , with the help of the "exponent matrix" M. the off- 
diagonal terms of which arc the exponents {M.)ij — n-ij (for i ^ j), whereas the 
diagonal terms are simply the exponents corresponding to the intra-component 
correlations, (M)u = m*. The zero-counting argument yields the matrix equa- 
tion 



relating the component filling factors to the exponents, and if M is invertible, 
all component filling factors are fixed by the inverse equation 



If the determinant det(AI) is zero and the matrix thus not invertible, not 
all component filling factors can be determined. In analogy with the two- 
component case this hints at underlying ferromagnetic states. A perfect SU(i£") 
ferromagnetic state is obtained when all components are equal odd integers, 
mi = riij = 77i, in which case one obtains again a simple (fully anti-symmetric) 
Laughlin wave function for all particles regardless of to which component they 
belong. For K = 4 and m = 1, this is just the SU(4) ferromagnetic state at 
v = 1 which we have already discussed in the context of the quantum Hall effect 
at vq = ±1 in graphene (Sec. 15. 2. 3D . 

Notice, however, that contrary to a two-component system, where one only 
needs to distinguish between an invertible and a non-invertible matrix, the sit- 
uation is much richer for K > 2. One may indeed have different "degrees" of 
invertibility that are described by the rank of the matrix. Consider, e.g., the 
fully anti-symmetric wave function with m t = n,-j = m. In this case, Eq. (|5.17[) 
actually consists only of one single equation relating the component filling fac- 
tors, i.e. 1 = m(y\ + ... + vk) = vnv, and all other lines of the matrix equation 
are simply copies of the first one. The rank of this matrix is 1, i.e. only the 
total filling factor is fixed, v = 1/m [SU(if) ferromagnet] whereas in the case 
of an invertible matrix the rank is K and the K lines in the matrix equation 
(|5.17p represent (linearly) independent equations. If the rank of an exponent 
matrix is smaller than K but larger than 1, the resulting state is neither a full 
SU(-ftT) ferromagnet nor a state with completely fixed component filling factors 
(or polarisations) - it is rather a state with some intermediate ferromagnetic 
properties. 




(5.17) 




(5.18) 
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As for two-component Halperin wave functions (|5.8[) . a generalisation of 
Laughlin's plasma analogy allows one to distinguish between physical (i.e. ho- 
mogeneous) and unphysical states (which show a phase separation of at least 
some of the components). Indeed, the exponent matrix M. must have only 
positive eigenvalues in order to describe a homogeneous state |104j . We finally 
mention that A4 encodes not only information concerning the filling factors 
(|5.18|) . but fully describes the quantum Hall state (|5.16|) . such as its topological 
degeneracy, the charges of its quasi-particlc excitations as well as the statistical 
properties of the latter |110j . 



Appendix A 



Electronic Band Structure 
of Graphene 

In this appendix, wc calculate the band structure of graphene in the tight- 
binding model the results of which we have summarised in Sec. 11.2.31 
Because graphene's honeycomb lattice consists of two distinct sublattices A and 
B, the electronic wave function 

^k(r) - a k ^ A) (r) +6 k ^ B) (r), (A.l) 

is a superposition of two wave functions, for the A and B sublatticc, respectively, 
where a k and 6k are complex functions of the quasi- momentum k. Both (r) 
and ip^ (r) are Bloch functions with 

1% } (r) = £ e lk - R '^)(r + S 3 - R), (A.2) 
Ri 

in terms of the atomic wave functions phi^'(r + 5j — R;) centred around the 
position R; — S j , where dj is the vector which connects the sites R; of the under- 
lying Bravais lattice with the site of the j atom within the unit cell. Typically 
one chooses the sites of one of the sublattices, e.g. the A sublatticc, to coincide 
with the sites of the Bravais lattice such that 8a = 0. 

With the help of these wavefunctions, we may now search the solutions of 
the Schrodingcr equation 

where H is the full Hamiltonian for electrons on a lattice, which is of the type 
(|2.2| mentioned in Sec. 12.11 Here, we have chosen an arbitrary representation, 
which is not necessarily that in real space0 Multiplication of the Schrodingcr 
equation by tp^ from the left yields the equation ip^Hip^ = ek^ipk, which may 

1 The wavefunction Vk( r ) i s i thus, the real space representation of the Hilbcrt vector -!/> k . 
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Figure A.l: Tight-binding model for the honeycomb lattice. 



be rewritten in matrix form with the help of Eqs. (|A.1[) and (IA.2I) 

(4, &» Hk ( ^ ) = e k (at, 61) (A.3) 
Here, the Hamiltonian matrix is defined as 




Hk= \ ^r B> Htp^ ] ,(b)* tt ,(b) |-W k , (A.4) 



and the overlap matrix 



accounts for the non-orthogonality of the trial wavefunctions. The eigenvalues 
e k of the Schrodingcr equation yield the energy bands, and they may be obtained 
from the secular equation 

det [H k - e£S k ] = 0, (A.6) 

which needs to be satisfied for a non-zero solution of the wavefunctions, i.e. for 
a k 7^ and fo k ^ 0. The label A denotes the energy bands, and it is clear that 
there are as many energy bands as solutions of the secular equation (|A.6|) . i.e. 
two bands for the case of two atoms per unit cell. 

From now on, we neglect the overlap of wave functions on neighbouring sites, 
such that the overlap matrix (|A.5[) simply becomes the one matrix 1 times the 
number of particles N due to the normalisation of the wave functions. The 
secular equation then tells us that the energy bands are just the eigenvalues of 
the Hamiltonian matrix (|A.4[) . Furthermore, one notices that because the two 
sublattices are equivalent from a chemical point of view, we have Htjj^ = 
ip^* Htp^\ and the diagonal terms therefore contribute just a constant shift 
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to the band energies that we may set to zero. The only relevant terms are then 

l£ B ee 4{ A) *H4{ B) = Nt? 



the off-diagonal terms in Eq. lfO|) . H£ B = tp^Hip^ = Nt£ B , with the 



hopping term 

t£ B ee ^ e lkR < / d 2 r <^>* (r - R k )H^ (r + <5 AB - R m ) , (A.7) 
Ri J 

where Sab is a vector that connects an A site to a B site. 

In order to obtain the basic band structure of graphcne, it is sufficient to 
consider a hopping only between nearest-neighbouring sites described by the 
hopping amplitude 

t = J d 2 r4, A *(v)Hcf> B (v + S 3 ), (A.8) 

where we have chosen Sab = ^3 (sec Fig. IA.1[) . Notice that one may also take 
into account hopping to sites that arc further away such as the next-nearest 
neighbours which turn out to be on the same sublattice and which would thus 
yield diagonal terms to the Hamiltonian matrix. However, whereas we have 
t ~ 3 eV, the hopping amplitude for next-nearest-neighbour hopping is roughly 
10 times smaller [5T| and only marginally affects the low-energy properties of 
electrons in graphene. 

If we now consider an arbitrary site A on the A sublattice (Fig. IA.1|) . we 
may see that the hopping term (|A.7j) consist of three terms corresponding to 
the nearest neighbours B\, B2, and B3, all of which have the same hopping 
amplitude t. However, only the site -B3 is described by the same lattice vector 
(shifted by #3) as the site A and thus yields a zero phase to the hopping matrix. 
The sites B\ and B2 correspond to lattice vectors shifted by 

\/3a , r- . , VSa r- . 

a 2 = — — (e x + v3e y ) and a3 ee a 2 — ai = — — {— &x + voe y ), 

respectively, where a = | ^3 1 = 0.142 nm is the distance between nearest- 
neighbour carbon atoms. Therefore, they contribute a phase factor exp(ik • a 2 ) 
and exp(ik-a3), respectively. The hopping term (j A.7|) may therefore be written 
as 

where we have defined the sum of the nearest-neighbour phase factors 

7k= l + e' ka2 +e ,ka3 . (A.9) 

The band dispersion may now easily be obtained by solving the secular 
equation (|A.6|) . 

e A (k)=Ak£ B |=At| 7k |, (A.10) 



and is plotted in Fig. 12.21 The band dispersion is obviously particle-hole sym- 
metric, and the valence band (A = — ) touches the conduction band (A) in the 
incquivalcnt points 

47T 
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which one determines by setting 7±k = and which coincide with the two 
inequivalent BZ corners K and K' . Because the whole band structure is half- 
filled in undoped graphene, as we have mentioned in Sec. 12. li the Fermi energy 
lies exactly in these points K and K'. 



Continuum Limit 

The low-energy electronic properties may be obtained by expanding the band 
structure in the vicinity of these points, and the low-energy Hamiltonian is 
obtained simply by expanding the sum of the phase factors (|A.9[) around K and 
K>, 

7p=7k=±K+p = 1 + e ±lK ' a VP^ + e ±iK.a 3e ip.a 3 

~ 1 + e ±a */ 3 [1 + ip ■ a 2 ] + e^' z [1 + ip ■ a 3 ] 
= ±(0), ±(i) 

By definition of the Dirac points and their position at the BZ corners K and 
K', we have 7p = 7±k = 0. We limit the expansion to first order in |p|a. 
Notice that, in order to simplify the notations, we have used a system of units 
with h = 1, i.e. where the momentum has the same units as the wave vector. 
The first order term is given by 



±m -v^a 
Ip ' = « 



= T^(p x ±ip y ), (A.ll) 



which is obtained with the help of sin(±27r/3) = ±\^3/2 and cos(±27r/3) = 
— 1/2. This yields the effective low-energy Hamiltonian 

flS=£>(Px<T a + £Py<r») > (A.12) 

in terms of the Fermi velocity 

3ia 

(A.13) 

The index £ = ± denotes the valleys K and K' , and one obtains at the K point 
the Dirac Hamiltonian mentioned in (12.41) 



H D =vp-<r, (A. 14) 

whereas the low-energy Hamiltonian at the K' point reads 

H' D = up ■ a* , (A.15) 

with <r* = [a x , —a y ). Both Hamiltonians yield the same energy spectrum which 
is therefore two-fold valley-degenerate. 
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Notice that if one prefers to avoid the complex conjugation in the Hamilto- 
nian (|A.15[) , one simply changes the representation by interchanging the A and 
B sublattices, in which case one may write the Hamiltonians for the two valleys 
K (£ = +) and K' (£ = — ) in a compact form, 



(A.16) 
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Appendix B 



Landau Levels of Massive 
Dirac Particles 

Mass Confinement of Dirac Fermions at B = 

Even in the absence of a magnetic field, electronic confinement in graphene turns 
out to be quite tricky because a simple-minded approach in terms of a potential 
^conf = V(y)t cannot confine Dirac electrons. This fact is due to an intrinsi- 
cally relativistic effect that is called the Klein paradox, according to which a 
(masslcss) relativistic particle may transverse a potential barrier without being 
backscattered |112) . This effect may be understood in the following manner: 
consider an incident electron in the region with V = the energy of which is 
slightly above the Fermi energy. In the potential barrier, the Dirac point is 
shifted to a higher energy that corresponds to the barrier height and the Fermi 
energy lies now in the valence band, where the electron may still find a quantum 
state (with the same wave-vector direction and the same velocity v) - instead 
of moving as an electron in the conduction band, it thus simply moves in the 
same direction as an electron in the valence band [Fig. IB.lf a)]. This is in stark 
contrast with quantum mechanical tunneling of a non-relativistic particle, for 
which the transmission probability through a potential barrier is exponentially 
suppressed because of a lacking quantum state at the same energy as that of 
the incident electron. 

The problem is circumvented by a so-called mass confinement 

V coni = V(y)**=( V M J {y) y (B.l) 

and we discuss first the simpler case of a constant mass term Mo z that needs to 
be added to the Dirac Hamiltonian. That this term yields indeed a mass may 
be seen from the Dirac Hamiltonian at B = 

H%=v V -* + Mo*=( , M . , v(Px ~J Py) V (B.2) 
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Figure B.l: (a) Klein tunneling through a barrier. An incident electron in 
the conduction band (CB) above the Fermi energy, which is at the Dirac point 
before the barrier, transverses the barrier as en electron above the Fermi energy 
in the valence band (VB). The valence band is partially emptied because the 
Dirac point has shifted to a higher energy corresponding to the barrier height. 
(b) Mass confinement. A gap opens when the particle approaches the edge, 
which becomes a forbidden region where no quantum state can be found at the 
energy corresponding to that of the incident electron. 



the diagonalisation of which yields the energy spectrum 

e x (p) = AV« 2 |p| 2 + A/ 2 , 

which is gapped at zero momentum. This is nothing other than the dispersion 
relation of a relativistic particle^ with mass m such that M = mv 2 . Qualita- 
tively one may see from Fig. IB.lT b) why a mass confinement is more efficient 
than a potential barrier. Indeed, when the particle approaches the edge with 
M(y) / a gap opens. An electron slightly above the Dirac point may then 
only propagate in the region with M = 0, whereas at the edge its energy lies in 
the gap which is a forbidden region, and the electron is thus confined. 

Similarly to the B = case, one may find the energy spectrum of the massive 
Dirac Hamiltonian (|B.2jl in a perpendicular magnetic field, which reads, in terms 
of the ladder operators a and aJ, 



Hi 



M 

v{n x + m y ) 



v(U x - Uly) 

—M 



M V2ffc 



(B.3) 



Its eigenvalues may be obtained in the same manner as in the M 
Sec. I2.3.2|) . and one obtains 



e Xn = AWM 2 + 2 



h 2 v 2 

1 B 



case (c.f. 



(B.4) 



1 The sign A = 



— corresponds to the anti-particle. 
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for the massive relativistic LLs, n 7^ 0. 

Special care needs to be taken in the discussion of the central LL n = 0, 
which necessarily shifts away from zero energy. The associated quantum state 
(|2.24p is zero in the first component uq, whereas the second component is given 
by vq = |0). In order to satisfy the second line in the eigenvalue equation 

r*^ _ ^ I V2%a\{ \ ( 



one needs to fulfil 

at u = (eo + M)«o = (e + M)|0), (B.5) 

such that the only solution is eo = —M. The relativistic n = LL is therefore 
shifted to negative energies and does no longer satisfy particle-hole symmetry. 
This effect is called parity anomaly and depends on the sign of the mass. 

In the case of graphene, we need to remember that there are two copies 
of the energy spectrum, one at the K point and one at the K' point. As 
we have disussed in Appendix [A] the Hamiltonian (|B.3[) describes the low- 
energy properties at the K point whereas we need to interchange the A and B 
sublatticcs at the K 1 point and add a global sign in front of the off-diagonal 
terms [see Eq. (|A.16[) ]. 




Naturally, the eigenstatcs of this Hamiltonian arc the same as those of the 
Hamiltonian (|B.3j) at the K point, but the eigenvalues change their sign. Due 
to the particle- hole symmetry of the levels (|B.4[) , the global sign does not affect 
the energy spectrum for n 0. However, the n = LL, which does not respect 
particle-hole symmetry, must again be treated apart, and one finds in the same 
manner as for the K point the condition corresponding to Eq. ()B.5I) , 



- V2^ a f u = (eo - M)v & = (e - M)|0) . (B.7) 

IB 

One notices that the n = LL level at the K' point shifts to positive energies as 
a function of the mass, such that the overall level spectrum for graphene, when 
one takes into account both valleys, is again particle-hole symmetric, but the 
valley degeneracy is lifted for n = 0. 

The case of a mass term that varies in the y-dircction, such as for the mass 
confinement potential, may finally be treated in the same manner as we have dis- 
cussed in Sec. 13.1.21 the system remains translation-invariant in the ^-direction, 
such that the Landau gauge is the appropriate gauge and the wave vector k in 
this direction is a good quantum number. Because this wave vector determines 
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the position of the eigenstate in the y-direction, yo = kl 2 B , the energy spectrum 
is given by the expression (|3.2ip . 



h 2 v 2 



eA„, yo ;? = AW M 2 (y ) + 2-p-n, (B.S 



for n 7^ and both valleys £ = ±, whereas the n = LL is found at 

en=o,„ ;« =-^M(y ). (B.9) 
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